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The crews... the supervisory staff ... the equipment... 
and the experience ... for the successful completion 


of long or short term contracts covering every phase of 


geophysical exploration work. 


s record of more than 14 years continuous service 


independent’ 
to a tong list of important oil producers in North and South 
America merits your confidence. You are invited to consult 


with us obout your oil exploration problems. 
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236 E. Courtland Street, Philadelphia 20, Penna. 


VIRGIL KAUFFMAN HOMER JENSEN 
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Compositing Panel 


The function of this compositing 
panel is to combine energy from 
each recording channel with an 
equal portion of energy from each 
of the adjacent recording channels 
to produce an increased directional 
discrimination against interfering 
superficial waves. In order to in- 
sure an equal input of energy from 
each recording channel to the 
compositing panel, this panel is 
introduced into the circuit between 
the amplifiers and the recording 
oscillograph, the automatic gain 
control of the amplifiers insuring 
an equal input of each recording 
channel to the compositing panel. 


The compositing circuit has been 
designed so that energy other than 
that desired in the compositing is 
less than one-tenth of one per cent. 
The switching arrangement pro- 
vides for the use of any or all 
channels either composited or in- 
dependent as desired. With this 
switching arrangement a seismom- 
eter spread with independent end 
traces and center traces may be 
made up of either 24 traces or 12 
traces. 


SPECIFICATIONS 
Weight—35 pounds 


Dimensions— 
12”x2014,’x11 5/16” 
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Precision Machine Shop 


po precision machine shop is maintained by National Geo- 
physical Company for the use of the research and development 
staff. This shop is completely equipped for all types of pre- 
cision instrument work, and staffed with instrument makers. 
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25-B Power Panel 


The power panel serves as a central power amplifiers is connected to the pov 
source for all units of the equipment. It by means of a power cable whic 
supplies voltage from the storage battery, disconnected at the plug near th 
300-volt voltage from a vibrapack, and a fier. This facilitates the use of cit! 
high frequency exciting voltage for the 18, or 24 amplifiers. The 240-amp 
operation of the automatic gain control. 6-volt heavy-duty storage batter 
A voltmeter is located on the power panel connected to this power panel f 
to facilitate the checking of the various the entire power necessary for th 


power supplies. Each amplifier case of 6 tion of the recording equipment. 


SPECIFICATIONS 
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SPECIFICATIONS 


COMMUNICATION UNIT 
Number of tubes—9 


Dimensions— 
RECEIVER 

W eight—781/,, pounds Power output—1 watt 

Power input: A.V.C. range—essentially flat 
Standby—32 watts from 1 to 10,000 microvolts 
Receiving—50 watts Sensitivity—1 microvolt to 
Transmitting—150 watts give 1 watt output with 5:1 

TRANSMITTER signal to noise ratio 
Power output—10 watts Frequency—30-40 mc. Fixed 


Frequency—30-40 mc. Crystal tuned. Crystal controlled 
controlled Number of tubes—8 


Power source—Because of the necessity for conserving battery 
power where the only source of supply is a six-volt storage 
battery, the current drain has been kept to a very low value, 
especially for the receiver which operates all the time com- 
munication is being maintained. When the transmitter is to 
be operated, a push-button on the microphone brings the 
transmitter into operation. 


Antenna—A special type of vertical half-wave antenna fed by 
a flexible co-axial transmission line is used. The antenna 
gives a low angle of radiation and an effective increase in 
power. The antenna is mounted on top of the truck on a gear 
mechanism so that it may be raised and lowered without 


climbing up on the truck. 


The transmitter and receiver are crystal controlled, and neither 
requires adjustment over long periods of time after installation. 
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ATLAS EXPLOSIVES 


designed especially for 
GEOPHYSICAL PROSPECTING | 


Atlas Petrogel No. 1 HV 60%. The outstanding explosive of the 
seismic prospecting field. A gelatinous type explosive made for 
the job. Combining high strength with high velocity and water 
resistance, it handles easily in the field—will force through mud 
in drill holes without bulging the cartridges—yet is easy to prime. 
Packed by weight. 

2 Ibs. x 2” diameter 

5 lbs. x 2144” diameter 


Atlas Twistite Assembly. A fast coupling device enabling the 
blaster to connect Atlas Petrogel explosives into continuous rigid 
columns. Each cartridge of explosive is threaded into a sleeve, 
which in turn, may be threaded over an adjacent cartridge. Con- 
venient to use. 


Atlas Technical Service. Call on the nearest Atlas Office for 
service on special explosives problems, Our trained representa- 
tives will be glad to help. 


Remember—You Get the Extra Safety of Manasite Too 


ATLAS 


Powder Company 
Wilmington 99, Delaware 


Please mention GrEopHysics when answering advertisers 


5 
& 
; 
: | 
a8! 
BATT 
FE 


GEOPHYSICS the Journal of the Society of oe Geophysicists 


the Varied Needs 
Seismic Prospectors 


packaged in 3 sizes, “Nitra- 
mon” S now has wider application than ever 
before! 

“Nitramon” S* is supplied in heavy metallic, 
watertight containers. These are easily and quickly 
joined together, forming rigid charges which can 
be readily loaded under the most difficult conditions. 
To meet varied requirements in the field, the follow- 
ing container sizes are offered: 


Diameter Weight of Units 
af 1 Ib. or 4 Ibs. 
= 5 Ibs. or 10 Ibs. 
16% Ibs. 


Used with Du Pont “SSS” Blasting Caps 
and ‘‘Nitramon” S Primer, “Nitramon” S 
provides the safest, most convenient and de- 
pendable combination yet developed for 
seismic prospecting. Ask your Du Pont field 
representative for complete information. 


E. 1. DU PONT DE NEMOURS & CO. (INC.) 
EXPLOSIVES DEPARTMENT 
WILMINGTON 98, DELAWARE 


*“Nitramon" S cannot be deto- 
nated by the strongest commercial 
caps, open flame, friction, falling 
objects or impact of rifle bullets. 
Yet when used with a combination 
of“ Nitramon” S Primer and either 
electric blasting caps or Prima- 
>», cord, it is readily detonated. 


DU PONT 


“NITRAMON” S 
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ite Shot. Hole Casing 


e In every phase of geophysical shot hole work, Fedralite Plastic Shot 
Hole Casing means big savings. It's “right" for swamp or prairie loca- 
tions. And that's why it was quickly accepted by the trade . . . because 
Fedralite does provide a strong, light weight, trouble-saving, economical 
shot hole pipe in the field. 


Fedralite Plastic Shot Hole Casing is easy to handle at the shot... 

convenient to transport. It provides plenty of strength for jetting 

. : and drilling operations. It is smooth-surfaced, with no burrs or sharp 

@ Light Weight edges to injure hands. It is factory threaded with firm, clean, smooth 

@ Strong-Tough threads to promote faster couplings, tighter joints. Couplings are 
attached at the factory. 


© Sole Fedralite Plastic Shot Hole Casing has been continually improved ever 
e Factory- Threaded since it was introduced. Federal Engineers have followed it into the 
P field, watched it work, checked results, added refinements. That's why 
© More Economical today Fedralite Plastic Shot Hole Casing is the best ever made .. . why 
it provides even faster, more economical shot hole work. That's why, too, 
it's making more trips to the field every day. Ask a present user about 
Fedralite Plastic Shot Hole Casing. You can get immediate delivery 

from Gulf Coast stocks. 


FEDERAL ELECTRIC COMPANY, INc. oF Texas 


700 Waugh Drive 
405 Velasco Street 701 | Grand 
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TYPE VR (left): A newly designed recording Base 


Station Magnetometer which is self-con- 


tained and permits visual observation 
while recording time, temperature and 
diurnal variation. No darkroom or auxiliary 
equipment necessary. 

Standard sensitivity: visual 10 gammas 


per scale division; recorded 10 gammas per 


millimeter. 


INSTRUMENT 
CORPORATION 


4607 MONTROSE BLVD. HOUSTON 6, TEX. 


TYPE V (right): A new, improved Schmidt type 
magnetic field balance. 
Standard sensitivity: 10 gammas per 

scale division. Scale range: 1200 gam- 


Ruska Instruments are of the latest design and 


are based on most recent developments. They 


are superior in precision, workmanship and ma- 
terial, are easy to operate and maintain, and 


have a fine appearance and a durable finish. 
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THE PLANNING OF A FOREIGN 
GEOPHYSICAL OPERATION* 


HERBERT M. DAWSONT 


ABSTRACT 


This paper contains a discussion of the advance preparation necessary for a 
properly conducted foreign geophysical operation. Consideration of the economic fac- 
tors involved emphasizes the necessity of having proper equipment and an efficient 
supply system. 


It is the purpose of this paper to outline the necessary steps which 
should be taken before any geophysical work is started in a new and 
more or less unknown area—at least, unknown to the party or parties 
planning the work. It is presumed that the company has already taken 
the necessary legal steps for doing business in the country, and has 
obtained a concession. This discussion is intended for companies en- 
tering into foreign exploration for the first time. The experience of the 
author has been in South America only, so the particular problems en- 
countered there may be unduly stressed, and the special problems of 
some other part of the world omitted, but the general outline should 
be the same for any foreign area. This discussion will necessarily be 
only an outline, as the limited space does not allow for details. This 
paper is more concerned with seismic work than any other type of 
geophysical operation, because it is by far the most expensive but 
the general steps should apply to any type of geophysical exploration. 

As every geologist and geophysicist knows, the first step is to 
have a competent and experienced man, thoroughly familiar with 
foreign geophysical operations, make a preliminary but complete 


* Read at the Chicago and Houston meetings of the society, 1946. 
} Petty Geophysical Engineering Co., San Antonio, Tex. 
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study of the area. This is many times more important in foreign work 
than in domestic work, as there are many more problems involved. 
Naturally, the physical aspects should be studied carefully to deter- 
mine the type of terrain the party or parties will have to work, and if 
it turns out to be virgin jungle country instead of cultivated farms and 


grazing lands, then the problems are multiplied many times. Also, if 


the area is unpopulated, or very sparsely populated, then all labor and 
food has to be brought from the outside. 

The accessibility of the area to the port of entry into the country 
and the means to be used to move the equipment to the area should 
be studied, as well as the proximity of towns and villages. The com- 
plete transportation system will have to be studied, and equipment 
ordered to arrive as soon as, and usually earlier than, the first geo- 
physical party. Most of the foreign areas having good roads are either 
not in possible oil producing territories, or have already been explored. 
It may be that a river will provide the cheapest and easiest access for 
supplies, in which case boats will have to be ordered if they are not 
locally available. It may be that mule transport will be the only pos- 
sible means of entry, in which case trails and pastures for the mules 
will have to be cleared months in advance of the actual beginning of 
geophysical operations. In cases it has been found that an amphibian 
plane provides the most economical means of moving and supplying 
the party. Whatever means is used—either plane or mule—it will 
take time to procure equipment, and the necessary units must be 
ordered well in advance. 

The above indicates the necessity of determining well ahead of 
time the means to be used to transport men, equipment and supplies 
into and out of the area. The whole system of transportation should be 
set up in advance of the arrival of the first party in order to facilitate 
the beginning of the work. Supplying a party of one hundred men, 
which is approximately what the seismograph party will carry, 
is a considerable operation. The transportation facilities will have 
to be set up to handle a minimum of twenty tons per month, and 
that is many mule loads. In all of the calculations about the transpor- 
tation system it is desirable to allow a wide margin of safety when de- 
termining the amount of transportation equipment which should be 
available to the party. The operation might require more men than 
is at first anticipated, and therefore more supplies, or it might con- 
sume larger quantities of dynamite or casing, or some other thing 
may show up which could not be originally calculated into the esti- 
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mate. It is definitely not economical to try to get by on a minimum 
of transportation facilities. 

A thorough study should be made of the climatic conditions in the 
region to be worked. If no accurate data are available, as is usually 
the case in isolated areas, then just talking to the local residents 
around the country will at least be of some help; in almost every place 
they know when the rainy season or the dry season is expected, and 
they usually do not miss it too far for the planning needed for geo- 
physical operations. Trying to get a new party started to work in the 
middle of the rainy season is a very costly and aggravating procedure, 
and it is not to be recommended unless extreme urgency is required, 
and a time limit set to complete a definite operation. If it is possible 
to determine in advance the parts of the area to be worked which can 
be most advantageously worked in the dry and the wet season then 
the work should be laid out this way in advance. In some areas where 
portable work is done with mule transport the work is not held up 
unduly by the seasonal rains; however, it would still be difficult for 
a new party to start operations in the middle of the rainy season, even 
in portable work. 

The availability of local labor should also be studied, and the 
sources of labor supply laid out in advance. In some areas this is ex- 
tremely important, as it may take several weeks from the time the 
operation of actually hiring labor is started until the laborers are on 
the ground at the place where the work is to start. Of course, it will 
probably have been necessary to send some labor into the area in ad- 
vance of the actual arrival of the party, to do preparatory work, and 
’ these men may be switched to the regular work of the party. 

In almost all areas the service of a qualified doctor is beneficial, and 
in some isolated areas it is an absolute necessity. Aside from the aspect 
of the safety of the men, a doctor with the crew can, by examination, 
decide whether or not a man is able to work. This may not seem im- 
portant up here, but in cases where men get paid for some sicknesses 
and slight accidents it is amazing how many supposed accidents and 
illnesses can come up if there is not a qualified doctor in camp. 

Naturally, the equipment needed for use in the work is definitely 
tied in with the type of transportation to be provided. Of course, it 
will be decided early in the investigation whether the work will be 
portable or with trucks, or by using both methods. Sometimes it is 
definitely advantageous to switch a truck crew to portable work dur- 
ing the rainy season. In either case, sufficient equipment of all types 
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that it is calculated will be needed should be provided so that the 
work will not be delayed. Most areas where trucks can be used during 
the dry season are very dry during that time, and the use of two drills 
and four water trucks on a single seismic party is highly recommended. 
In case portable work is to be done then certainly it is very important 
to include ample portable drilling equipment. Usually it is the drilling 
which slows down any seismic operation. Also, in supplying the trucks 
for the party (in case it can be truck operated or is to be supplied by 
truck), it is cheaper to furnish extra units than to have the party shut 
down at times awaiting supplies. It is much cheaper in any foreign 
operation to furnish more equipment than it is actually calculated will 
be essential, rather than to try to save on this by furnishing the least 
equipment that it is thought the party will need. 

The choosing of the proper technical equipment to fit the needs 
of the area is very important. If the equipment can be truck trans- 
ported at all times then this problem is lessened; however, most areas 
require at least a certain amount of portable work, even along with 
truck work. In portable work the size and weight of each object to be 
carried or transported has to be considered very carefully. In some 
areas it is impossible to use equipment which has single pieces larger 
than can be carried by one man, although this is unusual. This should 
be determined in advance, and if necessary, special equipment should 
be designed to fit the particular problem of the area. It is well to have 
sufficient equipment of various types on hand, so that the actual op- 
erator and user of the equipment will have some latitude of choice to 
adapt it to the area. Of course, it should be unnecessary to mention 
that a very large supply of every spare part imaginable for all equip- - 
ment should be carried along with the party. Much of the special 
equipment of geophysical parties is not available in the country of 
operations, and it will usually take at least six weeks, and probably 
longer, to get any pieces or parts of equipment in, even by air express, 
by the time it clears the customs and all the other existing im- 
pediments. Although some companies even think it desirable to carry 
extra units of the basic equipment, in case of an accident to one of the 
vital pieces, few, if any, actually do so. However, it would be advisable, 
especially if several parties are to be operated simultaneously in an 
area. 

Along with all the technical equipment comes the camp equipment 
for the men. Local conditions should be studied in advance in this 
respect, and it should even be determined what the labor law may 
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require the company to furnish with respect to labor and camp facili- 
ties on the party. The laws in this respect in some countries are be- 
coming rather stringent, and have to be obeyed to avoid trouble. As 
for the camp equipment for the technical men, it is strongly recom- 
mended that the best and most comfortable camp equipment avail- 
able, as well as the lightest and most portable, be secured. The well- 
being and contentment of the men is at all times a vital concern to the 
company doing the work, as contented men will turn out more work 
than dissatisfied men. At the present time the trend for geophysical 
camps, even in the middle of the jungle, is going somewhat toward the 
elaborate side. Electric light plants in a camp are common, and re- 
frigerators, even far out in the jungle, are becoming more common 
all the time. Naturally the old time geologists laugh at the “bunch of 
sissies”? working out in the jungle now, but they all enjoy a cold bottle 
of beer or a cold drink when they come by, if they can find one. It 
might be mentioned that a complete list of camp equipment will in- 
clude at least 200 items. This is quite a few to plan for in advance, 
but all must be considered very carefully. The list has to be complete 
down to the can opener, which is one of the most important items in 
a foreign geophysical camp. 

There are some few areas left in the world where even arms for the 
protection of the working personnel are necessary. This will not be 
true in very many cases, but where there are hostile Indians or other 
peoples it will be necessary to negotiate with the government involved 
for a permit to import arms. This should all be arranged in advance. 

A word of caution might be given here regarding packing and ship- 
ping equipment. It is very necessary to provide proper customs 
declarations for the country of destination on all shipments. Equip- 
ment should be carefully and securely packed, and the contents of 
each box carefully listed so that the customs men cannot find discrep- 
ancies. Quite a few geophysical parties have been delayed in starting 
work because the equipment was not properly packed and declared 
for the customs to enter the country where the work was to be done. 

The supplying of even one seismic party is a considerable operation. 
Naturally the supplies of dynamite and casing will have to be ordered 
considerably in advance of the beginning of operations, and a suitable 
storage space for this dynamite will have to be arranged before its 
arrival. This should be as close to the scene of operations as is avail- 
able and the government will permit. A study of the food available 
for purchase for the technical personnel should be made in the country, 
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and it may be found necessary to ship in food for the party from the 
States. Some companies have sent men in dhead of the start of opera- 
tions in a very isolated area to clear the jungle and start farms to 
raise fresh food for technical personnel and labor and probably also 
to raise corn for the mules. This not only provides a dependable sup- 
ply, but also cuts down on the cost of transportation. Naturally a 
study of supplies of all kinds available at local markets should be 
made, and also the relative cost of locally purchasing supplies or of 
shipping in from the States should be analyzed. It has become cus- 
tomary to ship out most geophysical parties from the States with 
enough supplies of all types on hand for one year’s operations. This 
allows the party to operate six months or so and determine what sup- 
plies are needed and used most. They can then re-order at this time 
and receive the first supply shipment necessary sometime toward the 
end of the first year. Of course, this takes some very good planning 
beforehand to see that all necessary items are shipped with the regular 
party. Otherwise the party will end up sending frequent frantic cables 
for forgotten items to be shipped air express from the States. This is 
not only costly, but very probably will delay operations. 

Certainly choosing the personnel for the party is as important as 
choosing proper and sufficient equipment. The limited space does not 
allow for discussion of this point, but it is highly important, as the 
success of the operation depends almost entirely on these men. Be- 
sides the technical qualifications involved, the men must be hard, 
steady workers of a fairly even temperament, good natured enough 
to live congenially with the other men in the closely confined camp 
atmosphere. The men must be ingenious to take care of unusual 
situations. Such men are not easily found. 

Most of the above has had to do with scouting of the area, order- 
ing of supplies and equipment and choice of personnel. Some actual 
field preparation should be done in the area before the arrival of the 
party. Even in case the area is accessible by road, a camp should be 
prepared in advance and the first lines opened up for the geophysical 
work, so the party arriving can move right to the point where work 
is to start. It is common practice in this respect to send two or three 
of the men involved in the geophysical operations on ahead of the 
rest of the men and supplies so that they can supervise the opening 
up of the area and be thoroughly familiar with it by the time the re- 
mainder of the group arrives. Usually the party chief is among these 
men. The laborers should be hired in advance so that the full comple- 


| 

| 

| 

| 

|| 


THE PLANNING OF A FOREIGN GEOPHYSICAL OPERATION 441 


ment of men for all operations will be on hand upon arrival of the 
party. Some labor will already be working before the party actually 
arrives, and this is very beneficial, as these men are slightly broken 
in to working with the personnel of the party, and will be able to get 
operations under way much faster. It is very costly to delay the entire 
party of laborers and technical men while the first camp is built and 
the first road or line opened up for the geophysical work. A man who 
_is good at timing and guess-estimating should be excellent to have in 
charge of the whole project, as proper timing means a saving of 
money. 

To illustrate better the necessity for all of the above planning and 
timing it would be well to consider the cost of a seismic party in 
foreign operations. Naturally, this cost varies considerably, depend- 
ing on country, area, and operator, but it will be an unusual area 
where the minimum cost does not amount to at least $25,000 per 
month. This can run much higher in certain areas, where the operation 
of one seismic party can cost as much as a hundred thousand dollars 
a month. However, the average cost of operation of a seismic party is 
probably around $45,000 per month. The way a conscientious seismic 
party man calculates this out in the field, it means that it is costing 
more than $250 for every hour’s time he loses in the work. This 
thought usually puts much urgency on the workers on one of these 
parties, but it is easy to see that they cannot do an efficient job unless 
they have sufficient and proper equipment. It would not take many 
hours delay in any month to pay for an extra truck or an extra boat or 
an extra drill, or whatever else was necessary to speed up the work. 
If an extra truck unit saves the operator one day in a month, then 
it would pay for itself by being there that day. If a party is tied up for 
a whole week for lack of some essential spare part for the equipment— 
and this has happened more than once in foreign operations—then 
the operator has lost enough money to pay for considerable spare 
parts to be carried on hand. In other words, with 25 days in the average 
work month, and $45,000 per month the average cost for a seismic 
party, the operators can afford to pay for considerable advance pre- 
paratory work and extra equipment to keep from delaying the crew 
even a few days. The cost of the technical seismic equipment and 
personnel usually is only a small part of the cost of the whole opera- 
tion, so it is not real economy to cut down on the number of qualified 
technical men or the amount of equipment to be furnished on the 


party. 
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In the old days when all of us had lots to learn about operation of 
foreign geophysical parties, it usually took anywhere from a month 
to six months to get a party into actual operation after arrival in the 
country. This was always costly, but in those days it was considered 
a necessity. Now we all know that with proper planning and timing 
considerable money can be saved in carrying out the geophysical 
operations of a company. Costly delays are not always accidental, 
but can often be avoided by proper planning. 
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NOTES ON SHOT POINT PROCEDURE* 
F, J. WILLIAMSt 


ABSTRACT 


Some characteristics of seismic explosives are reviewed, and an attempt is made to 
outline optimum shot point procedure from the standpoint of efficiency, safety and 
quality of data. Care in preparation of charges is stressed, with special emphasis on 
adequate cap firing currents and charge assemblies. 

Suggestions are made in regard to the choice of the location for the up-hole time seis- 
mometer. 

It is concluded that adherence to well-defined specific shot point procedures is es- 
sential to the procurement of highest quality recordings. 


INTRODUCTION 


This paper is derived from a course in shot point procedure for 
shooters. The principle sources of information were: explosives manu- 
facturer’s literature, seismic data, and practical field experience. Al- 
though these notes were originally intended for company shooters, 
it was later conceived that certain portions might be of general inter- 
est to the seismic field. 

Such a study tends to make one agree with Sharpe in his observa- 
tion that; “Of the three physical processes involved in seismic explora- 
tion: namely, the initiation of seismic waves, their propagation, re- 
flection, refraction, and dispersion, and the recording of some func- 
tion of the motion of the surface, we possess the least satisfactory 
understanding of the initiation process.’! One is impressed by the 
frequency with which shot point data enter into the simplest of seismic 
computations. For instance it has-been observed that of the ten 
computational items listed on a typical record six are related to the 
shot point; that is, charge size, charge length, charge depth, time 
break, up-hole time, and near-surface velocity. 

In this paper it is assumed that the usual routine of handling 
shot point data in computations is adequate, and discussion is con- 
fined to shot point procedures which, though apparently of a minor 
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importance, nevertheless merit attention from the standpoint of 
quality recordings as well as shot point efficiency and safety. As a 
basis for determining optimum shot point procedure it seems neces- 
sary first to review some of the characteristic of seismic dynamites 
and detonators. 


SEISMIC DYNAMITES 


While many different types of powder were utilized in pioneer days, 
the exacting requirements of seismic exploration have narrowed the 
selection down to a product developed specifically for the seismic 
trade by leading explosives manufacturers. 

Frequently, this special powder is called High Velocity 60% L. 
F.. Gelatin. Here ‘‘High Velocity” refers to the consistently high rate 
at which the detonation propagates through the charge; “60%” is 
the proportion of nitroglycerin; and “L. F.”’ means low freezing. This 
product is composed chiefly of a nitroglycerin-nitrocellulose gelatin 
mixed with carbonaceous materials and an oxidizing salt. Several 
characteristics make this explosive especially adapted to seismograph 
work, such as: 

(1) High explosive power. 

(2) Great water resistance. 

(3) Effective detonation under great water pressure. 

(4) Suitable plasticity for convenient handling. 

(5) Density of 1.5 which makes it sink in water. 

(6) Rigid packaging which makes the resulting charge convenient 
to assemble and stiff enough to allow forcing through muck 
and sand in the shot hole. 

(7) Stable physical properties permitting fairly rough handling in 
the drill hole without undue hazard. 

An exploded charge releases energy in several different forms, 
and only a small part of the total energy initiates seismic waves. A 
portion of the energy is expended in crushing the shot hole walls and 
blowing water and mud out of the hole. Unfortunately from the seismic 
standpoint, most of the energy is converted into heat which is dissi- 
pated in raising the temperature of the material surrounding the point 
of shot. A number of gases; such as carbon dioxide, steam and a small 
volume of poisonous carbon monoxide are generated suddenly within 
the charge boundaries, and the heat of reaction raises these gases to a 
temperature of the order of 5000°F. Some idea of the enormity of this 
temperature may be gained by comparing it with the melting point 
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of steel, 2500°F. approximately, and the operating temperature of the 
electric-arc, about 6300°F'. 

When detonation transforms a solid stick of dynamite into a 
variety of gases, enormous pressures are suddenly developed in the 
immediate vicinity, for the gases naturally seek to occupy a much 
greater volume than a solid of equal weight. It it were possible to 
confine the resulting gases within the same dimensions as the original 
charge, pressures in the order of 700,000 pounds per square inch would 
be developed by a standard grade seismic dynamite. 

The detonation does not occur instantly, but proceeds along the 
column of explosive at a rate which may vary from about 6000 to 
30,000 feet per second, depending on the explosive and on the condi- 
tions under which it is detonated. In addition to the propagation of 
the detonation along the cartridge, another factor involving time is 
the momenta of the gases. If the explosion occurs in a restricted air- 
filled space, the gases produced by the detonation do not fill the whole 
space instantly, but a series of complex phenomena take place. High 
speed photography has revealed that a shock or compression wave 
first proceeds from the explosive, propagated by the air medium sur- 
rounding it, and that this is followed by the products of detonation, 
traveling at a lower velocity.” 

In seismic prospecting the explosive is, of course, ordinarily fired 
under water and the above conditions do not apply exactly. However, 
some investigations of explosion pressure waves in water have also 
been made.‘ In these experiments a dynamometer was placed at a 
distance of one meter from a charge under water. When the charge 
was fired two distinct pressure waves were recorded, having approxi- 
mately equal amplitudes, but occurring about .o5 second apart. These 
two pressure waves were attributed by Nash and Martin to the shock 
wave and the expansion of the products of detonation, respectively.® 

Additional evidence that the ground motion generated by the 
conventional discharge may be of a rather complex nature, and of 
comparatively long duration, has been provided by actual field ex- 


2 Nash, H. E. and Martin, J. M., “Significance of Some Fundamental Properties 
of Explosives, with Special Reference to Geophysical Prospecting,” Gropuysics, Vol. 
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3 Lawrence, R. W., “Mechanism of Detonation in Explosives,” Gropuysics, Vol. 


IX, No. 1, pp. 1-18, January, 1944. 
4 Marshall, Explosives, Vol. 2, P. Blakiston’s and Sons, Philadelphia (1917). 
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periments described by Sharpe. Among the results of these experi- 
ments he observes “‘a low frequency secondary wave of amplitude com- 
parable to that of the primary wave... .’” He finally concludes that 

“|. the true ground motion which results from the detonation of a _ 
charge of explosive in a drill hole is not a simple pulse, but a relatively 
long train of almost sinusoidal waves... .’”’ Constructive treatment 
of such effects of an explosion is clearly beyond the scope of this 
paper, but the significance of these findings is self-evident. Further 
experimentation relative of this phase of shot point kinetics is in order, 
with improved seismic resolution and interpretation as the objectives. 

Returning to the propagation of detonation within the charge it- 
self, it is noteworthy that the so called “high velocity” dynamite 
detonates consistently at a high propagation rate (about 20,000 feet 
per second) even when confined under water heads of about 350 feet. 
The importance of using large diameter cartridges is emphasized by 
the results of investigations conducted some time ago by Hercules 
Powder Company.’ In these tests it was found that the cartridge di- 
ameter affects both the rate of detonation and the completeness of 
propagation along the length of the column. For instance, 1-}"'X8” 
cartridges failed to detonate under a 200-foot head. By way of com- 
parison, 3” X16” cartridges detonated consistently at the high propaga- 
tion rate even though the water head was increased to as much as 
346 feet, and the same strength cap (No. 6) was used throughout. 

In view of the high temperatures generated, it is not surprising 
that under certain conditions dynamite elects to “burn” instead of 
carrying its normal function to completion. Experience indicates that 
“burning” is usually due to one or more of the following causes: 

(1) Use of a type of powder which is not designed to withstand 

the water pressures encountered. 

(2) Use of inadequate detonator. 

(3) Improper priming methods. 

(4) Preloaded charge left under great water pressure until soaked. 

(5) Detonator becomes accidentally separated from charge. 

(6) Firing of charge in one hole causes preloaded charge in ad- 

jacent hole to burn. 


6 Sharpe, Joseph A., “The Production of Elastic Waves by Explosion Pressures, IT. 
Results of Observations near an Exploding Charge,” GEopnysics, Vol. VII, No. 3, pp. 
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(7) Reloading too soon after previous shot. Especially is this 
likely in dry holes; this alone is sufficient reason for water 
tamping. 

Burning powder is best avoided by selecting an explosive that will 
perform properly under the water pressures encountered, and em- 
ploying proper priming methods. 

The time lag to allow between placing a charge in a hole following 
the shooting of a previous charge is somewhat indefinite because of 
the various factors involved. The Institute of Makers of Explosives 
Pamphlet, page 64, lays down the following rule: 


“Don’t explode a charge to spring or chamber a bore hole and 
then load anoth r charge into it before it has cooled sufficiently.” 


The above statement is necessarily vague, since the results depend 
to a large extent on the quantity of explosives used and the medium 
in which it is shot. In general it is well to beware of rock or dry holes. 
Whereas a dry rock hole will sometimes retain enough heat from the 
previous shot to detonate the following one, this is very unlikely for a 
clay hole filled with water. The importance of water as a heat dissipator 
in rock holes likewise becomes apparent. When the water ceases to 
boil and bubble, the hole temperature has certainly dropped below the 
boiling point of water. 

Several cases have been reported where the firing of a charge in 
one hole set off a charge in an adjacent hole. In some instances the 
charge in the adjacent hole “burned” and the generated gases forced 
the charge out of the hole before the cap was fired. Obviously such an 
occurrence can become extremely dangerous to personnel. 


SEISMIC CAPS 


The electric blasting cap used in seismograph work comprises an 
electrical igniting element or bridge wire surrounded by a priming 
composition and the “‘cap load,” all sealed compactly in a waterproof 
copper container or capsule. Two insulated copper wires are brought 
out from the igniting element through a water-proof seal at one end 
of the capsule for connection to the “firing line.”’ With sufficient elec- 
tric current the bridge wire heats up, igniting the priming composi- 
tion, which in turn sets off the cap load. There are two prevailing cap 
loads or strengths; No. 6 (medium strength) and No. 8 (high strength). 

For some purposes the function of the cap is “to impart to the 
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explosive an initial impulse of sufficient strength to set up a detonat- 
ing wave which will be propagated throughout the charge.’’* That is, 
the cap supplies the initial impact or quick blow needed to set off 
the dynamite charge. For quarry blasting, the above definition of the 
function of a cap is sufficient. In such cases the electric cap is simply 
a device enabling the shooter to stand at a safe distance while setting 
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Fic. 1. Curve showing the average firing characteristics of a seismic type cap. 


off a charge. The electric cap as used in seismic work, however, has 
another highly important function; namely, it must consistently fire the 
charge within very small time limits. Therefore the complete function 
of the seismic cap should be EFFICIENTLY TO DETONATE 
THE CHARGE WITHIN DEFINITE TIME LIMITS. 

Since the accuracy of calculations depends upon precisely timed 
seismic events, in same cases starting from the instant current is applied 
to the cap bridge wire, it is very important that the time characteristics 
of individual caps remain uniform or consistent throughout a given 
seismic survey. Available manufacturer’s tests indicate that at least 
some of the modern caps have been perfected to fulfill the above 
requirements, provided they are fired with a uniform value of electric 
current. This point is illustrated by Fig. 1, which shows the current- 


8 California Cap Company, Detonators for High Explosives, Manual. 
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time curve for a modern seismic cap.® According to the manufacturer, 

this type of cap is unique in that its detonation coincides precisely 

with the rupture of the bridge wire. Because of this desirable feature it 

. has largely displaced some of the older type caps where detonation 
lagged the rupture of bridge wire by a variable amount, depending on 
the firing current used. 

Upon referring to the graph, it will be noted that the cap detona- 
tion time in seconds has been plotted against the firing current in 
amperes illustrating that for any given firing current there is a cor- 
responding time-lag between the instant the current is applied and 
the resulting detonation. Thus for a firing current of 1.0 ampere, the 
curve indicates a firing time of about .0033 second. By way of com- 
parison, a firing current of 5.0 amperes corresponds to a firing time 
of only .coo4 second, or .0029 second less firing time than for the 
smaller current value. This is, of course, a serious time discrepancy 
in seismic work, where inaccuracies in the order of .oo05 to .cor second 
are the goal. It will be noted, however, that at about 4.0 amperes the 
cap curve flattens out and indicates a constant firing time for currents 
up to 10.0 amperes. For accurate results the firing equipment should 
therefore be designed—and maintained in adequate condition—to 
deliver consistently currents of at least 4.0 amperes to the cap. 
Some of the precautions which may be taken to insure good cap per- 
formance should be discussed. 

As all of the elements of the firing circuit are ordinarily connected 
in series and the sum of their individual resistances represents the 
total resistance of the circuit, it is a comparatively simple matter to 
design a blaster which will fire the cap at the proper current value. 
However, it should be stressed that in such calculations it is assumed 
that all splices have been carefully made and that no high resistance 
contacts are in the circuit. Poor splices and bad contacts can easily 
reduce the firing current well below the minimum value required for 
consistent detonation time. This undoubtedly is a major cause of 
misfires and the attendant wasted time, effort, and expense. It is a 

good rule, regardless of whether the lead wires are enameled, to scrape 
the ends of wires bright before making a splice. Loosely twisted, un- 
scraped joints may introduce surprisingly high resistances into the 
electrical circuit. 

There is also the “short circuit” type of defective joint which may 


® Copied from Hercules Powder Company “‘Vibrocap” firing curve. 
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prove just as detrimental to good cap performance as the high re- 
sistance joint, as where complete or partial contact is unintentionally 
established between the two leads. There may be actual contact be- 
tween wires or bare untaped splices may be exposed to the con- 
ductance of mud and water. The obvious solution is to make sure 
that all joints are properly taped and well separated. Every precau- 
tion should be taken to fire all caps within close current limits set 
high enough to insure uniform time delays. 


PREPARATION OF CHARGES 


One cap manufacturer stresses the importance of priming as fol- 
lows: “In making ready for a blast, there is no one step more im- 
portant than the proper preparation of the primer, by which is meant 
the cartridge of explosives into which the blasting cap is inserted and 
secured. The added assurance against misfire, not to mention the in- 
creased efficiency of the blast, will well repay extra time or effort 
expended in the careful preparation of this item.’ The detonation 
of the cap itself does not set off the whole body of the charge, but 
acts only upon the powder with which it is in intimate contact. Ex- 
perience has shown that proper preparation of the primer reduces 
misfires, increases the efficiency of the shot, and provides an addi- 
tional safety factor. 

The hole for the cap should be started near the middle of the 
cartridge and so angled that the cap lies as nearly parallel to the 
cartridge axis as possible. The cap should be entirely embedded, 
with the hole in the powder of correct size to allow the cap to enter 
snugly but freely. A bronze skewer or punch of proper size and with 
gentle taper is preferable. A skewer with too much taper will leave 
an air pocket around the cap which reduces the impulse to the powder. 
The hole should never be directed through the cartridge at right 
angles to its axis because in such position the cap is least effective in 
detonating the charge and it is also apt to become damaged or even 
disloged in loading operations. Also there is real danger of a premature 
detonation if the cap should happen to strike a hard obstacle in the 
hole during the loading process. 

The cap should be securely fastened in place so as to not become 
dislodged from the cartridge during loading operations. It is believed 
to be a questionable practice to try to save time by using poorly 


10 California Cap Company, op. cit. 
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prepared charges in which too much reliance is placed upon half- 
hitches. In such cases the wires are likely to get broken or at least the 
insulation may be damaged at the point where the wires cross in the 
half-hitch. Moreover, the cap leads are substantially weakened at 
the half-hitch so that the full tensile strength of the copper is not 
utilized. Another important point is that the cartridge shell not be 
discarded, since it provides additional rigidity for the cartridge and 
protection for the cap. ‘Hercules Spiralok” provides one method of 
eliminating half-hitches." First, remove the tubular shell or sleeve 
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Fic. 2. Single stick charge assembly in which half-hitches are eliminated. 


from the stick to be capped, insert the cap near the middle of the 
cartridge and then screw the shell back on over the cap wires, as 
shown in Fig. 2. For larger charges, where more than one stick of 
powder is used, the top shells will join at about the point where the 
cap enters the cartridge. In this case the top shell can be cut half in 
two and the charge assembled as indicated in Fig. 3. ; 
Large charges should ordinarily be assembled in tandem with the 
“primer” cartridge at the top end. With the cap placed in the top 
part of the charge, the main force of the explosion proceeds down- 
ward from the cap location so that the maximum amount of the 
charge is confined between the point of initial detonation and the 
bottom of the hole. In other words, with detonation starting at the top 
of a long series of end-to-end cartridges there is less chance of some 
of the sticks being thrown apart by the shock wave in the early stages 
of the explosion before the detonation has propagated to them. 


1t McAfee, J. C., “A New Shell for Seismic Shooting,” The Explosives Engineer, 
April, 1942. 
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Partially detonated charges are not only inefficient but also very 
dangerous, for some of the unexploded portions may be blown from 
the hole and detonate among personnel at the surface. 

This ‘“‘pushing aside” is most likely to occur at poor junctions be- 
tween sticks, a possibility which should be kept in mind when as- 
sembling long charges. The shooter should strive for a compact 
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Fic. 3. Multiple stick charge assembly without half-hitches. 


‘ assembly without dirt or other foreign matter between the cartridge 


ends. Incomplete explosions may also occur where the primer is care- 
lessly prepared or where the cap is too weak. For long charges in 
deep holes with high water pressures the No. 8 strength cap is recom- 
mended instead of the usual No. 6 strength cap. Under these condi- 
tions it is also advisable to use the largest practicable diameter 
cartridges, as discussed previously. 

For deep holes of 300 to 700 feet in depth, it has been often 
found easier to pull the charge down than to push it. For such condi- 
tions tamps are abandoned in favor of one of the various weight type 
loading devices. One loading arrangement consists essentially of a 
long iron weight (approximately 500 pounds) suspended on a cable 
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from an A-frame on the shooting truck and operated by a winch 
with brakes. As shown in Fig. 4, the load is arranged tandem fashion 
and taped along a suitable length of about No. g iron wire. The wire 
is bent into a hook at one end for insertion in the bottom end of a 
pipe extension on the lower end of the weight. Since the hole is some- 
times fouled by cap wires on the first shot, it is a good practice to clean 
the hole between shots. This is done by attaching a bull-plug in place 
of the loading pipe and running the weight down and up the hole. A 
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“spoon type” loading device (Fig. 5) is sometimes used for deep holes. 
This is similar to the above described weight arrangement except that 
the weight consists of a slotted iron cylinder into which the assembled 
charge is inserted. The charge is planted in the same manner as before. 

Primacord should be used to reduce the hazard of fouling high 
voltage wires when shooting holes near electric power lines. Fatal acci- 
dents have occurred in some instances” when the cap leads or firing 
line blew out of the shot hole and made contact with high voltage 


12 G, M. Kintz, Informal Discussion of Explosives Hazards on Seismograph Crews. 
Geropnysics, Vol. XI, No. 2 (April 1946), pp. 148-164. 
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lines. Primacord Detonating Fuse consists of an explosive core having 
a detonation rate of about 20,000 feet per second contained within a 
waterproofed textile covering. It is a powerful detonating agent for 
dynamite, and has been found especially useful in seismic work as a 
means for reducing the length of the firing leads. For instance, for an 
80-foot shot hole near a power line a 60-foot length of primacord can 
be attached to the charge at one end and capped at the other end, 
thus eliminating approximately 60 feet of dangerous firing leads. The 
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Fic. 5. Charge assembly using “spoon” loader. 


cap should be far enough below the surface so that the primacord or 
cap blast will not interfere with the normal functions of seismometers 
near the shot hole. 

In making up a primacord charge, a hole should be punched 
laterally through the top end of the cartridge. One end of the prima- 
cord is threaded through this hole and tied in a solid knot as shown 
in Fig. 6. The open end of the primacord should be wrapped with 
friction tape to exclude moisture, and the opposite end should be cut 
off squarely and kept dry. The cap may be attached to the primacord 
by means of one of the regularly supplied unions in accordance with 
the manufacturer’s instructions. It is very important that the prima- 
cord be cut off squarely and that the cap be seated firmly against it, 
since an air gap of only %” may cause a misfire. The assembly should 
be well taped at appropriate points to exclude moisture. The exact 
length of primacord used should be noted, of course, so that the proper 
correction for time-lag can be made. 
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UP-HOLE SEISMOMETER LOCATION 


The shot point seismometer has been termed “the most important 
seismometer in a reflection seismograph setup” by some authorities. 
The practice of placing this seismometer on the surface at the shot 
point to obtain the travel time from charge to surface is a common 
procedure. But since this up-hole time is one of the corrective factors 
used in near-surface time adjustments and since many subsequent 
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Fic. 6. Charge assembly using “‘primacord.” 


calculations for each record are affected by the location of this instru- 
ment, this detail should be discussed. 

The important thing to keep in mind at all times is that the up-hole 
time measured at the shot hole should be representative of the travel 
time through the corresponding near-surface interval directly under 
the regular recording seismometers near the hole. In other words, 
usually the shot point seismometer is placed directly at the top of 
the shot hole; and it is assumed that the up-hole travel time meas- 
ured applies under the recording spread as well; and the accuracy of 
time adjustments depends on how well this assumption holds true. 
Choice of location is a pertinent factor. For example, if the recording 
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spread seismometers are placed on a road and the shot point is located 
a short distance away in a field, an up-hole time recorded through the 
soft soil at the shot point is probably not truly representative of 
conditions along the hard packed road where the spread seismometers 
are resting. A more accurate measurement would result if the up-hole 
selsmometer were placed along the edge of the traveled roadway. 
Also, in some instances the shot point may be located on consider- 
ably higher or lower ground than the road bed where the spread is laid 
out so that an up-hole time measurement at the hole obviously fails as 
a measurement of conditions at the spread. Here too the up-hole seis- 
mometer should be placed on the road. 


CONCLUSIONS 


It is believed that adherence to the procedures outlined herein 
will improve the quality of seismic data to some degree and will 
increase the over-all efficiency of a field crew through the reduction of 
misfires. 
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ABSTRACT 


- The present paper is an attempt to take account of the lateral variations in velocity 
that are found usually by shooting a number of wells. More particularly, we show how, 
by tilting the linear distribution of velocity with depth in the proper way, we can tie 
the various wells together, accounting approximately (and usually to a sufficient ex- 
tent) for all lateral variations of velocity. The problem arising when three wells are 
shot is discussed in some detail. It is also shown that the lateral change of velocity can 
be determined approximately by shooting one well alone, and an approximate method 
of determining the tilt of the datum plane is discussed. 


INTRODUCTION 


The first refinement over the use of a simple constant velocity, 
with its associated straight line ray travel paths, is the linear increase 
of velocity with depth. In many cases in actual practice, this linear 
increase in velocity is sufficiently good so that it serves as a very 
excellent approximation to the actual situation. A very impor- 
tant advantage of the use of the linear increase in velocity with depth 
is that we are not limited in our interpretation procedure to the 
calculation of small dips. It is quite evident that where constant aver- 
age velocities are used, dips are limited to 90°; whereas, with the line- 
arly increasing velocity-depth relationships, we can deal with dips 
greater than go°. It also seems important in some cases to have pro- 
cedures set up that are uniform for the small dips and for the large 
dips, so that no judgment may enter in the problem in transferring 
from a low dip problem to the high dip problem. If the linear distribu- 
tion of velicity with depth is used, the transition from low dip prob- 
lems to high dip problems is automatically taken care of. 

However, some criticism has been made of the use of linearly in- 
creasing velocities with depth, or the circular arc ray path, on the 
grounds that full account of lateral variations is not taken care of. 
In the present paper, a more or less complete approximate solution 
to this problem is presented. The solution is approximate, but it has 
the great advantage that it is very simple to apply. 

We discuss first the case where we have the velocity surveys from 


* Read by title at the annual meeting of the Society, 1945. 
+ United Geophysical Co., Pasadena, California. 
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three wells, located on three corners of a triangle in a given area. 
We also discuss the case of two well surveys and of one survey. It is 
quite possible, using one survey alone, to obtain data by which the 
tilt of the datum plane can be calculated as well as the datum velocity 
and the rate of change of velocity perpendicular to this datum plane. 

For the case where the dips are all very gentle, the method pre- 
sented by Daly! can be recommended, but for steep dips a method 
similar to the present one must be resorted to. Where both steep and 
gentle dips are involved, the present method seems to be required 
in order to preserve the continuity of the interpretation. 


CASE OF THREE WELLS 


We suppose that the three wells do not lie in a straight line so that 
they determine a triangle of appreciable area. We suppose also that 
these wells have been shot, and that from the data we have obtained 
values of the datum velocity, Vz, and of the rate of increase of velocity 
with depth, (a). These values for the three wells are (Vai, Vae, Vas} a1, a2, 
a3). It is desirable for our purpose that the values of a and of Va be 
determined by the method described by Legge & Rupnik.? We then 
take as our a-value, the average value of three a’s, that is to say, 
a= (a+ 

Legge & Rupnik,? show that the datum velocity, Va, is related to 
the a-value through the observed values Z; and #;, where /; is the 
vertical travel time from the datum to the depth Z;. It is also found 
that Va is almost a linear function of @ in the neighborhood of the 
“best fit”’ value. Thus, if we adjust the a-values by replacing the in- 
dividual a-values by the average, a, then we have to adjust also the 
datum velocity values according to the above relationship. Suppose 
these adjusted datum velocities are respectively Vin, Vae, Vaz, which 
are equal respectively to Va+a61, Va+ade, Vat+ad3, where 6;+ 52+ 63=0. 
Then for each of the three wells, the velocity is equal to Vat+aé;+<aZ, 
where 1=1, 2, 3. Thus, 6; is, for the i well, the change in elevation 
required to pass from the level datum up or down to the elevation 
at which the velocity is the constant value Va, where Va=(VaitVae 
+Vuas)/3. When 6; is positive, the datum point at this well must be 
elevated in order to ‘‘climb up” to the place where the velocity has 
the value V,. When 6; is negative, we must “dig down’ to the place 
where the velocity is the datum velocity Va. 


1 John W. Daly, Bull. A.A.P.G., Vol. 28, No. 5 (May, 1944), pp. 615-628. 
2 Geopuysics, Vol. VIII, No. 4 (Oct., 1943), pp. 356-361. 
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We then have for the three wells the co-ordinates (x;, y;, 6;) where 
i=1, 2, 3. The x and y co-ordinates are then say the north and east 
co-ordinates, whereas, the 6 co-ordinate indicates the elevation of the 
datum at the well. Thus, in this way, three points are determined 
at the three wells through which the datum plane can be drawn. 
Along this datum plane, we may assume the velocity has the con- 
stant value Vz. Now, unless all the 6’s are zero, this datum plane will 
be tilted, but, of course, if two of the 6’s are zero, then the third must 
be zero since their sum is zero. It is, however, quite possible for one 
of the 6’s to be zero and the other two oposite in sign. In this case, the 
strike of the dipping datum plane passes through the point for which 
5 is zero, and cuts the line joining the other two wells in the middle. 
In case none of the 6’s are zero, two will always be of one sign and the 
other will be of the opposite sign. Suppose, for example, that 6; is 
negative, but that 62 and 63; are both positive. Then, in passing from 
well #2 to well #1, we must decrease the elevation of the datum plane 
due to its tilt and cross the level datum plane. Likewise, in passing 
from well #3 to well #1, we also decrease the elevation of the tilted 
datum plane and pass through the original level datum plane. These 
two points where we intersect the original level datum plane con- 
stitute two points on the strike line of the tilting datum plane. 
These two points are readily found by graphical or other means, and 
thus the strike of the tilting datum plane is determined. When we 
have determined the strike, we can very readily determine the dip 
angle of the datum plane which we shall designate from now on as d. 

The measurement of a has been referred to a vertical line. We now 
have to consider the fact that the direction is no longer vertical, but 
is also tilted at an angle, d. Now we consider that the angle, d, of the 
tilt of the datum plane is very small, and so the correction of the 
a-value is likewise small; however, if such a correction is to be made, 
we have only to pass from the a-value we obtained from the vertical 
reference, which is from a to a/cos d. 

We have first found a tilting datum plane along which the 
datum velocity has a constant value Vz, and for which the a-value is 
a/cos d. We next correct all reflection times to this tilting datum 
plane; from here on we carry out our computations in exactly the 
same manner as if this datum plane were level. The words “depth 
below datum” must now be replaced by the term “‘distance from the 
tilting datum plane.” It is fairly evident that this adjustment will 
alter both the calculated depth and dip of the reflector. 
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CASE OF N>3 (WELLS) 


In such a case, we join the wells together by a series of straight 
lines forming a triangular network. Unless there is too much variation 
_in the a-values over the area, it is desirable to use as an a-value, the 

average of all the a-values for the set of wells. In this way, we obtain 
a datum surface which is a series of triangular planes all fitting con- 
tinuously together. It is, of course, important that the angles at the 
edges of these triangles between the planes be small, otherwise, too 
large discontinuities are introduced into the interpretation and the 
consequent adjustments have to be made at all the edges. The process 
for each triangular segment is exactly the same as that described 
above in the case where we have three wells only. 


CASE OF TWO WELLS 


In case we have only two wells, we have to make some sort of an 
assumption with regard to the strike of the datum plane, as it is 
obvious that two wells are not sufficient to determine its strike. If 
we have no evidence whatsoever, it seems best to assume the strike 
is perpendicular to the line joining the two wells. The discussion is 
entirely similar to that given above for the three well case. 


CASE OF ONE WELL 


In the case where only one well has been shot by the ordinary 
methods, of course, we have no information whatsoever with regard 
to the lateral variations of velocity in the section. We may, of course, 
assume a hypothetical variation based upon our knowledge of the 
regional geology, but this is likely to be somewhat risky. 

It is possible to shoot a single well in such a way that the tilt of 
the datum plane can be found from the data if the well is shot from 
several shot points surrounding the well. It is found that for each 
depth we can suppose that the well geophone has been changed into 
a source instead of a receiver, and that geophones are placed at the 
surface in place of the shots, then the distribution of times due to 
the tilting datum plane will be such that the minimum travel time 
from a given depth will not be exactly at the well, but will be offset 
to one side somewhat. Let this offset amount be indicated by 6. 
Then it is not difficult to show that we have the following approxi- 
mate relation: 


sind = 26V 
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where 


d=dip angle 
Va,=datum velocity 
a,=increase of velocity with depth 
Z;=depth of well geophone. 
Va. and a, are found approximately by removing the effect of dip by 
averaging values taken on all sides of the well. 6 is measured more or 
less directly from the data. Z; is given by the depth of the geophone 
in the well, thus leaving only the dip angle, d, to be determined. 

The above method of handling the lateral variation at one well 
of course involves some rather rough approximations. However, it 
seems sufficiently accurate in view of the inaccuracies of the data 
likely to be obtainable in such a case. 

It is suggested that in shooting such a well, the shot points may 
be located on rays issuing from the well making angles of 120° with one 
another. Probably a minimum of two shot points on each of the three 
rays would be required. Such a survey, of course, would be con- 
ducted only for a well in a rather complicated area where strong 
lateral variations of velocity are very likely to occur. ; 


OPTICS OF RECORDING OSCILLOGRAPHS* 
D. H. WEINSTEINt 


ABSTRACT 


This note reviews the results of A. C. Hardy and applies them to the design 
of lens and mirror systems of recording oscillographs, such as those used in seismic 
prospecting. Requirements are given for maximum intensity of illumination in the re- 
cording paper consistent with small moment of inertia of the moving mirror. 


In many of the moving mirror oscillographs used at present, the 
physical dimensions of the mirror are of primary importance. The 
trend has been to make galvanometers smaller and more sensitive, yet 
the necessity of maintaining adequate illumination on the recording 
paper makes it essential for the designer to consider the relation- 


Source Stop Mirror Paper 
f 

| 

! 


Fic. 1. Oscillograph optical system of first class, using spherical 
galvanometer mirror. 


ships between physical size of the mirror and illumination. This sub- 
ject has been well set forth by Hardy.! The purpose of this note is 
to review his results and show how they may be applied with ad- 
vantage. 

It is convenient to separate galvanometers into two classes accord- 
ing to their optical systems. The first class consists of a centered 
system composed of spherical lenses. The second class consists of a 
system composed of cylindrical lenses. 


* Read at the Chicago and Houston meetings of the Society, 1946. 

+ The Superior Oil Co.; Houston, Texas. 

1 A. C. Hardy, JOSA and RSI, 14 pp. 505, 1921. The Principles of Optics, Hardy 
and Perrin, McGraw-Hill, 1932. 
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All instruments of the first class may be represented as a spherical 
galvanometer mirror. Such a system is shown in Fig. 1, where the 
mirror is represented in the system as a stop and two lenses. A source 
stop is also shown, as the source itself may be too large. 

The illumination in this case is given by 


I, = BA/om (1) 


where B is the brightness of the source, A the area of the mirror, and 
Ym the distance from the mirror to the paper. 


Gyi.Lens Paper 
Source Mirror 


ALN 
Wi 


Fic. 2. Oscillograph optical system of second class, using cylindrical 
galvanometer mirror and cylindrical lens. 


If B and »v,, are fixed, then the illumination varies with the area of 
the mirror. To obtain the greatest area for a given moment of inertia 
the mirror should be longer than it is wide. Systems of the first class 
are not widely used in recording because they do not afford the il- 
lumination possible with systems of the second class. 

All instruments of the second class may be represented by a spheri- 
cal galvanometer mirror and a cylindrical lens. It is simpler, however, 
to represent it by the system shown in Fig. 2, which consists of a 
filament source, a cylindrical mirror whose axis is parallel to the axis of 
the source, andacylindrical lens whose axisis at right angles to theaxis of 
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the mirror. For convenience the mirror is shown as an aperture and 
two cylindrical lenses. The cylindrical lens is focused on the mirror, 
so the dimension of the mirror along its axis does not act as a limiting 
stop. In other words, from the viewpoint of the cylindrical lens, the 
source of light is the mirror itself, and so changing its length merely 
changes its apparent size and does not affect its brightness. The width 
of the mirror still acts as a stop. If this procedure were also applied 
to the width, rotating the mirror would no longer move the image. 

The filament is assumed long enough to fill the cylindrical lens 
with light. If it is not, then a condensing lens may be used in con- 
junction with the source so as to enlarge it. The image of the source 
is, of course, focused on the mirror, and care must be taken to make it 
larger than the mirror. 

The expression for the illumination with a system of the second 
class is given by 


= B(wm/Vm)(We/V) (2) 


where B is the brightness of the source, wm the width of the mirror, 
Ym its distance from the paper, w, the width of the cylindrical lens, 
and 2, its distance from the paper. 

We see that the length of the mirror /,, does not enter into the ex- 
pression for J, for the reason given above. Nevertheless, the length 
1, is of importance in that it is desirable to make the image or: the 
recording paper very nearly square. The dimensions of the spot are 
given by 

Wi = 


where w,= width of source or source stop. 

From (3) it follows that if J, is decreased and all the other optical 
constants held fixed, /; is proportionally decreased. This decrease in 
the length of the image spot will result in the recording of the undevi- 
ated line to appear as though illumination has been lost by an amount 
proportional to the decrease in J. This is simply due to the fact that, 
if the paper is moving at a constant speed, the time during which any 
point on the film is exposed will be decreased. However, when the 
galvanometer trace moves at right angles to the undeviated line, then 
the trace will be just as dark as when /,, is not decreased, but of course 
the trace will be thinner. 

It must be mentioned also in this connection that when /,, is de- 


(3) 
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creased there will be a proportionate decrease in the illumination of 
the cylindrical lens. This is perfectly proper since the image is reduced 
in proportion and thus the illumination of the image is maintained as 
previously stated. 

The ratio of /m/wm in many of the galvanometer mirrors now in use 
varies from about 1.5 to 7, with /,, running from about .o80 to .250 
inches. Assuming that the spot definition is satisfactory with such 
mirrors, it may still be an advantage to reduce /,, to say, one fourth 
its value. To maintain the same image size, it would then be necessary 
to increase 2, so that v./(%m—) is multiplied by four. 

In a typical case, ».=1 and v,=11. Then v, would be increased 
from 1 inch to 3.15 inches. Of course, the ratio w./v., which is ordi- 
narily about 3, must be kept constant. This calls for a cylindrical lens 
of about 2 inches aperture. This increase in aperture will exactly com- 
pensate the decrease in illumination on the lens and thus produce an 
image of exactly the same size and illumination as in the previous case. 

Where the moment of inertia of the mirror does not need to be re- 
duced, the illumination may be increased. This can be done by reduc- 
ing /,, and making the appropriate increase in Wm. 


AN OBSERVATIONAL METHOD TO OVERCOME ZERO 
DRIFT ERROR IN FIELD INSTRUMENTS* 


IRWIN ROMANt 


ABSTRACT 


By using a four-step oscillating method of observation, four values are obtained 
for the increment in the measured value between consecutive stations. If the drift or 
“zero creep” of the instrument is regular, even when large, or not linear, this furnishes 
a means of improving the reliability of the readings and permits the detection of er- 
ratic readings. In areal surveys, a primary net of values may be established by this 
oscillating method and this net used as a base into which to adjust chains of observa- 
tions. Intermediate values may sometimes be desirable, such values being adjusted 
into the primary net by the oscillating method. An example shows the application of 
the method to a gravimetric survey. 


INTRODUCTION 


In certain types of “equilibrium” measuring instruments, the pre- 
cision of results is affected by a gradual creep or drift of the reference 
base. Among geophysical instruments, two instruments are especially 
subject to such uncertainties, the magnetometer and the gravimeter. 
The drift is usually the result of numerous effects, such as those due 
to temperature, earth tides, or mechanical shocks. For many pur- 
poses this drift is not important, as the measured differences are large 
compared with the amount of drift or the drift may be roughly approx- 
imated by some standard set of values. In many surveys corrections or 
intrinsic errors of observation are of the same order as the drift un- 
certainties, and interpretations of readings are not aided by an at- 
tempt to correct for drift. However, for some surveys the fineness of 
the observation may justify the use of control methods to eliminate 
drift or diurnal variations. It is for this type of survey that the present 
method was evolved and judicious considerations must be used in 
deciding when the extra time and cost are justified by the increased 
accuracy of the results. 

A secondary advantage of the present method lies in its ability to 
detect obvious errors and often to permit a form of adjustment of 

* Published with permission of the Director, U. S. Dept. of Interior, Bureau of 


Mines. Read at the Chicago meeting of the Society, April, 1946. 
{ Geophysicst, U. S. Dept. of Interior, Bureau of Mines. 
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data, keeping clearly in mind that the original data may not justify 
further adjustment. 

The method is applied to the original readings before corrections 
are made for predictable effects and before the readings are converted 
to physical values by means of the sensitivity factor. 

The method will be referred to as the “four-point oscillating” or 
more simply as the “oscillating” method. A “primary net” of closed 
loops is observed by the oscillating method and the closure error dis- 
tributed uniformly over the links of each loop. After values have been 
“established” for each station on the primary net, a “‘secondary net” 
is observed by the oscillating method, each line of this net beginning 
and terminating at a station of the primary net. The difference of 
values at the ends of each line is then adjusted uniformly over its 
links to agree with the difference as established on the primary net. 
Finally, “chains” are observed, by non-oscillating readings along paths 
between end points, each of which lies on either a primary or a second- 
ary net. Assuming the values of the end points to be established by 
the primary or secondary nets, any discrepancy is adjusted uniformly 
over the “‘links” of the ‘“‘chain.”’ 

After values are established for each point of the survey, correc- 
tions are made and the readings are reduced to the desired physical 
units by means of the sensitivity factor. 


THE OSCILLATING METHOD 


The principle of the oscillating method is illustrated schematically 
in Fig. 1. Roughly, the stations are observed in groups of four, witha 
return to the second station of a group as the first station of a new 
group. This procedure is modified slightly at the beginning and end of 
the line in order to gain the desired objective. If the initial stations 
are A, B, C, D, these stations are observed in the order: 


AB, ABC, ABCD, BCDe, CDef, D, etc., 


the lower case letters indicating that the stations are not in this group. 
For the intermediate stations K, L, M, N, the order is: 


K, ijKL, j]KLM, KLMN, LMN,, MNp}q, N. 
At the end of the line, for stations W, X, Y, Z, the order is: 
W, uwWX, WXY, WXYZ, XYZ, YZ. 


In planning a line, especially when the observations are to be 
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terminated for the day, the end readings should be considered care- 
fully, as extra stations are of little value. 

For a loop in the primary net, the terminal point Z coincides with 
the initial point A. For a line in the secondary net, both the initial 
point A and the terminal point Z should lie on the primary net. For 
chains, not observed by the oscillating method, each end point of the 
chain should lie on either the primary or the second net. 

In the oscillating method, a specific link is observed four times. 


ERMINAL 


INTERMEDIATE 


INITIAL 
vu vw x 2 
| 
IJK MN P 


Fic. 1. Order of observation for 4-point oscillating system. 


Thus, the step LM occurs in the sequence LMKLMNLM four times, 
with each of the first two M observations between two L observations 
and each of the last two L observations between two M observations. 
If the drift curve of the instrument is known with sufficient accuracy, 
we may interpolate the reading which the instrument would show at 
one station of the pair at the time the other station was actually read. 
From the actual reading at one station and a fiducial* value for the 
other station, the difference in readings between the two stations of a 
pair can be determined. In the present sequence, if the drift curve is 
known, six values can be determined, one for each of the three ob- 
servations at L and one for each of the three observations at M. Since 


* Fiducial is used in the mathematical sense of a reference position not observed; 
it is fixed for each specific set of observations. 
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the drift of the instrument is usually not known, we may assume that 
it is linear between successive observations at the same station and 
determine the fiducial value by linear interpolation. In the present 
sequence four values can be determined, ene for each of the first two 
observations at M and one for each of the last two observations at L. 
In this method, successive readings at a station must be made with 
reasonably short elapsed time. . 
Numerically, consider the set of readings LML read as follows: 

at time o, the reading at L is Li. 

at time ¢,, the reading at M is M. 

at time ¢:+¢, the reading at L is Le. 

Assuming that the instrument has a uniform drift between the 
readings L, and Le, we may interpolate linearly for the value L’ which 
the instrument would show at time 4, when M is actually read. This 
value is 


t tol, + 


ty + te ty + be 


With this fiducial value at ZL for the time when M is read, the 
change in reading from L to M is determined as M—L’. Graphically, 
the same result may be obtained. If each reading is plotted as a func- 
tion of the time of observation, successive readings at each station 
may be connected by straight lines. The ordinate of one line at the 
time another station is read gives the fiducial value at the former sta- 
tion, and the vertical distance from a point of reading to another line 
gives a determination of the difference in values between the two 
stations concerned. 

After the successive differences have been determined, these are 
added to the accepted value at the first station to furnish an accepted 
value for each of the other stations of the set. Refinements, such as 
adjustment of residuals, closure adjustments, error corrections and 
smoothing may be applied to the differences before accepted values are 
determined. 

From the accepted values at each station, the value of the drift 
may be determined for each observation by subtracting the accepted 
value from the observed value. This drift curve serves as a check on 
the behavior of the instrument as well as the basis for detection of 
errors and for smoothing of the data. 
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ILLUSTRATIVE EXAMPLES 


To illustrate the method, consider primary loop B of a specific 
gravimetric survey, the observations of which are shown in Table 1 
and Fig. 2. The reduction of the data is shown in Table 2. In the first 
column of Table 2 are shown the initial station (coded as a) and the 
terminal station (coded as b) of each step in the loop. The time of each 
observation in the group is shown. The column headed “AT” is the 
time interval in minutes between the consecutive readings shown in 


TABLE 1. OBSERVATIONS—PRIMARY Loop B 


Station Time Reading Station Time Reading 
fe) 9:20 5 11:44 16.3 
I 9:27 26.5 6 11°49 21.0 
° 9:34 4 30.7 
I 9:41 26.4 5 £2502 16.5 
2 9:50 12.8 6 12205 2052 

7 $2322 
° 10:02 1.0 5 T2235 15.6 
I 10:08 26.3 6 12:30 2022 
2 10:16 12.5 7 12:42 30.5 
3 10:22 8.2 8 12:48 37-5 
I 10:33 28.0 6 1257 21.0 
2 10:42 14.2 7 1:03 30.9 
3 10:49 8.7 8 1:09 37:0 
4 10:57 30.8 ° 1:18 251 
2 11:08 5 7 1:28 30.7 
7.9 8 1°34 
4 11:20 30.9 ° 1243 
5 16.2 
4 11:38 30.8 ° 2:01 t.7 


the next column. The fiducial values are the values interpolated 
linearly between the adjacent lines of the column “readings” for the 
time of the line of the entry. The primes denote this interpolation. 
The interpolated values have been carried to an extra decimal to 
allow for forcing errors. Under the columns ‘‘d-a”’ are the individual 
values of the increment in reading from the initial to the terminal 
point of the step, as determined from reading and fiducial values by 
direct subtraction of a from b’ or a’ from b, and the arithmetic mean 
value for the group. Under “Deviation” are shown the individual ex- 
cesses of the increments from the mean of the group and the root- 
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mean-square of each group. The drift is indicated in the final column 
and is obtained by subtracting for a station the calculated value ob- 
tained as described below from the actual reading at that station. The 
“calculated” value at Station o is the first reading at that station. The 
calculated value of Station 1 is the calculated value of Station o in- 


Ob. 


READINGS 
nN 


\ 
dn 


| ---------4 


TIME 


Fic. 2. Observation—primary loop, B, of table 1 and Fig. 7. 


creased by the determined value of the mean increment from o to tr. 
In this way the calculated values are found to be as shown in Table 4. 

If the increment from o to 8, —34.97, is added to 8, the value for 
o is found to be 2.06 instead of 1.40, indicating a closure error of 
+0.66 scale divisions. If this error is acceptable, it may be distributed 
uniformly over the steps. Any remainder is distributed over the steps 
with the largest RMS deviation. Thus, the steps O-1, 2-3, 4-5, 6-7, 
7-8, 8-0, are each decreased by 0.07 and the steps 1-2, 3-4, 5-6. are each 
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TABLE 2. REDUCTION OF DATA—PRIMARY Loop B 
At Fiducial Value b—a Deviation 
Sta. | Code} Time ie. | Sua Drift 
Div. || Code | Reading || Individ. | Mean || Individ. | RMS 
° a 9:20 7 : ° 
I b 9:27 7 26.5 a’ 1.25 +25.25 —0.06 —o.21 
a 9:34 7 5.2 b’ 26.45 +25.35 | +25.31 |} +0.04 —0.30 
b 9:41 21 26.4 a’ 1.08 +25.32 +o.o1 | 0.04 |} —o0.31 
6 b’ 26.32 +25.32 +0.01 —0.40 
10:0 26.3 —o.4I 
I a 0:41 9 | 26.4 || 
2 b 9:50 18 12.8 a’ 26.37 —13.57 +o. 28 —0.06 
a | 10:08 8 | 26.3 b’ 12.80 —13.50 | —13.85 |} +0.35 
b 10:16 17 12.8 a’ 26.84 —14.04 —o.19 | 0.32 || —0.06 
10:33 9 28.0 b 13.72 —14.28 —0.43 
10:42 14.2 1.34 
2 a | 10:16; 6 | 12.8 | 
3 b 10:22 | 20 8.2 a’ ¥3.22 — 4.92 | +0.40 +0.66 
a 10:42 7 14.2 |) 8.57 95.32 |||) 
| b 10:49 19 8.7 |i a’ 14.01 — 5.31 | +0.01 | 0.26 || +1.16 
11:08 6 13.5 || 8.09 — 5.41 —0.09 
| | i 0.3 
3 | ro:49 ; 8 | 8.7 
4 b | 10:57; 17 | 30.8 | a’ 8.44 |, +22.36 —0.52 +0.38 
| T1314 | 6 7.9 || 8 30.87 |, +22.97 | +22.88 || +0.09 
| | 23 7.84 |; +23.06 +o.18 | 0.31 || +0.48 
| 11:33 | 5 | 77 | 30.83 | +23.13 +o. 25 
| 11:38 | | 30.8 | | 0.3 
4 | @ | 11:20 5 | 30.9 | | | 
| 8 16.2 a’ 30.87 | —14.67 —o.19 +o. 26 
a | 11:38 6 | 30.8 |: b’ 16.27 || —14.53 | —14.48 || —0.05 | 0.15 
b | 11344 13 | 16.3 a 30.77. || —14.47 +o.o1 +0.36 
| | 11:57 6 | | b’ | 16.44 |, —14.26 +0.22 
12:03 1 | 0.5 
5 a | 11:44 5 | 16.3 || | 
6 b 11:49 14 | 21.0 a’ 16.35 + 4.65 —0.34 +0.07 
a 12:03 12 16.5 b’ ci ae § + 4.61 | + 4.99 || —0.38 | 0.40 
b 12:15 16 21.2]; @ 10.11 + 5.09 +o.10 +0.27 
5 15.6 b’ 21.20 + 5.60 +0.61 
12:3 21.2 0.27 
6 a 12:15 7 2.21 
7 b 12:22 14 | 31.2 || a’ 21.20 +10.00 +0.33 +0.60 
a 12:36 | 6 21.2 |) b’ 30.71 + 9.51 | + 9.67 —0.16 | 0.25 
| 6b 12:42) 1§ | 30.5| + 9.36 | —0.31 —o.10 
| 6 21.0 |; b 30.79 + 9.79 +o.12 
| 1:03 | | 30.9 |} 0.30 
f a 12:42 6 | 30.5 | 
8 b 12:48 a’ 30.61 + 6.80 +0.46 +0.47 
a 1:03 6 | 30.9 |! b’ 37-14 + 6.24 | + 6.43 || —o.19 | 0.29 
b 1:09 tq. | 37-0 | a’ 30.85 + 6.15 —o.28 —0.03 
: 1:28 6 | 30.7 |) b 37-15 + 6.45 +0.02 jae 
37-2 |I 6.17 
354 1:09 37.0 |! 
° b 1:18 9 2.1 || a’ 37-07 —34.07 ° +o.70 
a 1:34 | 16 | 37.2 || b’ 2.10 —35.10 | —34.97 || —0.13 | 0.08 
b 1:43 @ 37.00 —34.90 +0.07 +0.70 
a 36.8 b 1.90 —34.90 +0.07 —0.23 
b 2:01 | +co.30 


decreased by 0.08. This leads to the “closed” values shown in Table 3. 
The closed values may not be more accurate than the calculated 
values but eliminate the closure error, a process justified if the closure 


error is permissibly small. 
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Inspection of Table 2 shows that the readings vary from 1.0 to 
37.5, a range of 36.5. The step differences vary numerically from 4.99 
to 34.97, the total of the numerical values being 137.90, distributed 
over nine steps, or an average, without regard to sign, of 15.32 per 
step. The deviations of the individual step determinations from the 
mean of each corresponding group vary from zero to 0.61. The root- 
mean-square of the deviations in each group varies from 0.04 to 0.40 
with a mean value of 0.23. The closure error represents about 0.5 per- 
cent of the aggregate step differences, about 2 percent of the range of 
the readings, and about 4 percent of the average step difference. 


TABLE 3. CALCULATED VALUES—Prmary Loop B 


Calculated Closed 
Station Value Value 
° 1.40 1.40 
I 26.71 26.64 
2 12.86 12.7% 
3 7-54 7 +32 
4 30.42 30.12 
5 15-94 15-57 
6 20.93 20.48 
7 30.60 30.08 
8 37.03 36.44 
° 2.06 1.40 


The values shown for this loop are typical. Adjustment of the 
readings has not indicated sufficient improvement to warrant the ef- 
fort. However, in some nets, obvious errors are indicated. In Fig. 3, 
the values for station o read between 1 and 2 o’clock all seem to be 
high. If they were lowered by one whole unit, the drift curve would 
be slightly smoother and the closure error would be changed from 
+0.66 to —0.34. As the observations at 7 and 8 are not affected, the 
discrepancy would seem to represent a reading error, not an instru- 
ment drift. On the other hand, the hump between 10 and 11 o’clock 
is indicated on stations 1, 2, 3 and also probably at 4, which may have 
been read one unit too low at 10:57. Thus, an instrument disturbance 
is indicated at this part of the survey. This hump seems typical, as 
it has a generally sinusoidal shape, such as would result from a shock 
with sufficient energy to displace the reference system. The irregular- 
ity between 12 and 1 o’clock would seem to indicate an interval of 
seismic disturbances, such as those due to nearby mine equipment in 


= 


474 IRWIN ROMAN 


motion. However, no general trend or apparent errors are indicated 
with clarity and the irregularities may be considered as inherent. 

It is sometimes possible to detect gross errors, as in loop number 3 
of the same survey as the one just discussed. The drift curve shown in 
Fig. 4 indicates several minor errors and a gross error. The gross error 
occurs at the terminal readings at station 0, and it is evident that the 
readings should properly have been two whole scale divisions larger. 
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Fic. 4. Drift curve of primary loop, C. 


Correction of this error decreases the closure error from — 2.22 to 
—o.22 and shifts the drift points from the circles to the crosses in Fig. 
4, improving the appearance considerably. 

In each of these two examples the last station of the loop was 
affected. This may have been due to chance or may have a psycho- 
logical explanation. If an adjustment is made at a station other than 
the initial or the terminal, the effects are more complicated. Each 
such station adjustment would affect two links and have variable ef- 
fects on the deviations. Also, all subsequent calculations of station 
values would be affected, thus shifting the drift curve, as well as modi- 
fying a portion of it. 

In the two loops just discussed the drift was relatively small. In 
Fig. 5 the drift curve is shown for loop 1 of a different gravimetric 
survey involving a larger drift. Except for the observations at 2:53 
and 3:39, the drift curve is fairly smooth. As the loop was not closed, 
no closure error was obtainable. The individual deviations have a 


a lic. 3. Drift curve of primary loop, B. 
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Fic. 5. Drift curve of primary loop, J. 


maximum of 0.32 scale divisions with a root-mean square of each 
group having a maximum of 0.25 scale divisions. The drift curve as 
calculated, showed a maximum value of 5.54 divisions with the drift 
of the last station 4.60 divisions. 


SECONDARY NET 


To determine the degree of reliability of gravimeter readings, 
several loops have been analyzed for comparison of the differences in 
readings along primary loop and along a secondary loop. The loops 
selected may be considered as typical and the conclusions will con- 
stitute a guide rather than a statistical summary. 

A secondary loop is measured from one point of the primary net to 
another and is usually not closed. However, we may consider it as 
closed by way of the primary net, which has been established by closed 
loops. 

Secondary Loop 1. Some of the important data of secondary loop 1 
are shown in Table 4. It may be seen that on the basis of a single 
reading at each station, the step C6-5 may have values varying be- 
tween +6.1 and +3.5. Based on the average readings on a particular 
run, the step may have a value of +3.60 or +5.56. Based on the 
oscillating method, the step has the value +5.07 on the primary net 
and +4.28 on the secondary net, representing a closure error of —0.79 
instead of a value between —1.2 and — 2.6as based on single readings. 

Secondary Loop 2. Values for secondary loop 2 are shown in con- 
densed form in Table 5. For this link, the closure may have values 
from +o0.1 to +1.3 based on single readings, with an average of 
+0.66. Based on the oscillating method, the closure is —0.65. On this 
loop, the readings were very consistent. 

Secondary Loop 3. In Table 6 the closure for secondary loop 3 is 
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TABLE 4. SECONDARY Loop 1 


Initia] Station: C6 
Terminal Station: C5 


Primary Loop: 1 step 
Secondary Loop: 8 steps 


Station 
Item Method -~ Readings Aver- | Maxi-; Mini- 
Step age | mum | mum 
Primary | Single C6 35-2 | 35-2 | 35-2 | 34-7] 34-92 | 35-2 | 34-7 
Loop Reading C5 40.7 | 40.8 | 40.2 | 40.2 | 40.48 | 40.8 | 40.2 
Co-5 — | +5.56 |+6.1 |+5.0 
Oscillating | C6-5 — | +5.07} — 
Secondary] Single C6 34-4 | 34.6 | 34.6] — 34-53 | 34-6 | 34-4 
Loop Reading Cs 38.2 | 38.1 | 38.1] — 38.13 | 38.2 | 38.1 
C6-5 — | +3.60 |+3.8 |+3.5 
Oscillating | C6-5 — | +4.28} — 
Closure | Single C6-5 —1.96 |—2.6 |—1.2 
Reading 
TABLE 5. SECONDARY Loop 2 
Initia] Station: C5 Primary Loop: 1 Step 
Terminal Station: C4 Secondary Loop: 8 Steps 
Saieal Method Val Primary Secondary Cl 
erva e alue Loop Loop osure 
Single Minimum — 20.9 —20.4 +0.1 
Reading Maximum —21.7 —20.8 +1.3 
Average — 21.23 —20.57 +0.66 
Oscillating — 21.43 — 22.08 —0.65 
TABLE 6. SECONDARY Loop 3 
Initial Station: C4 Primary Loop: 1 Step 
Terminal Station: C3 Secondary Loop: 8 Steps 
Loop 
Interval Method Value Closure 
Primary Secondary 
Single Minimum +14.8 +15.5 —o.1 
Reading Maximum +15.8 +15.7 +0.9 
C4-3 
Average +15.13 +15.63 +0.50_ 
Oscillating +15.28 +15.92 +0.64 
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shown to vary from —o.1 to +0.9, with an average of +0.50 for 
single readings, and to have a value of +0.64 by the oscillating 
method. 

Secondary Loop 4. In secondary loop 4, the initial station C6 is on 
primary loop C and the terminal station B2 is on primary loop B. 
These must be connected through the base Bo=Co. The summary is 
shown in Table 7. The closure by the single reading method has values 
from —4.4 to —8.8 with an average of —6.18. The oscillating method 
shows a closure of — 5.21. Accordingly, there is a marked discrepancy 
between the step as determined on the secondary and on the primary 


TABLE 7. SECONDARY LOOP 4 


Initial Station: C6 on P Loop C, 
Terminal scenes B2 Secondary Loop: ro steps 
Loop 
Interval Method Value Closure 
Primary Secondary 
Single Minimum —2.4 —10.8 —4.4 
Reading Maximum —6.4 —II.2 —8.8 
C6-B2 
Average —4.79 —10.97 —6.18 
Oscillating —5.15 —10.36 —5.21 


loops. Whether the discrepancy is due to instrumental irregularity or 
to an operational mistake cannot be determined from the data. As 
loop C indicated an error of two whole divisions, it is possible that 
the instrument was set up over the wrong location, or that there was a 
sudden unexplained shift in the instrument. In any event, primary 
loop C should be checked by further observations. 


CALIBRATION 


In calibrating an instrument over a known gravity difference, vari- 
ous methods are used. One method is to compute the difference for 
consecutive readings and to average the results. If the instrument has 
a monotonic shift, this results in a reading which may be too high or 
too low if the readings are paired before subtraction. If the differences 
are taken for each consecutive pair of readings, the result may be in 
error. If the readings at each station are averaged and the difference 
then calculated, the error will usually be small. In Table 8 and Figure 
6 are shown the data for a calibration loop. A distinct drift downward 
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is noticeable. The results for various methods of computing the differ- 
ence are shown in Table g. If the readings for station 1 are reduced by 
204.1, the results shown by crosses in Figure 6 combine with the read- 
ings for station 2 to form a drift curve for the loop. This curve indicates 
a low reading at 11:13 and a high reading at 1:20, with lesser ir- 
regularities elsewhere. Starting with pairs at 10:40 results in a larger 
difference than starting with 10:57. The average of all pairs falls 
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Fic. 6. Observation on calibration loop. 


TABLE 8. CALIBRATION Loop 


Station 1 Station 2 

Time Reading Time Reading 
10:40 349-3 10:57 144.3 
347.2 11:29 143-3 
11°45 346.5 12:00 141.8 
12:16 345-9 12231 FAT 
12:47 345-7 1:03 140.8 
1:20 345.6 1:39 140.7 
1:56 344.3 2212 140.1 
2332 343.9 2:48 140.2 
3°05 343-9 3°23 139-5 
3:38 343-5 3°54 139.1 
Av. 345.58 Av. 141.15 
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reasonably close to the value obtained by the oscillating method, as 
does the average of the pairing with the even and odd starts. 

After the gravimeter was provided with a better temperature con- 
trol, a calibration loop with four readings at each of two stations pro- 
duced oscillating differences varying from 164.65 to 164.75 with a 
mean value of 164.69 and a root-mean-square deviation of 0.03 divi- 
sion, corresponding to about 2 microgals over a time interval of about 
an hour. The direct readings varied by 0.2 at one station and 0.3 at 


TABLE 9. COMPARISON OF CALIBRATION METHODS 


Single Reading Method 

Successive pairs 

From From All pairs | Average Method 

10:40 10:57 

reading reading 

Minimum difference 203.7 202.9 202.9 _ 203.40 
Maximum difference 205.0 204.8 205.0 _ 204.85 
Difference range 1.9 1.45 
Average difference 204.43 203.79 204.13 204.43 204.09 
RMS of deviations 0.42 0.52 0.57 — 0.38 


the other, so that single readings might lead to values 0.5 divisions 
apart. The difference of the average readings was 164.66 and the 
average of all consecutive pairs was 164.69. 


PLANNING A SURVEY 


If the variations in a survey are small enough to justify or require 
the reduced speed incidental to multiple readings, the four-point oscil- 
lating method is recommended, as it leads to more accurate values 
than arithmetic averaging and as it helps to detect gross errors and 
sometimes permits improvement by smoothing readings beyond the 
least count of the instrument. 

Such a survey is indicated in Fig. 7. A base is selected near the cen- 
ter of the area to be surveyed. Primary loops, indicated by capital 
letters in Fig. 7, are selected to cover the area, each closing at the base. 
For small areas or preliminary surveys, these primary stations are con- 
nected by chains adjusted uniformly for closure on the primary net. 
For larger areas or closer surveying, a secondary net is planned to con- 
nect primary stations by the oscillating method. The closure on the 
primary net is distributed uniformly over each loop of the secondary 
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net. Finally chains are used to connect points, each on one of the two 
nets and adjusted uniformly for closure on these nets. In Fig. 7 the 
secondary loops are indicated by numbers, the last digit of which 
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Fic. 7. A survey plan. Three primary loops, A, B, C, are shown by heavy dashed 
lines and solid circles; secondary loops by dotted lines and open circles. 


indicates the secondary station on the loop, and the preceding figures 
indicating the number of the loop. 


SUMMARY 


This paper has discussed a four-point oscillating method of ob- 
serving with an instrument subject to an appreciable zero-creep or 
drift. Examples of the method and comparisons with other methods of 
interpreting the data are given. For instruments reliable within the 
permissible least count, single readings are recommended as being 


A 
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faster. However, the present method is recommended when the in- 
strument is reliable except for the drift. 


APPENDIX 


For those interested in the mathematical aspects of this method, 
the following analysis is included: 

Since time measurements involve only intervals, the initial time is 
arbitrary and the constant term in the drift function cannot affect the 
final results. Since the method is equivalent to determining vertical 
distances between curves, it is apparent that for linear drifts, the 
vertical distances are constant and hence the linear term of the drift 


TABLE 10. ALGEBRAIC REDUCTION 


Step Difference 
Station Time Reading 
Individual Mean 

a ty a 

b te by ay f D 1 

a ts a2 Do 

b ty bo D3 

a ts a3 Ds 

b te bs 


function should not affect the final results. The quadratic term of the 
drift function does not disappear from the final results for the general 
case, but for the technique of observations selected for this method 
and for readings made at equal intervals, the contribution of this 
term is nil. For conditions approximating these conditions the con- 
tribution will be small and for actual surveying, the quadratic drift 
may be neglected. The cubic term of the drift function contributes to 
the final result but because of balancing terms, this effect is also 
negligible for actual surveying. 

Quantitatively, these conclusions are justified by the following 
analysis. Consider a set of readings as shown and reduced in Table 10. 

For linear interpolation and arbitrary reading times, 


ay =a+ (a2 Di, 
tz — ty 
, ts , 
by = 6b, + (b2 — dy) Dz = be! — ae 
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tg — bs 
= a2 + (a3 a) Ds; = bz — 
5 
bs = be + (bs — be) = — a3 
6 — la 


Do + Dz + Ds; + D,). 
The drift error is E where 
4E = 4|Do (bo ao) | =Di+De+D3+ Ds 4(bo ao). 


where do and bo are the true values at the two ends of the link referred 
to a common reference value. 
If the drift function is 


S(t) = qo + git + gol? + gat? + --- 
and we write f/, =f(t,), we may write, for the technique discussed above 
4 =a+ fi bi = bo + fo 
a2 = ao + fs bo = bo + fa 
a3 = ao + fs bs = bo + fo. 


Then 

By = Di (by ~ = + fo 
t3 — le 

E, = Dz — (bo — a) = (fs — fe) 
ta — bs 

Es; = Ds — (bo — ao) = + (fa — fs) — ier (fs — fs) 
ts 

Ex = Ds — — ao) = Us fi) 


which may be formulated as 
trai — ft 
— 
Thus 
E = + Es + E,). 


i 
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For simplicity, it is convenient to calculate the contribution to E 
of each term of the drift function separately. This is possible because 
fi; occurs only in the first order in E. 

Since the drift function occurs in E; only in differences of pairs of 
values, the go term makes no contribution to E. 

The contribution of git in to E; is 


(<2 Ee 


— be 
— te 
so that the linear term of the drift function makes no contribution 
to E. 
The contribution of got? in f(t) to E; is 
— 


ty 


— 
= (— — te?) — (tear — te) (tere + 
= (— — tepitere + — tates] 
Assigning values 1, 2, 3, 4 to k and adding, the contribution of the 
quadratic drift to E is 


Q= = {ti(ts — te) — — te) + ta(ts — ts) 


By a similar analysis, the contribution to E of the cubic term in the 
drift, gst*, is found to be 


4 \+(ts—ts)(ts— te) (tr 


If we assume, as a reasonable approximation to standardized oper- 
ating technique, that the operating times have the constant value « 
and that the travel time between adjacent stations is v, the same 
for the entire set of readings, the contributions Q and C may be 
evaluated. The time # includes the complete operating time from ar- 
rival at to departure from the station. For this assumption, inspection 
of Fig. 1 leads to observation times shown in Table 11. 

For these observation times, Q vanishes and 


II 
C= (~ ue + = s) for the initial and 
4 2 


2 
‘ 
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terminal links and 


I 
C = v® for intermediate links. 
2 


As a numerical example of the magnitudes involved, assume a 
pure cubic drift of five divisions in three hours. Then 


= qst® 
where 
: 
(180)8 


If the operating time is «= 5 minutes per station and the travel time 


TABLE 11. OBSERVATION TIMES 


Link 
Initial Intermediate Terminal 
Time . 


° ° 
ty u+ v v u— v 
ts 2u-+ 20 3u+4 3u+4v 
bs 3u+ 30 sv sv 
ts 5u-+6v 6u+8v 5u+6v 
te 6u-+7v 6u-+ 70 


is ¥=3 minutes between adjacent stations, the drift error for the end 
links is 


3 3 
55 (=) + = (=) = 0.0004 scale divisions 
4 2 \180 


and for intermediate links is 0.0006 scale divisions. 


DISCUSSION 


The foregoing paper should be of interest to geophysical field 
operators and technical supervisors, especially in gravimeter work. In 
this connection, some supplemental considerations may be of sufficient 
interest to warrant their publication. 

Dr. Roman’s paper recommends four successive readings at each 
observation point in a line of stations to eliminate the influence of 
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instrumental drift. Accepted gravimeter field practice!:?*-* seems to 
have been to make only three successive readings at each observation 
point. This naturally suggests a question about the optimum number 
of successive readings for a pre-assigned precision. It is easily shown 
(assuming linear drift between two successive readings at one point, 
equal precision of all readings, and equal time intervals between all 
readings) that the probable error of the determination of the value 
(e.g., gravity difference) of a single “‘leg”’ is 


> (1) 


where c is the probable error of a single reading and 1 is the total num- 
ber of successive readings at the two endpoints of the leg. However, 
the minimum number of readings which will give a single determina- 
tion of a leg is three (two readings at one station to establish the in- 
strumental ‘‘zero reading” and drift and one intermediate reading at 
the other station). The probable error of the observed difference be- 
tween the two points from such a set of three readings (calculated 
directly from the Theory of Errors, since Eq. (1) does not holdforn< 4) 
may be defined as the probable error of a ‘‘single observation” and is 
given by 


3/2. (2) 


Introducing this precision unit, Eq. (1) becomes 


8n — 22 
na (n 2 4). (3) 


Therefore, solving Eq. (3) for m, the necessary number of successive 
readings per leg can be pre-assigned for any desired precision, if the 
precision factor (ro) of the instrument is known. 


1L. L. Nettleton, “Geophysical Prospecting for Oil,”” McGraw-Hill, 1940, p. 38. 

2 S, Hammer, “A Gravity Profile Across the Central Appalachians . . . ,”” American 
Geophysical Union Transactions of 1941, Part IT, p. 353. 

3 E. V. McCollum and Andrew Brown, “Use of the Gravity Meter in Establishment 
of Gravity Bench Marks,” Gropuysics, VITI, 1943, p. 381. 

4G. P. Woollard, “Transcontinental Gravitational and Magnetic Profile of North 
America and Its Relation to Geologic Structure,” Bulletin of the Geological Society of 


America 54, 1943, P- 752. 
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Conversely, this formula can be used to estimate the instrumental 
precision. Applying Eq. (3) in the Law of Propagation of Errors, we 
find that the expected (“Normal”) probable error of closure of a closed 
loop consisting of m legs is 


= To | (n = 4) (4) 


where n, is the total number of readings at both the two endpoints for 
any one of the legs, each leg being observed independently. For the 
less general but more common case in which all of the m; are equal, 
Eq. (4) becomes 


To 8n 
(n 2 4). (4a) 


3 


The formula for the probable error of closure when the whole closed 
loop of m legs is observed by a single continuous looping operation 
with n’ readings on each of the stations (with a correspondingly 
smaller number of readings on the two endpoints) is 


r + (8n' — 15)m 
(n’ 2 3). (5) 


Therefore, the instrumental precision factor (7o) can be estimated 
from the observed errors of closure with either type of operation by 
applying the appropriate formula (4), (4a), or (5). 

Knowing 7o, formulas 4, 4a, and 5 also provide bases for assigning 
“tolerable” loop closure errors which may be used to decide whether 
the observations are normal or subnormal, requiring repetition. For 
example, if the instrument performance is “normal,”’ the chance is 
only one in six that the observed error of closure exceeds twice the 
probable error of closure calculated by these formulas. 

The relative efficiency of the two operating methods represented 
by Eqs. (4a) and (5) may be judged in the following manner. Equating 
the two formulas, we can calculate the relationship between the num- 
ber of readings m and n’ for equal precision with the two methods. 
However, Eq. (4a) represents a total number of readings in the loop 
of m legs 


N = nm (6) 


i 
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and Eq. (5) similarly represents a total number of readings 

N’ = n'm+1. (7) 
Using the proper values of ” and n’, the ratio of Eqs. (6) and (7) gives 
the relative number of total readings by the two methods to yield 
equal precision for a given number of legs. The results are indicated in 
the following table: 


TABLE I. RELATIVE NUMBER OF TOTAL READINGS FOR EQUAL PRECISION 


m 4 6 8 10 15 20 


N/N' for 
N/N’ for 
N/N' for 
n'=7 E 1.078 48.100 4.136: 2.344. 


It is evident that the continuous looping method* requires fewer net 
setups per leg for a given precision and therefore, in the sense implied 
here, is the more efficient. 

Another idea which is brought to mind by Dr. Roman’s paper is the 
relative rate of progress along a line observed by ‘“‘loops”’ or by “‘loop 
and fill-in,” in which the former designation means that all stations 
in a line are observed by looping back and forth and the latter means 
that the line of subbases only is observed by looping, and the inter- 
mediate stations between the subbases are observed only once on a 
single return trip along the line. (See references 2 and 4 above.) The 
ratio of the total times required to observe a line of stations by these 


ort. I ) to + d 
m+1 N+1 : 


* The “continuous looping” or more simply “looping” method is similar to Roman’s 
“oscillating”? method, but with a minor difference in the sequence of observations. If 
successive stations are designated by the letters a, b, c,d, e, - - - , then the sequence of 
station observations in Roman’s “oscillating method” is 

ababachchdcdce - 
while the sequence of observations in the continuous looping method (for three setups 
per station, i.e., n’=3 in Eq. (5) above), considered in this analysis, is 
ababchedcde 


(8) 
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for the case of three readings per station in the looping operation, 
where 

T, =time required with the “loop” method 

T;r=time required with the “loop and fill-in” method 

N =Total number of stations in the line 


n =Number of subbases in the line 

m =Number of intermediate stations per subbase interval (assumed 
uniform) 

d =distance between stations (assumed uniform) 

V =average speed of travel between stations (assumed uniform) 

fo =time required for a reading (length of stop-assumed uniform) 


Plotting curves from formula (8) with various values of the parameters 
shows that the loop and fill-in method is the faster for closely-spaced 
stations, high speed of travel, and many intermediate stations per 
subbase interval, but that there are conditions under which the loop 
method would be faster. 

There is, of course, a difference in precision in the two methods 
which should be considered in choosing the operating method. With 
the loop method, all stations have equal and relatively high weights; 
whereas. with the loop and fill-in method, the intermediate stations 
have lower weight than the subbases. 

SIGMUND HAMMER 
Gulf Research and Development Co. 
Pittsburgh, Pa. 


AUTHOR’S REPLY 


The discussion by Dr. Hammer should prove a valuable adjunct to 
the present paper in evaluating the efficiency of the oscillating method. 
Reviewing it leads to several comments which may assist the reader: 

1. In the paper the objectives were to present the method in suffi- 
cient detail to permit its being adopted by other operators and to 
give an indication of the magnitude of the error introduced into the 
final calculation by the drift of the zero point of the instrument. 

2. All readings were assumed to be exact. It was recognized that 
the drift error represents only one source of error, but the overall error 
was not considered pertinent to the main objectives. 

3. Although valid in large populations, many of the formulas of 
sampling statistics are frequently misleading when applied to popula- 
tions of small extent, such as occur in geophysical surveys. They are 
especially likely to be misleading when applied by readers not trained 
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in their basic assumptions in the same manner as reduction of data 
by least squares will sometimes lead to less satisfactory values than 
graphical analysis. 

4. One important advantage of the oscillating method is its ability 
to detect abnormal errors in otherwise good data. By taking four read- 
ings instead of three, as previously, four differences are obtained in- 
stead of two. If one difference is outstandingly different from the other 
three, an error may be suspected and the reading omitted or adjusted. 
The drift curve shows this discrepancy at a glance. In some instances, 
the nature of the error may be deduced and the data corrected. 

5. The value of the method increases as the maximum drift in- 
creases with respect to the reliability of the readings, especially for 
non-linear drifts. If the maximum drift is small, chain readings are 
quicker and sufficiently reliable. For linear drift, a uniform distribu- 
tion of the closure error will prove useful and reliable. However, the 
chain and one-step oscillating method provide no information on the 
drift behavior as the closure error may be due to many causes and is 
not always due to linear drifts. 

6. In measuring quantities which are subject to large fluctuations, 
the oscillating method may be misleading. It leads to an accurate 
measure of the differences between the end-point values at the times 
of reading but such precision may be meaningless. In gravimetric 
surveys, the gravity values may be assumed to stay constant over 
long intervals so that a specific gravity value may be associated with 
each station. Hence, extra precision is frequently warranted. In mag- 
netic surveys this may not be the case. Since the earth’s magnetic field 
is continuously changing, a highly precise determination of the abso- 
lute value of the vertical or horizontal component at a given point, as 
a reference value, has little meaning since such a value is valid only at 
the moment of measurement. On the other hand, a mean value of a 
series of such measurements, similar to mean sea-level for ocean ele- 
vations, would be of interest. However, if a survey is made for the pur- 
pose of detecting local anomalies, precise measurement of magnetic 
differences may prove highly desirable. The present oscillating method 
can remove errors which affect the entire survey uniformly, provided 
the closure error is sufficiently small. If the closure error is large, a 
local sudden disturbance is usually to be suspected. 

7. If the errors of individual readings are comparable to the errors 
due to drift, the oscillating method furnishes no advantages over 
chain readings with several independent readings at each station. Dr. 
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Hammer’s discussion does not indicate whether the analysis is based 
on the latter scheme or considers the distribution in time of the read- 
ings as well as the number of readings at each station. 

8. A primary criterion as to the desirability of increasing the pre- 
cision of determination is the magnitude of the error relative to the 
range of the measured quantity. Thus, for a range of 20,000 mi¢rogals 
there would usually be no advantage in trying to make.a determina- 
tion to 10 microgals. On the other hand, in an area for which the range 
of readings is 200 microgals, an attempt to determine gravity differ- 
ences or values to a few microgals may lead to detection of an anomaly 
if the readings are sufficiently uniform and if such an accuracy is at- 
tainable by repeating observations. 

9. Under ordinary conditions, the oscillating method corrects for 
the drift, whether due to internal or regional changes, surprisingly 
well. The numerical example verifies this conclusion. However, it 
makes no pretense at correcting for individual reading errors except 
when one or a few readings have values obviously not in agreement 
with the other readings. The method can detect such discrepancies but 
the interpreter must assume the responsibility for adjusting or dis- 
carding readings and for requesting any necessary repeat measure- 
ments. 

IRWIN ROMAN 
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A RESISTOLOG SURVEY OF THE LOMA ALTO- 
SEVEN SISTERS AREA OF McMULLEN 
AND DUVAL COUNTIES, TEXAS* 


THOMAS S. WEST{ anp CLARENCE C. BEACHAM} 


ABSTRACT 


Additional Resistolog field data are shown. This survey is in the Loma Alto-Seven 
Sisters area of McMullen and Duval Counties, Texas. The subsurface geology of this 
area is also shown along with an electric log cross section to which Resistologs have been 
added. Four of the apparent resistivity curves employed for calculating Resistologs are 
included for demonstrating the relatively great extent to which apparent resistivity 
may be influenced by shallow inhomogeneities along a traverse of electrodes which are 
moved for securing a variation in electrode separation. Several cases of direct detection 
of oil and gas saturation and successful structural correlations are presented. 


Additional field results obtained with the resistivity Increment 
or Resistolog procedure are presented. The area covered by this sur- 
vey is in the Jackson trend of South Texas. It includes the northern 
edge of the Seven Sisters Field in Duval County, the Loma Alto 
Field in‘McMullen County and the intervening area. These data are 
not, however, more recent than those previously described but repre- 
sent material which was omitted from the previous paper because of 
space requirements. Although the method of direct current resistivity 
prospecting has been previously described! a brief review of the pro- 
cedure should be of interest. 

The Resistolog method is characterized by the elimination of the 
influence of shallow inhomogenities from the final data. This has 
made possible the determination of the various subsurface electrical 
anomalies by direct measurement with a degree of accuracy which, 
under average field conditions, seems to be limited only by the ac- 
curacy with which electrical measurements may be made. Various 
special tests by which features were proved to be of subsurface origin 
have also been described. 

Unfortunately this procedure is still in the research stage insofar 


* Presented at the Houston, Texas, meeting, 1945. 

¢ Blanco Oil Company and Al Buchanan, San Antonio, Texas. 

t Somerset, Ohio. 

1 Thomas S. West and Clarence C. Beacham, Precise Measurement of Deep Elec- 
trical Anomalies, Gropuysics, Vol. IX, No. 3, pp. 494-539 (Oct. 1944). 
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as equipment is concerned, except for use on a semi-experimental 
basis. The required accuracy of measurement of potential difference is 
barely within the scope of the best measurement procedure which 
has thus far been employed and then only when in the hands of an 
operator having unusual skill and experience. Only with the exact 
measurement procedure and equipment previously described has it 
been possible to obtain Resistologs which are positively correlatable 
with relatively large electrode separation. Positive correlation of 
Resistologs was not possible with a number of seemingly equivalent 
procedures, apparently because the maximum attainable accuracy 
of measurement of potential difference under field conditions was 


CALIBRATED 


AUXILIARY 


POT. 


Fic. 1. Schematic diagram of circuit employed. 


slightly but critically reduced. The primary source of difficulty is not, 
however, due to lack of sensitivity of electrical instruments but to 
the disturbing influence of natural earth current fluctuations. It is 
believed that this difficulty may ultimately be overcome by the de- 
velopment of photographic recording equipment suitable for obtain- 
ing, in the presence of earth current fluctuations, the very high de- 
gree of accuracy of potential difference measurement required. The 
procedure employed for obtaining these and previous data! makes 
possible a measurement of potential difference which is independent 
of natural earth potentials, without the use of commutating devices. 
The possibility is thereby avoided of eliminating any electrochemical 
potential which might require a relatively long time to reach maximum 
value. Fluctuations of natural earth potentials do, however, increase 
the error of observation making necessary the averaging of a large 
number of values for securing the required accuracy of measurement. 

Fig. 1 shows a schematic diagram of the electrical circuits em- 
ployed. An auxiliary potentiometer is used for bringing galvanometer 
deflection to zero when no current is flowing in the calibrated poten- 
tiometer working current circuit and to the current electrodes. The 
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Fic. 2. Illustration of the resistivity increment (Resistolog) method to eliminate 
the effect of shallow inhomogeneities. 


latter two circuits are then closed and galvanometer deflection again 
brought to zero by adjusting the calibrated potentiometer. If the 
galvanometer deflection is still zero when the calibrated potentio- 
meter working current circuit and circuit to the current electrodes 
is opened, the calibrated potentiometer setting is recorded. Otherwise 
this cycle of operation is repeated until a stable zero galvanometer 
deflection is obtained for this period of time. 
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Fig. 2 illustrates the influence of shallow inhomogeneities on elec- 
trical resistivity measurements and the method employed for reduc- 
ing their influence to a negligible value. In Fig. 2-A, the influence of a 
shallow inhomogeneity along the potential electrode traverse is illus- 
trated. Direct current is passed between electrodes 1 and 3. The 
potential difference resulting from this current flow is measured be- 
tween two potential electrodes at E. Electrodes 1 and 3 remain sta- 
tionary for a given station while measurement of potential difference 
is made for a series of positions of the potential electrodes. The mag- 
nitude of current flow, the resulting potential difference, and separa- 
tion between electrodes provides data from which an apparent re- 
sistivity value may be calculated for each position of the potential 
electrodes. A curve may be plotted from these data showing the varia- 
tion of resistivity with electrode separation, ‘‘a.” If for otherwise 
homogeneous conditions, potential electrodes are moved across a 
shallow inhomogeneity a corresponding feature is produced on the 
resulting apparent resistivity curve. 

Fig. 2-B shows the apparent resistivity curve which results if 
potential electrodes remain stationary on a shallow inhomogeneity, 
while electrode separation is varied by moving current electrode 1. 
This curve is approximately a horizontal straight line except for rela- 
tively small electrode separation. 

The Resistolog procedure? is illustrated in Fig. 2-C. An additional 
current electrode, 2, is placed in the vicinity of electrode, 1, preferably 
along the line determined by 1 and the potential electrodes. For each 
position of the potential electrodes two measurements of potential 
difference are made, one with current flow between 1 and 3 and an- 
other for current flow between 2 and 3. These two potential differ- 
ence values along with corresponding values of current flow and elec- 
trode separations provide data from which two apparent resistivity 
curves, R; and R2, may be plotted. For electrode separation “a” of 
interest in deep prospecting, switching from electrode 1 to 2 results 
in varying electrode separation by the distance between these elec- 
trodes along the flat portion of the curve of Fig. 2-B. Obviously no ap- 
preciable change in the influence of any inhomogeneity which might be 
present in the vicinity of potential electrodes resultsfrom this operation. — 
A ratio is then taken between the two resistivity values for each posi- 
tion of the potential electrodes. This provides data for a resistivity in- 
crement curve R,/R2 which is approximately independent of shallow 


2 U.S. Patents No. 2,231,048 and No. 2,217,780. 
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inhomogeneities along the potential electrode traverse. This curve 
will not however, be independent of the influence of the horizontal 
subsurface beds lying below the current electrodes since such beds 
will produce an influence on apparent resistivity values which varies 
with electrode separation ‘“‘a.” The resistivity increment curve R,/R2 
thus indicates the rate of variation of apparent resistivity with elec- 
trode separation. This resistivity increment curve or Resistolog, will 
therefore have maxima or minima which correspond to inflection 
points on apparent resistivity curves. Correlation and interpretation 
of Resistolog are based entirely on these maxima and minima. 

Electrode separations employed in the Loma Alto-Seven Sisters 
area which will be described subsequently are as follows: Between 
potential electrodes and between current electrodes 1 and 2, 600 ft.; 
between electrodes 1 and 3, from 3,000 ft. to 4,000 ft. The distance 
‘‘a”” was, of course, variable. The increment by which potential elec- 
trodes were moved was 75 ft. in all cases. The location of station em- 
ployed is the point midway between electrodes 1 and 2. As discussed 
in our previous paper,! this point is established by reversed traverses 
across subsurface boundaries known to overlie the subsurface area 
prospected. The average elevation of electrodes 1 and 2 is used as the 
elevation of the station. i 

The formula for calculation of apparent resistivity may be de- 
rived from the differential form of Ohm’s law. The derivation of 
this formula has been described in numerous publications.* For the 
electrode arrangement employed with the Resistolog procedure the 
formula for calculation of apparent resistivity is: 


R, = 2E/I[1/a — 1/(a + 


in which J is current in amperes; E is potential difference in volts; a 
is the designated variable distance; 6 is the distance between potential 
electrodes; c is the distance between electrode 1 and 2; and, d is the 
distance between electrodes 2 and 3. 

Figs. 3, 4, 5, 6, and 7 are photographs of some of the equipment em- 
ployed in this survey. 

Figs. 8 and g illustrate the extent to which apparent resistivity 
curves may be affected by shallow inhomogeneities even under ap- 
parently uniform surface conditions and with relatively large elec- 


3D. O. Ehrenburg and R. J. Watson, Mathematical Theory of Electrical Flow in 
Stratified Media, etc., A.I.M.E., T.P., No. 400, 1931. 
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Fic. 3 Fic. 4 
Fic. 3. Inside view of truck employed for potential difference measurement 
showing galvanometer and lamp and scale reading device. 


Fic. 4. Non-polarizable porous pot type potential electrode employed 
and portable potential electrode line reel. 


Fic. 6 FIG. 7 


Fic. 5. Current electrode Loma Alto Field. Two or three 3” iron pipe in drilled and 
mudded holes 8’ to 10’ in depth were connected in parallel for each current electrode. 
Where rock occurred near the surface pipe or cable was buried in shallow trenches. 


Fic. 6. Power truck panel showing carbon rheostat employed for regulating cur- 
rent flow, magnetic switch and sensitive relay by which current was switched on and 
off by remote control from the potential truck, telephone for communication between 
trucks, miscellaneous switches and meters and ammeter light. Ammeter for current 
electrode circuit not in place when photographed. 


Fic. 7. Inside view of power truck. Storage batteries which were used 
as a current source are below elevated floor. 
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Fic. 8. Resistolog at station No. 45 with corresponding apparent resisitivity 
values plotted against potential electrode position. 
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Fic. 9. Resistolog at station No. 22 with corresponding apparent resistivity 
values plotted against potential electrode position. 
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trode spreads. No rock outcrops, gravel deposits, or other obvious in- 
homogeneities were in evidence along the potential electrode traverse 
in either of these cases. The surface soil was a very clayey loam. Yet 
the predominating features on apparent resistivity curves, R, and Re, 
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Fic. 10. Comparison of Resistologs at stations Nos. 45 and 22, also comparison of 
corresponding apparent resistivity values plotted against the separation between cur- 
rent and potential electrodes. 


very evidently occur at a given value of ‘‘a” and therefore at a given 
position of the potential electrodes rather than at a given separation 
between current and potential electrodes. Inspection of the corres- 
ponding Resistologs in each case reveals that at least the major fea- 
tures, which are obviously associated with shallow inhomogeneities 
in the vicinity of the potential electrodes, have been eliminated. The 
extent to which the apparent resistivity curves were influenced by 
shallow inhomogeneities is further emphasized by Fig. 10 in which 
apparent resistivity is plotted against separation between current 
electrodes (1 or 2) and the nearest potential electrode, rather than 
against the distance “‘a.” As is evident from inspection of the elec- 
trode arrangement, the location at which curve R,; was obtained at 
both Resistolog stations 22 and 45 differed from that of Re only in 
that electrode 2 was used as a current electrode rather than elec- 
trode 2. When plotted against electrode separation the apparent re- 
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sistivity curves at the same Resistolog station have little if any 
similarity when compared at the same electrode separation. Likewise, 
the apparent resistivity curves at station 22 are not similar as to 
maxima and minima or other features to those at station 45, even 
though the well at which station 45 was located was structurally only 
75 ft. deeper than the well at station 22. Evidently any features of sub- 
surface origin present on these apparent resistivity curves are small 
in proportion to the influence of shallow inhomogeneities. Compari- 
son of adjacent apparent resistivity curves for detecting maxima or 
minima which occur at a given electrode separation in both reveals 
that only the sequence of values from an electrode separation of 2,500 
to 2,575 at Resistolog station 45 results in any degree of similarity, 
while no such features appear to be present at Resistolog station 22. 
In contrast the Resistologs at station 45 and 22 have definite similar- 
ity. Furthermore, correlation of these Resistologs indicates a struc- 
tural difference which agrees to within 11 ft. with the true structural 
difference as determined by correlation of the logs of the wells at these 
stations. As shown by the map and cross section of Figs. 11 and 12, 
stations 22 and 45 were approximately 2} miles apart and differ in 
depth structurally by 75 ft. (i.e., without correction for difference in 
surface elevation). Such a degree of similarity of Resistologs along 
with the close agreement with well data as to structural position is 
rather definite evidence that features are associated with correlatable 
subsurface strata. This conclusion seems well supported by other 
Resistologs in the area both as to structural correlation and to the 
detection of oil or gas saturation. 

To determine the suitability of the Resistolog procedure for use 
in the Seven Sisters-Loma Alto area, a preliminary traverse was made 
across the eastern edge of the Seven Sisters Field. At the point chosen 
the water-oil contact in both the Government Wells sand at 2,500 
ft. and the Chernosky sand at 2,200 ft. were determined by drilling 
to be approximately coincident. Resistologs were taken at stations 
only 600 ft. apart on opposite sides of the water-oil contact. Not only 
did the Resistolog on the oil saturated area have features of the type 
resulting from layers of higher than normal resistivity which were 
absent in Resistologs off this area, but these features also differed in 
electrode separation by an amount corresponding to the actual dif- 
ference in depth. These Resistologs along with associated well logs 
have been previously published.’ 

Fig. 11 is a subsurface map of the area between the Seven Sisters 
and Loma Alto Fields contoured from well data on the basis of pres- 
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ent information. This shows the locations of the Resistolog stations 
of the cross section. The indicated cross section (Fig. 12) shows plotted 
well logs and, where available, the resistivity portion of electric logs. 
Resistologs have been plotted against true depth on this section by 
use of an empirically determined curve giving relation between. 
depth of prospecting and electrode separation.' Key Resistologs are 
also shown. Fig. 13 shows Resistolog correlations and comparisons 
with well logs. 

Potential electrodes were moved by 75 ft. increments in all cases, 
this survey having been made before it was discovered that significant 
additional detail could be secured by use of smaller increments. Cur- 
rent values ranged from 6 to 10 amperes. The accuracy of measure- 
ment of potential difference was such that a given group of potential 
difference values had a probable error of less than one microvolt. 
At least two such groups were used in calculating final data. 

Potential difference measurenent was made approximately 10 sec- 
onds after current flow was initiated. From the standpoint of resist- 
ance, inductance and capacitance. this should have resulted in a steady 
state measurement. It is not known whether this conclusion holds in- 
sofar as electrochemical phenomena are concerned. In fact, observa- 
tion of galvanometer deflection suggests that steady state conditions 
may not have been actually attained. 

Three separate oil saturated areas are established in the area 
covered by this-map—Seven Sisters, the designated Argo sand area 
which thus far has proved non-commercial, and Loma Alto. The pro- 
ductive area shown at Loma Alto is from a depth of approximately 
2200 ft. i.e., from the Chernosky zone. The deeper and more recently 
discovered Loma Alto-Wilcox field is about three-fourths of a mile 
north of the indicated shallow area. 

The traverse from Parr Oil Company No. 1 Sling northward was 
made before the discovery of the north edge of the Seven Sisters Field 
by drilling Parr Oil Company No. 2 Sling. The existence of the edge 
of the field at this point was indicated by comparison of Resistologs 
at stations 29 and 31 as well as by other nearby Resistologs. It will 
be noted from Figs. 12 and 13 that a high resistive feature was pres- 
ent in the Resistolog at No. 1 Sling and absent in the one at No. 2 
Sling. This feature occurred at an electrode separation having the 
same relation to actual depth as that established by the previously 
mentioned traverse across the east edge of the field. The net anom- 
aly resulting from the replacement of oil saturation by water 
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saturation in the top portion of the Government Wells sand is indi- 
cated by the shaded area of Resistolog No. 29, Fig. 13. The form of 
this difference is probably altered to some extent by a stray oil sand 
of low permeability immediately above the Government Wells in No. 
2 Sling, the anomaly attributed to this sand being indicated by the 
designated shaded area of Resistolog 31, Fig. 13. 

The indicated Argo sand area of Figs. 11 and 12 is a case in which 
a definite Resistolog anomaly results from an apparently non-com- 
mercial shaly oil sand. The thickness of this section in the only well 
drilled within this area, Grimes No. 1 Pursch, was 30 ft. This sand was 
not present in any of the surrounding wells although Schoolfield and 
Fraser No. 1 Pursch had about 8 ft. of oil saturated sandy shale in 
this zone. The resulting Resistolog anomaly from the above 30 ft. 
of shaly oil sand is as large as that obtained on 10 ft. of permeable 
oil productive sand at approximately the same depth in the Loma 
Alto Field. ; 

Resistolog 58 when compared with Resistologs 45 and 22 clearly 
shows the resulting anomaly. The shaded area of Resistolog 58, Fig. 13, 
indicates the magnitude of this anomaly when compared with Re- 
sistolog 45. Obviously these two Resistologs definitely correlate ex- 
cept that the feature corresponding to the Argo sand is absent in 45. 
The magnitude of the anomaly in Resistolog 62 at Grimes No. 1 
Pursch which is attributed to the Argo sand when compared with 
Resistolog 45 is shown by the shaded area of the former Resistolog 
Fig. 13. The similarity of Resistologs 45 and 22 is quite striking, 
these Resistologs being 2} miles apart. The indicated correlation of 
Resistologs shows 45 to be 102 ft. lower structurally than 22 as com- 
pared with an actual structural difference of 113 ft. as determined 
by correlation of the logs of the wells at which these Resistologs were 
located. The indicated structural position at Resistologs 58 and 62 
was similarly close to the actual. Since the increment by which po- 
tential electrodes were moved was 75 ft., the possible error in cor- 
relation was 75 ft. if considered solely on the basis of detail. This 
increment of electrode separation corresponds to a depth increment 
of 68 ft. according to the previously mentioned curve! showing rela- 
tion between depth and electrode separation. However, by basing 
correlations on Resistologs having the greatest degree of similarity, 
the actual error was considerably less. Resistologs which differed in 
depths by multiples of 68 ft. almost always had a higher degree of 
’ similarity than those which did not approximate this relation. 


: 

4 


502 THOMAS S. WEST AND CLARENCE C. BEACHAM 


The series of stations immediately east of Darby No. 1 Hagist 
were for determining whether the sand lens producing in the Loma 
Alto Field extended in this direction. This traverse was made before 
the drilling of Schimmel No. 4 Hagist. The Loma Alto sand was miss- 
ing in both Darby No. 1 Hagist and Texio No. 2 Atkinson. Inspec- 
tion of Resistolog 46, the location at which Schimmel No. 4 Hagist 
was subsequently drilled, shows a well defined high resistive anomaly 
in the Loma Alto zone (Figs. 12 and 13). The magnitude of this anom- 
aly by comparison to Resistolog 45 is more quantitatively illus- 
trated by the shaded area of Resistolog 46 in Fig. 13. This area 
represents the difference betweén Resistologs 45 and 46 in the Loma 
Alto sand section. The validity of this Resistolog interpretation was 
apparently verified by Schimmel No. 4 Hagist, which encountered 8 
or 10 ft. of permeable oil sand in the Loma Alto zone and was com- 
pleted for an initial oil production of 100 bbls. per day (Fig. 11 and 12). 

The anomaly in Resistolog 45, opposite the top portion of the 
Government Wells sand, apparently results from the greater resistiv- 
ity of the top member of the Government Wells section in Darby No. 1 
Hagist as compared with this sand in Texio No. 2 Atkinson. This 
sand was not cored in No. 1 Hagist, and as a result it is not known 
whether this increase in resistivity is caused by an increase in hard- 
ness or by oil or gas saturation. 

All Resistologs in the Loma Alto-Seven Sisters area are shown in 
Fig. 14 in the correlated position. The reference lines provide a means 
of determining the electrode separation ‘‘a” for any point on any 
Resistolog. Surface elevation for each Resistolog is also indicated. 

Generally the procedure in this area was to base interpretation 
principally on key Resistologs which covered somewhat greater than 
average subsurface interval. Unusual care was also exercised in taking 
data at these stations. Interpretation was based on a comparison of 
these with Resistologs at the various wells surrounding the area under 
consideration. These key Resistologs and standards of comparison are 
shown in Fig. 12. Additional Resistologs usually covering a relatively 
small subsurface interval were then taken for filling in details between 
key stations. 

With reference to the group of Resistologs within the Argo sand 
area outlined on the map of Fig. 11, the feature interpreted as indi- 
cating oil saturation in the Argo sand is definitely present in each 
case as is evident by reference to Fig. 14. The structural position in- 
dicated by these correlations, particularly if confined to Resistologs - 
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Fic. 14. Showing all Resistologs taken in the Loma Alto-Seven Sisters area in the 
correlated position. Surface elevation and reference lines for determining electrode 
separation of any Resistolog point are also shown. 
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having a high degree of similarity, agrees with that indicated by sub- 
surface data. 

The anomaly attributed to the Argo sand is definitely absent from 
all Resistologs outside the above Argo sand area except Resistolog 
46 which as previously stated was over an area which later proved to 
be productive in the Loma Alto sand. This latter sand occurs only 
about 40 ft. stratigraphically above the Argo (Fig. 12). It was, how- 
ever, possible to differentiate between these sands in interpreting 
Resistologs. This was based on the assumption that Resistologs which 
most perfectly correlated were at stations at which the depth to 
subsurface markers differed by an even multiple of the increment by 
which electrodes were moved. 

No attempt was made in this area to apply the previously pro- 
posed! procedure for distinguishing between higher than normal re- 
sistivity due to oil saturation and that due to other causes. With this 
procedure any area is interpreted as oil saturated which has all the 
following characteristics: 

(1) Higher than normal resistivity in a given stratum as deter- 
mined by correlation of Resistologs. 

(2) A down dip edge which follows a given structural contour as 
determined by Resistolog correlation, and— 

(3) which is indicated by Resistolog correlation to be associated 
with a structural condition which is favorable for oil accumulation. 

The Argo sand area and possibly the Loma Alto Field of Fig. 11 
represent cases in which steps 2 and 3 of the above procedure would — 
not be applicable since oil accumulation is controlled by lensing out 
of the sand in all directions. Normally, however, strand line accumu- 
lations do have a water-oil contact on the down dip edge of the pro- 
ducing area in which case steps 1 and 2 and to some extent 3 would 
apply. 
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ON THE USE OF ELECTROMAGNETIC WAVES 
IN GEOPHYSICAL PROSPECTING 


C. W. HORTON* 


ABSTRACT 


An approximate analysis of the behavior of electromagnetic waves in a conducting 
medium is. given. The approximations consist of replacing pairs of electrodes by dipoles 
and of using only the first order images in the case of a layered earth. It is shown that 
under typical conditions one can measure the depth of an electrical interface 6,000 feet 
deep by means of electromagnetic waves. It is further shown that even a thin layer of 
salt water or oil bearing sand at a depth of 6,000 feet gives an effect that is easily 
measureable. 

The response of the earth to a d.c. step function is analysed for the case in which the 
displacement currents are negligible. This approximation is valid for all geological 
formations. The shape of a commutated direct current is shown as a function of the 
distance of travel. 


INTRODUCTION 


There has been surprisingly little use made of electromagnetic 
waves in prospecting for oil. This seems to be due primarily to the 

lack of a theory that is capable of predicting even the order of mag- 
nitude of effects that are to be expected. A rigorous theory of the be- 
havior of a two layer earth has been developed by Riordan and Sunde! 
but their formulae are difficult to apply. In this paper approximate 
formulae are developed on the basis of the theory of images. 

It is naturally of interest to inquire if the depth of an electrical in- 
terface can be determined from electrical reflections, but a possibility 
of greater importance is that of locating oil directly without regard 
to structural conditions. This possibility exists since oil is a good in- 
sulator and usually it is found adjacent to a layer of salt water which 
is a good conductor. 

The simplest technique would be to start a d.c. pulse in the earth 
and to measure the reflection time of the successive arrivals as is 
done in ionosphere measurements. Unfortunately this is impossible 
because the conductivity of the earth produces an attenuation of the 
high frequency components and a dispersion of the Fourier spectrum 

so that the pulse becomes broadened and distorted. 


* Defense Research Laboratory, The University of Texas. 
1 Riordan and Sunde, B. S. T. J. vol. 12, pp. 162-77, 1933. 
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This leaves the possibility of steady state measurements as an 
alternative. It is possible to use only two electrodes and measure the 
input or driving point impedance of the earth. However, the use of 
four electrodes to measure the transfer impedance enables one to 
make a null measurement which is more accurate. For this reason 
an analysis is made of the four electrode arrangement of Fig. 1. The 
potential difference, V, between the electrodes 3 and 4 is determined 
as a function of the current, J, flowing in the electrodes 1 and 2. 
This is computed as a transfer impedance V/J. 

The formulae of Riordan and Sunde have been applied by West? 


I 
a a 
x 


lic. 1. The arrangement of the current and potential electrodes. 


for the special case of a homogeneous earth when « of Fig. 1 is equal 
to 2a. 
PERIODIC DIPOLE SOURCE 


In order to develop an approximate solution the current electrodes 
1 and 2 will be replaced by a dipole at O, and the potential difference, 
V, will be set equal to the product of the electrode separation, a, and 
the horizontal component of the electric field intensity at P. 

In the case of a two layer earth, the interfaces can be represented 
by an infinite set of images. Since the air may be considered as a 
perfect insulator, the strength and phase of the images in the air 
will be the same as the corresponding images in the ground. The 
electric field strength of a dipole varies as the inverse cube of the 
distance, and in a conducting medium the attenuation produces an 
exponential decrease of the field with distance. Therefore the effect 
of all but the first pair of images can be neglected. 

On the basis. of these approximations the electric field due to a 
current dipole at the surface of a two layered earth can be developed 
in a closed form. 

It is shown in books on electromagnetic theory* that the electric 
field, E, is given by 


2 West. Geopuysics, Vol. VIII, No. 2 (Apr. 1943), pp. 157-164. 
3 W. R. Smythe, Static and Dynamic Electricity, McGraw-Hill Book Co., 1939, 
p. 440. 
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— 
E = curl curl Z (1) 


=> 
where Z is the Hertzian vector. For the case of a conducting medium 
Z is the solution of 


WZ = + (2) 


Here »=permeability, e.m.u. 
o=electrical conductivity, e.s.u. 
x=dielectric conductivity, e.s.u. 
c=velocity of light, c.g.s. 
The problem of a cosinusoidal dipole source can be solved by put- 
ting 


Z = pof(r) exp (jut) 


where w is the angular frequency of the source, Po is a constant vector, 
and r is the radial vector from the source to the point of observation. 
Upon substituting Z in (1) one finds 


f(r) = (1/r) exp {— wkr/o — jonr/c} (3) 
where 

(ux/2) [(x + — 1] (4) 

n? = (ux/2)[(x 1]. (5) 


These expressions are given in Smythe.‘ 
Let 6 be the angle between the vectors po and 79 so that po- ro = Po 
cos 6. When E is evaluated for the case of a dipole, one has 


Eraa = cos 6/r?) [1 + (wr/c)(k + jn) | 


-exp [— (wkr/c) + — nr/c)] (6) 
Etan = (po sin 0/r*)[1 + (wr/c)(k + jn) + (wr/c)*(k + jn)?| 
-exp [— (wkr/c) + jo(t — nr/c)]. (7) 


These equations show that the electric vector is transmitted away 
from the origin with a velocity »=c/n. The medium has an attenua- 
tion constant of 27k/dm per unit distance. It will be shown that- for 
the earth k/n is very nearly equal to one so that the attenuation con- 
stant per unit distance becomes 27/d. 


4 Smythe, of. cit., p. 452. 
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Equations (6) and (7) can be derived from the solution of a dipole 
radiating into free space if the velocity of light, c, be replaced by the 
complex velocity c/(n+ 7k). 

If wo, (6) and (7) become 

Eyaa = (2p0/r?) cosié = Evo 


Etan = (po/r*) sin = Epo 


(8) 


which is the field due to a d.c. dipole. For convenience the fields at 
zero frequency will be denoted Ey and Ey respectively. 

It is easily shown that if a direct current, J, flows between two 
electrodes separated a distance a at the surface of the earth, then for 
small a, 


po = (9) 


where o is the conductivity of the earth. 

Before proceeding further it is desirable to illustrate these formu- 
lae with some numerical values. In sedimentary formations o will 
be of the order of 10° e.s.u. and K will be of the order of 10. One may 
take yu to be one e.m.u. For a medium having these values (4) and (5) 
may be approximated by 


n= k=~V/uo/f (10) 


for frequencies less than one megacycle. 
Therefore for sedimentary layers one finds the velocity of trans- 
mission, v, and the wave length, A, to be 


v = cv/f/uo cm/sec; = c/Vuof cm. (11) 
If o is measured in (ohm cm)~! and if y=1, then 
vy = V10°f/o cm/sec; = V'10°/of cm. (12) 


Table 1 shows the velocity and the wave length for different 
frequencies. 


TABLE 1. VELOCITY AND WAVE LENGTH OF ELECTROMAGNETIC WAVES 
o=10° €.S.U.; p~30 ohm-feet 


Frequency Wave Length Velocity 
(cycles per second) (feet) (feet per second) 

0.1 98,430 9,840 

1.0 31,120 31,120 

10 9,840 98, 430 


100 3,110 311,200 


/ 
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For media of sufficient conductivity that »=k, Equations (6) and 
(7) may be more simply written in terms of the dimensionless param- 
meter £=wkr/c=2mr/d. One has 


Eraa = + &) + exp [— & + j(wt — 8)] (13) 
Etan = + + + exp [— + j(wt (14) 


The wave length, A, is a significant measure of the penetration of 
the waves in the earth. When the wave has travelled a distance equal 
to \ the attenuation alone reduces the amplitude by a factor of 0.002. 
This is addition to the effects of spherical divergence. 


THE DIRECT EFFECT 


In order to avoid excessive dissipation in the earth one must use 
a low frequency. This means that the measurement of the potential 
can be accomplished in a fraction of one period so that it may be 
considered as an instantaneous measurement. If this measurement is 
made at ¢=o, then the real parts of (13) and (14) are 


= + cos + sin exp (— &) (15) 
Etan = Ewl[(1 + cos + + 2£) sin exp (— (16) 


If however the measurements are made at ‘=7/4, where T is the 
period, then the real parts of (13) and (14) are 


Eyaa = E,o[(t + sin — cos exp (— &) (17) 
Eun Ewl(1 + £) sin — + 2€) cos exp (— &). (18) 


Equations (17) and (18) have geen given without derivation by 
W. B. Lewis. Figure 3 of his paper shows some experimental measure- 
ments of the go° out of phase voltages which can be compared with 
the theory. 

It is obvious that if one wishes to observe electrical reflections by 
measuring potential differences along the polar axis, the measure- 
ments should be planned so that the direct effect of the dipole is mini- 
mized. A comparison of (15) and (17) shows that this is achieved by 
measuring at ¢=T/4 instead of at t=o. 

On the polar axis 8=o and the tangential component vanishes. 
If the potential electrodes are only a fraction of a wave length the 
current source, § <1 and Equation (17) becomes 


© 2 
5 William Bradley Lewis. GEopuysics, Vol. VIII, No. 1 (Jan. 1945), pp. 63-75. 


a 


510 C. W. HORTON 


exp(—£) has been set equal to 1. Now #=47°yafr?/c? and by making 
use of (g) one finds 


apf, / rc 


and the potential difference, V, between two electrodes separated a 
distance, a, is given by 


V/I = ohms. (19) 


This formula agrees very well with the experimental data given by 
Lewis in his Figures 3 and 4. These figures show that V/Z is propor- 
tional to f/r as indicated by Equation (19). The constant of propor- 
tionality is in agreement with the observations. 

One important aspect of Equation (19) is that the conductivity 
does not enter. As this field is the background against which reflec- 
tions must be measured, it is of great importance that this background — 
be independent of variations in the conductivity at different points of 


observation. 
REFLECTIONS 


The next step is the evaluation of the effect of an interface. Sup- 
pose that there is a surface layer of conductivity o; and thickness / 
and that the material below has a conductivity o2. Then a dipole of 
moment fo at the surface will have an image dipole at a depth of 2h 
whose moment is Rpo where R is the reflection coefficient (01—<«2) 
/(o1+¢2). There is an equal dipole at a height of 24 above the ground. 
The arrangement of the images is shown in Fig. 2. 

If his large, 6 is near go° and the horizontal component of the field 


intensity for the two images at the point P, Figure 2, is 
E, = 2R(p0o/8h*)[(x + &) sin — + 28) cos exp (— &) 
where* 
= 


If the potential electrodes are separated a distance, a, then the 
potential difference, V, is given by 


V/I = (a?/8mrh*o,) RF(é) (20) 


* If the interface is shallow so that the angle @ cannot be set equal to go°, both 
radial and tangential components contribute to the horizontal component of the 
electrical field. These tend to cancel each other but over the range of practical values 
this tendency is more than offset by the decrease in the /’ term in the denominator. 
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where 
F(é) = + &) sin  —.&(1 + 2€) cos exp (— &). 


If the ‘‘transfer impedance” V/J is measured as a function of fre- 
quency one obtains a curve that is proportional to F(é). This func- 
tion is plotted in Fig. 3. and shows a maximum at approximately 
£=2.6 with a Q of 0.8. This means that if fo is the frequency at which 

RB 


Lh 


Surface 


RP. 


Fic. 2. The location of the first order images of a dipole at the 
surface of a two layered earth. R=(o; —o2)/(o1 +02). 


the transfer impedance is a maximum, then the depth, 4, is given by 
h = (2.5/4m)\/109/of cm = 215\/p/fo feet 


where p is in ohm-cm. Or 


h = 1,185\/p/fo feet (21) 


if p is in ohm-feet. The value of p can be determined from electrical 
well logs or by independent measurements. 

As an example suppose the resistivity of the upper layer is 30 ohm- 
feet and reflection gives a maximum at fy=1 cycle per second. Then 
the depth of the interface as given by Eq. (21) is 6,490 feet. For these 
values the accuracy with which fp must be determined is 0.04 cps for 
each 100 feet of error in the depth that can be tolerated. 

It is necessary to make certain that the amplitude of the reflec- 


|_| 
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tion given by Equation (20) is large enough to be measured. As an 
example consider an interface 6,000 feet deep above which the re- 
sistivity is 30 ohm-feet and below which the resistivity is 60 ohm-feet. 
Suppose that the current and potential electrodes are each 1,000 feet 
apart. From Equation (20) one sees that the maximum amplitude of 
the refiection is 4.5 microvolts per ampere. By using a null method it 
is possible to measure variations of this order. 

Thus one reaches the conclusion that reflections can be observed 


“lou 


70 


Fic. 3. The dependence of the transfer impedance on # 
which is proportional to the frequency. 


from an interface 6,000 feet in depth with only a 2:1 resistivity 
contrast. 

It must be admitted, however, that the measurement of the trans- 
fer impedance as a function of frequency is time consuming. If the 
net result is simply a determination of a few depths, the results do not 
justify the method. 

Equation (20) suggests a different approach to the problem. Con- 
sider an area in which a strong electrical interface exists. After one or 
two transfer impedance curves have been measured, the value of the 
frequency fo corresponding to the maximum of F(é) can be determined. 
If V/I is measured for the single frequency fy at a large number of 
stations, the variations will be due mainly to the reflection coefficient 
R since F(é) at its maximum is relatively insensitive to variations in 
depth. 

The reflection coefficient is superior to the depth as an index for’ 
finding oil. In most oil fields the oil is surrounded by a body of salt 
water which means that the reflection coefficient varies from —1 over 
the salt water to +1 over the oil. For the sample computation given 
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above, the change in V/J on going from oil to salt water will be 27 
microvolts per ampere. This value is more than double the value of 
the first example because R was taken as one third before. 

The effect of several layers can be eliminated in the following man- 
ner. Find the value of fp corresponding to the layer of critical interest. 
Measure V/I at this frequency and at another frequency which is 
away from the maximum of F(é£) for this layer. By making a suitable 
choice of this second frequency the effects of all layers but the critical 
one are nearly the same for the two measurements so that their 
difference is peculiar to the layer of interest. 

In the discussion of the anomaly to be expected from an oil sand 
it has been tacitly assumed that the oil sand is infinitely thick. It is 
necessary to inquire to what extent this assumption is valid. 

Suppose that the resistivities of the earth above and below the oil 
sand are the same. Let / be the depth to the base of the oil sand and 
d be its thickness. Then the first order images due a dipole source of 
moment Po at the surface is a dipole of moment Roo at the depth 
2h—2d and a dipole of moment —(R— R°) po at the depth 2h. These 
may be considered as a dipole of moment R*o at depth 2h _and a 
quadrupole. This quadrupole consists of the two dipoles +Rpo and 
—Rpo at depths 2h—2d and 2h respectively. At the surface of the 
earth over the images the ratio of the electric fields due to the dipole 
to that due to the quadrupole is R’r/2d. In practice this ratio is much 
larger than one so the dipole predominates. For a porous sand filled 
with oil or water R is very nearly +1 so that R*~R and Equation (20) 
is still valid, provided it is remembered that h is the depth to the base 
of the oil sand. 


THE STEP-FUNCTION DIPOLE SOURCE 


The analysis has shown that the dissipation of energy is the pre- 
dominant factor in the propagation of electromagnetic waves in the 
earth. This means that the real currents are more important than the 
displacement currents, and, in consequence, Equation (2) may be 
approximated by 


~ 
/dt = (22) 


where a?=c?/4myo. This has the following solution corresponding to a 
dipole source: 


Z = (po/ryerfet (23) 


{ 

4 
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where erfc is the conjugate of the ¢ error function and ¢=r/2a/t, Upon 
substituting this expression for Z into Equation (1) and reducing the 
expressions to radial and tangential components one has 


Erna = E,olerfeg + cerf’é | (24) 
= Evolerfeg + + (25) 
where erf’¢ indicates the first derivative of the error function. 
— 
0.5 
Oo 
fo) S Ve? /O 


Fic. 4. The response of a homogeneous conducting medium to a step-function 
dipole source as a function of 1/# which is proportional to the time. Eyaq is the radial 
coniponent and Ean is the tangential component. 


As >, ¢-0 and erfc ¢-—1. Thus the field reduces to a dipole 
source. 

In order to visualize the variation of the field components with 
time, the bracketed expressions of Equations (24) and (25) are plotted 
against 1/¢?, which is proportional to the time, in Fig. 4. These curves 
show the slow build up that is characteristic of diffusion phenomena. 
It will be noticed that the tangential field reaches a maximum at 
¢=1. This corresponds to a time of 


t = myo(r/c)?. (26) 
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From Equation (11) one sees that the maximum travels with a 
velocity w given by 


w = 


where 2 is the velocity of a one-cycle-per-second periodic disturbance 
in the same medium. Thus the maximum has a high velocity near the 
source where the high frequency components are present and a low 
velocity in the distant regions where only the low frequencies have 
reached. 

In practical units of resistivity in ohm-feet and distance in feet, 
Equation (26) may be written (u=1) 


t = 3.05%(r?/p) X seconds. 


If a step function is reflected from the bottom of a 5,000 foot thick 
layer of resistivity 30.5 ohm-feet, the received pulse increases from 
one-half to the full peak value in 0.2 seconds. This means that the 
time of the maximum will be difficult to measure even if there are 
no disturbances. 

Ina short d.c. pulse the maximum electrical field will occur at very 
nearly the same time as the pulse is turned off. Thus the transient 
effects in the receiver caused by the ending of the pulse occur just at 
the time when the maximum of the reflected pulse must be accurately 
timed. This coincidence tends to ruin the measurement. 


COMMUTATED DIRECT CURRENT 


In actual field measurements it is more satisfactory to use com- 
mutated direct current instead of a sinusoidal source. Therefore, it 
is of interest to examine the variations in the shape of the square 
wave as it travels through the earth. 

A square wave can be resolved into a Fourier series but the series 
converges so slowly that the representation is not satisfactory for 
numerical calculations. A more satisfactory representation is ob- 
tained by considering the commutated waves as the result of a series 
of the step functions discussed in the last section. 

Thus if we define the voltage that appears at the second pair of 
electrodes at the time / after a step function of unit current at the 
primary electrodes as the indicial impedance E(#), then the voltage 
E due to a commutated unit current is 


= (— + mT/2) + (— +NT/2) (27) 


m=0 
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Fic. 5. The effect on a square wave of traversing a conducting medium. The ampli- 
tude of the tangential component of a wave generated by a square-wave dipole source 
as a function of time for various values of the parameter 47/# which, for a given fre- 
quency, is proportional to 1/r?. 


where T is the period of the commutation. If we denote by ¢» the value 
of ¢ corresponding to ¢+mT/2, then 


I 

Tut 

where é is the dimensionless parameter introduced in Equations (13) 
‘and (14). 

In order to illustrate the shape of the square wave after transmis- 


sion through the ground it is necessary to compute E(é) as a function 
of ¢/T for several values of 47/é. The series (27) converges very 


== (t+ mT/2) = (4x/€*)(m/2 + t/T) (28) 
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rapidly although the sum alternates by +1 as successive terms are 
added. The correct sum is the average value of these two limits. 

In Fig. 5, the voltage E(#) is plotted against t/T for several values 
of 47/é. The factor fo/r* is omitted so that the effect of spherical 
divergence is not allowed for. As the distance is increased, the phase 
lags more and more and the shape of the voltage variation becomes 
more and more sinusoidal. 


ELECTRIC FIELD OF AN OSCILLATING DIPOLE 
ON THE SURFACE OF A 
TWO LAYER EARTH* 


ALFRED WOLFft 


ABSTRACT 


The electric field of a low frequency oscillator placed on the surface of a two layer 
earth is determined in two special cases, namely, the case in which the conductivities of 
the two layers are nearly equal, and the case in which the lower layer is a perfect 
insulator; in the latter case, only terms of zero and first order in frequency are con- 
sidered. It is shown that, when the upper layer is sufficiently thin or is very thick, the 
mutual inductance of two wire elements on the surface of a two layer earth has the 
same value as for a homogeneous earth. In the case of an insulated layer, it is shown that 
the maximum departure of the value of mutual inductance of two collinear wire ele- 
ments from the corresponding value on a homogeneous earth is 35%. 


Most of the extensive literature of electrical prospecting inter- 
pretation is devoted to electrical resistivity work, and does not con- 
sider alternating current phenomena. The purpose of the present paper 
is to supply the solution of some problems in low frequency electrical 
exploration, without entering into evaluation of various methods of 
electrode arrangement or measurement procedure. The problem 
chosen as the most fundamental and, at the same time the simplest 
in this field, is the theoretical determination of the electric field of a 
low frequency oscillating dipole placed on the surface of a two layer 
earth. Although it is possible to obtain formulae in this case which 
are generally applicable, they are in the form of integrals over oscil- 
lating functions which do not lend themselves either to numerical 
evaluation or to any simple interpretation. The problems considered 
are therefore restricted to two special cases, namely, the case in which 
the conductivities of the two layers representing the earth differ 
only little, and the case in which the lower layer is a perfect insulator; 
in the latter case, only zero and first order terms in frequency are 
considered. In both of these problems it was found possible to develop 
explicit formulae for computation. 

The methods of dealing with the problems of oscillating electric 
fields in conductors were developed by A. Sommerfeld.'* They were 


* Read by title at the Chicago meeting of the Society, April 1946. 
+ Geophysical Research Corporation, Tulsa, Oklahoma. 

1 A. Sommerfeld, Ann. d. Phys., 28, 1909, p. 683. 

? A. Sommerfeld, Ann. d. Phys., 81, 1926, p. 1135. 
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applied in the field of telephone communcation by R. M. Foster? 
and by J. Riordan and E. D. Sunde.* Foster determined the electric 
field of a dipole on a homogeneous earth, Riordan and Sunde ob- 
tained the general formulae given in Section I of the present paper, 
but their special formulae are applicable mainly to infinite lines. A 
derivation of the general results obtained by the above authors will 
also be given in this paper by a method which neglects displacement 
currents from the beginning, and is therefore simpler to employ. 

In the following Section I, a precise mathematical formulation of 
the problem will be given, and a formula will be derived which is ap- 
plicable to the most general case of a low frequency oscillating dipole 
on the surface of a two layer earth. 

In Section II, the case of two layers of nearly equal conductivity 
is considered. The most interesting results are relative to terms of 
first order in frequency which may be called mutual inductance terms. 
It will be shown that a discontinuity in the value of the conductivity 
of the earth has no effect on mutual inductance when the plane of 
discontinuity is either very deep or very shallow. 

In Section III, the case of an insulated layer is treated to first 
order terms in frequency. The mutual inductance of two collinear 
wire elements on the surface of the layer is shown plotted as a function 
of the depth of the layer. 

In Appendix I, a certain integral is evaluated which is required in 
Section III. Appendix II deals with the evaluation of certain sums of 
Bessel functions. 

Throughout the paper, references will be given to G. N. Watson’s 
“Theory of Bessel Functions,’”’ Cambridge 1944, which will be 
quoted as Watson. 

It is assumed throughout the following that the magnetic suscepti- 
bility of the earth is zero. The electrical units employed in this paper 
are absolute electromagnetic units. It should also be noted that all 
integrals are to be extended from 0 to ~ unless otherwise specified. 

In several instances, the same symbols are used for different 
physical quantities or mathematical constants, but only where it is 
thought they will cause no confusion. 


I, FIELD EQUATIONS IN A GENERAL TWO LAYER EARTH 


Let an electric oscillator of frequency w/2m be placed on the sur- 
face of the earth, with the axis of the oscillator in the horizontal 


3 R. M. Foster, Bell Syst. Techn. Journ., 10, 1931, p. 408. 
4 J. Riordan and E. D. Sunde, Bell Syst. Techn. Journ., 12, 1933, p. 162. 
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direction. A system of rectangular coordinate axes is introduced, 
with its origin at the oscillator, the +. axis in the direction of the 
axis of the oscillator, and the +2 axis vertically downward. Physically, 
the oscillator is realized by two current electrodes on the x axis close 
to the origin of the coordinate system, the positive electrode on the 
positive side when the phase of the current is regarded as positive. 


2g 
Surrace 


\ 
r 
4 
Ex 
Ey Conpuctivity Ki 
h / Kea 
/ 
R, 
/ 
/ 
2h ¢ 
z 


Fic. 1. Diagrammatic view of dipole and coordinate system. 


Fig. 1 shows the arrangement of the electrodes. If the current in the 
dipole is Jo, and the spacing of the electrodes is /, the strength of the 
dipole is 


C = Il exp (jot) (1) 
with 
t= time 


The oscillators considered in this paper are supposed to be of a 
frequency low enough to permit neglecting the displacement currents 
in the conducting earth. This means that the assumption is made that 


w K (2) 


where k is the conductivity of the earth in absolute electromagnetic 
units, € is the dielectric constant of the earth, ¢€9 is the dielectric 
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constant of vacuum, c is the velocity of light in vacuum. Inequality 
(2) is fulfilled by an ample margin in most exploration work, except 
possibly at radio frequencies or higher. 

When the displacement current is neglected, it follows from Max- 
well’s equations for a conductor that the current density J and the 
magnetic field H can be derived from a vector potential A by means 
of the equations 


H = curl A 4nl = — 7y?A+VdivaA (3) 
where 
7? = 4njwk 
k = conductivity in abs. e.m.u. 


and the vector potential A satisfies the diffusion equation 
V?A = 7A. (4) 


Since the displacement of currents are neglected, an insulator 
has to be treated somewhat differently from a conductor. Because in 
an insulator both y?=o and J=o, it is seen from Equ. (3) that the 
additional condition must hold 


divA =o for k=o. (5) 


At boundary surfaces between two conductors, where there is a 
discontinuity in the value of conductivity, it follows from Maxwell’s 
equations that the tangential components of H and of I/y? must be 
continuous across the boundary. At a boundary between an insulator 
and a conductor, the second of these conditions becomes meaningless 
since both y” and J vanish in the insulator; the boundary conditions 
at such surfaces require only the continuity of the three components 
of the magnetic field H. 

The earth is regarded as a conductor of zero magnetic susceptibility 
with a horizontal plane of discontinuity in the value of conductivity at 
the depth h. The value of conductivity from the surface to the depth 
h(z<h) is denoted by k:, the conductivity of the lower layer (z>h) 
is ke. Similarly, the values of y in the two layers are denoted ; and 72 
respectively. 

It is required to find solutions of Equ. (4) which satisfy all bound- 
ary conditions, and are finite everywhere except at the origin of 
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the coordinate system, where the x component of the vector poten- 
tial A, becomes infinite as 
C/R (R? = x? + y? + 2°). (6) 
It will be seen that all conditions can be satisfied with 
A, =o (7) 


which will be postulated throughout. The boundary conditions at the 
plane of discontinuity (s=/) then assume the relatively simple form 
I 
A,, 0A ,/02, A:, — div A, continuous for z = h. (8) 


At the surface of the earth, the last of these conditions becomes 
meaningless since both y? and div A vanish in air which is an insulator; 
consequently, only the first three conditions have to be satisfied. 

It follows from the work of Sommerfeld? that all requirements 
can be fulfilled by the following choice of expressions for the vector 
potential components A,, A, in the various media (A, =o) 

In air (<0) 


A f (C + (ga) 


f C+ (ob) 


In the upper layer of the conducting earth (0o<z<h), where the 
value of conductivity is ky, 


A,= f [(AC/m1 + + ret? |Jo(Ar) -dd (10a) 


dy 
A,=-— + |Jo(Ar) — - (rob) 
Ox 
In the lower layer (>), where the conductivity has the value ko, 


dy 
A,=- f (11tb) 
Ox 
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with 

= N+ Mo” = (12) 
Jo is Bessel’s function of the first kind, r?=2?+y’, the y’s are functions 
of the variable of integration \, and must be determined from bound- 


ary conditions; all integrals are extended from 0 to ~. 
* In writing out Equs. (9)—(11), use is made of the two integrals 


ff 3/R (13a) 
I 

f = —e (13b) 
R 


of which the first is well known,' and the second is discussed in Som- 
merfeld’s paper.! It should be noted that derivatives of both integrals 
(13) with respect to z vanish at z=0, as follows from symmetry; this 
fact is made use of in the derivation of boundary conditions. 

There are three boundary conditions at the surface of the earth 
(s=0) corresponding to the first three conditions (8) 


C+ Wo = + fi + (14a) 
Avo = — 2) (14b) 
C+ V4 (74¢) 


and four boundary conditions at the plane z=/ separating the two 
regions of conductivities k; and k» respectively 


(AC/mi + Wid f + = vs (15a) 
(AC + war) f — = pos (15b) 
V3f + = ve (15¢) 


[(AC/ur + Wf + + — 
= 1/y22[Ws + (15d) 
with 
f = exp (— wih) 


a total of seven linear equations from which the seven unknown func- 


5 Watson, Art. 13.2 (8). 


fi 
| | 


524 ALFRED WOLF 


tions Yo may be determined. In the following, only ¥1, 
will be required. From Equs. (14a, b), (15a, b) the functions y, 2 
can be calculated, and then used to obtain ¥3—y, from the remaining 
boundary conditions. Thus 


= — A) + we) + (ur — we) (16a) 
= — be) (16b) 


W=2(- I (—+ ) 4 ] 


with 
A = (ui + + me) + (ui — — (17a) 


I pe If be M1 


The two components of current density at the surface of the earth 
are, according to Equ. (3), 


I I “(4 + fa ~~) 
x x x 
4m 0x? 


z 
4m Oxdy 0z 


since A,=o. Therefore, denoting 


I @A; 
A, = F(r) (4s +f dx- )- G(r) (18) 


it follows that the electric field at the surface of the earthis given by 


E, = — joF(r) + —G(r) = — joF + G” —- — +— (19a) 
Ox? r dr\r 


a? xy d/G’ 
E, = (19b) 
Oxdy 


NAL \y2? 
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because of the relation 
In Equs. (19) we have, according to Equs. (10), 


G(r) = f Kez (vs — ‘dd. —(20b) 


47 Ry 


Substituting Equs. (16), (17) into (20), the following expressions are 
obtained after a little calculation 


AC 


G(r) = o(Ar)-dd (2tb) 


— (Rime — Romi) f?]A 
The equivalent of the general Equs. (20), (21) were derived by 
Riordan and Sunde‘ who applied them to various problems. 

In the following Section II, Equs. (21) are much simplified by 
assuming that the difference between the conductivities of the two 
layers is much smaller than the average conductivity, and the inte- 
-grals are evaluated. In Section III, the substitution is made k2:=o0, 
and the integrals are calculated to first order in frequency. 

The electric field of a dipole on a homogeneous earth may be ob- 
tained from Equs. (19), (21) by substituting 


ki = ko = 


C fi ki — ke)(A+ m2) f? 
[ 


For details of calculation reference may be made to the paper of 
Foster ;* the result is 


r3 
G(r) = (22b) 
ankr 


from which it follows that 


E,= [x + (1 + — 3y?/r?] (23a) 


an kr? 


2CT x 1+ yr 
/ : 
/ 
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E, = (3xy/r?). (23b) 


kr? 


It is seen that E, is independent of frequency for a homogeneous 
earth. If the assumption is made 


| yr| KI 
and Equ. (22a) developed to first order in w, the result is obtained 
C 
E,= [(2a? — y?)/r?] — jwC/r. (23¢) 
an kr? 


The mutual inductance of two short insulated wires of length / and 
I’, with their ends grounded on a homogeneous earth, is seen to be 


Ll’ cos a/ri 


where a is the angle between the two wires, and 7; is the distance 
between their centers.® 

It should be noted that the series which follow mutual inductance 
terms, in developments in powers of frequency, begin with a term 
proportional to the 3/2 power of w. Hence, the mutual inductance 
calculated from such developments is of value only when it is rela- 
tively small. The same limitation applies to similar developments in 
series in the next two sections of this paper. 


II. TWO LAYERS OF NEARLY EQUAL CONDUCTIVITY 
When the values of conductivity of the two layers representing 
the earth are nearly equal, that is when 
Ak = ko ky (ke + k;)/2. 


Equ. (21) can be simplified considerably if only terms of the first 
order in Ak are required. The expressions for the functions F(r), 
G(r), and for the electric field E can be written as sums of two terms, 
one (subscript o) independent of Ak and identical with the correspond- 
ing expression for the homogeneous earth, the other term (subscript 1) 
proportional to Ak. Thus, 


F=F,+ F, G=G)+Gi E=£+sk 


Fo, Go, Eo are given by Equs. (22), (23) of the last section, in which k 
may be regarded as the average conductivity. 


6 G. A. Campbell, Bell Syst. Techn. Journ., 2, 1923. 
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There is no difficulty in deriving the expressions 


Fi(r) = — f (24a) 
C Ak 
Gi(r) = — — — dx. (24b) 
ark k 


The integrand on the right hand side of Equ. (24a) is transformed by 
multiplying both numerator and denominator by (u-A)? which gives 


at 
4mjwk k 


It is seen that both F; and G; can be expressed in terms of derivatives 
of integral (13b) and of the new integral 


U= f To(dr) «dd 


1(2* (= *) (g>0). (2s) 


2 2 2 


There is now no difficulty in showing 


{( 0? I 1 0?U 
k 02? 2 R r 


r or \r z=2h 


C Ak xy(d/1 0U 

amk k NOr\r Ordz/) 

Since U and its derivatives can be expressed by means of Bessel func- 


tions with an argument proportional to 1+j (j=+/—1), for which 
tables are available, Equ. (26) may be used for numerical computa- 


* This integral is given by Foster* without sufficient references. It may be verified 
by observing that both sides satisfy Equ. (4), and that, for z=o, Equ. (25) is correct ac- 
cording to Watson, Art. 13.6 (3). 
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tion. The subject will not be pursued any further, but explicit repre- 
sentations will be given for the case 
| yR,| KI 
where 
Ri? = 4h? + + y? 
to first order in frequency. 
Employing relations given by Watson, Arts. 3.7 and 3.71, it can 
be shown without difficulty that, to second order in y 
I I 3 


r drag R® 2 R(R+32) 


and, after a few purely algebraic transformations, 


C Ak 
En = — — (4h? + y? — 2x7)/RE 
ark k 
Ak 2h P2hRi-—r®? y(2Ri+ 2h) 
joc = | em 
k R3 R + 2h (Ri 4- 2h)? 
C Ak 3xy | Ak 2hxy 2Ri + 2h 
amk k Rk Ry (Ri+ 2h)? 


From general considerations, it is obvious that the plane of dis- 
continuity can have no effect on the electric field at the surface when 
the discontinuity is at very great depth. Equ. (27) adds the important 
result that the quadrature component of the electric field is not 
affected appreciably by discontinuities in the value of electrical con- 
ductivity very close to the surface (ko). 


III. OSCILLATING DIPOLE ON THE SURFACE OF 
AN INSULATED LAYER 


To obtain formulae applicable to the calculation of the electric 
field on the surface of an insulated conducting layer of conductivity 
k and thickness h, it is only necessary to substitute ke=o in Equs. 
(21); at the same time write & in place of k, for the conductivity of the 
layer. The resulting expressions are still quite complicated. As stated 
in the Introduction, only zero and first order terms will be considered 
in this Section. 

Now, the equation giving G(r) contains an integral which, while 


ELECTRIC FIELD OF AN OSCILLATING DIPOLE 529 


convergent for any frequency other than zero, is diveregnt for w=o. 
The reason for this condition is the fact that G(r) contains the poten- 
tial of a point charge in an insulated layer, which is well known to be 
infinite at zero frequency. Since, however, the electric field depends 
only on derivatives of G, a term independent of r can be added to 
make it possible to develop the integral in powers of w. This is done 
by substituting the factor 


Jo(Ar) Jo(Aa) 


for Jo(Ar) in the expression for G(r); the constant a may have any 
value, since, clearly, it has no influence on the derivatives of G. 

Carrying out these operations, the development of the electric field 
to first order terms in w may be written 


— jwC/r + jwC (28a) 
E, = — joC (28b) 
ark dxdy 
with 
Gi = f coth Ah[Jo(Ar) — Jo(Aa)]-dd (28c) 


G; = f [Jo(Ar) — r]-dd (28d) 


in which, moreover, a was taken equal oin Equ. (28d) since this makes 
the integral sufficiently convergent and easier to evaluate. 
Substituting 


I 
coth x = —+ x/(x? + 
x 
and making use of Watson, Art. 13.51 (5), it can be shown that 
I a 
G, = + tog —+2)> Ko(nmr/h) — 2>> | (29) 
r 
a well known formula, in which the summation. is extended to all 
values of from n=1ton=. 


The evaluation of G2(r) is more difficult, and is treated in Appendix 
IT; the result is 


2 
C 0°G 
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2h 
= r — —logr — 
3 


nar 
K,| — 
h 
— 2r >, ——_——— + const. (30) 
in which the constant is independent of 7, and in which the sums, 
as all summations in this Section, are extended from n=1 ton=. 
With the aid of formulae in Watson, Art. 3.71, it is now easy to 
derive 


I 2 


= = = > 
+ (31a) 
dG, 2 2 
+ (31b) 
and 
2h 2 


+ 2h/ner)Ki(nmr/h)| (32a) 
dG2 
“(+ =) [(8/nmr? + 
dr\r dr 


+ (4/rh)Ko(nmr/h)]. (32b) 


The electric field E may be decomposed into two components 
E=£,+£:, where £, is in phase with the applied current or C, and 
E, is the quadrature component. The following equations are then 
obtained from Equs. (28) 


y? “( I 
a 
dr? r dr\r dr 


C xy df1 dG 
Ey =) (33b) 


amk r dr\r adr 


| 
i 
i 
t 
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and, for the quadrature components, 


r dr? r dr\r ar 
ay 
Exy = — (34b) 
r adr\r dr 


The mutual inductance M of two short collinear wires of length / and 
I’ respectively, is consequently, given by 


M = Il'[1/r — 


Fig. 2 shows the expression (M—Mo)/Mo with My=/I'/r plotted as a 
function of h/r. It is seen that the mutual inductance on the surface 


1.0 2.0 3.0 40 5.0 6.0 
h/r 


—-0.3 7 


Fic. 2. Relative values of mutual inductance of two collinear wire elements 
on the surface of an insulated layer. 


of an infinitely thin layer has the same value as the mutual inductance 
on a homogeneous earth. It is also seen from Fig. 2 that the maximum 
departure of M from Mp is of the order of 35%. 


APPENDIX I 


By a simple substitution, Equ. (28d) is transformed into 


I — — 264 Jo(u) —1 


G.(r) = ar f en2eu du (35) 


(1 
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with o=h/r. This may be written 
ou |Jo(u) —1 


r sinh? u? 


By partial integration 


J 


40 u? 2u 


because* 


— Jo(u 
a- 2. 


ur 


Further, since* 


we obtain 


=1+ 2 f (cou ou — + (36) 
r 2u 


ou u? 


Next, consider the following integral A in the complex z plane 


A= “(cost oz — =\=" = + (37) 
2% 


Po(z) = Hi(2) = 


where 


H,', H;' being defined as in Watson, Art. 3.6. The contour of integra- 
tion is supposed to be on the real axis in the complex z plane, except 
for a half circle of radius ¢ above the real axis in the neighbourhood 
of o. The integrand has poles at the points jrn/o, where v is an integer, 
and tends to the limit zero exponentially at infinity above the real 
axis. In the neighbourhood of z=o we have, according to Watson, 
Art. 3.61, 3.62, 


* These two integrals are evaluated by differentiation under the integral sign of 
corresponding integrals which include the exponential exp(—amu) under the integral 
sign, and then making a—o. Use is made of Watson, Art. 13.2 (8). 


J —_ J I 
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coth oz — — = — 


Potlogz Pi B+ I 
with 
B =log2—-y+3 (y = Euler’s const. = 0.5772..) 
It follows that the semicircle ¢ contributes 
on? /6 — jrBo/3 


to the integral A. Further, since coth oz—1/oz is an odd function of z, 
and because 


Po(— 2) = — jaJo + Po(z) Pi(— 2) = — jadi + Piz). 


as follows from Watson, Art. 3.62, we obtain 


A = jr f (cots ou — + 
ou u 2u 


+ on*/6 — jrBo/3. 


This, in turn, must be equal to the sum of the residues of the integrand 
of Equ. (37) above the real axis, multiplied by 27j7. Making use of the 
relations 


Po(jx) = Ko(x) jPi(jx) = Ki(x) 


we obtain, after a little calculation, 


1\[Jo(u) —1  Jyi(u) 
(cotn ou — je + 
= Bo/3 — 20D) (log n)/n*x? — (4/3) log (x/a) 
— 26>, — >> (38) 


in which the sums are to be extended from n=1 to n=. From 
Equs. (36) and (38) the value of G2(r) is calculated. 


APPENDIX II 


The evaluation of G2(r) by means of infinite series is quite simple 
when the arguments of the Bessel functions are great, since the series 
are then rapidly convergent. For small values of the arguments rr/h, 
we note the following formulae 


| a 
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I 
2V 2 41 2 
I ny 
2y> 
I ¥9 y 
n 2 4 8 4 
ny? 
+ (41) 
647° 
with 
I 
y = Euler’s const. = 0.5772 n= 


The series (40) is obtained by differentiating (39). The series (39) and 
(41) are derived from a consideration of the integrals 


I I dz 
f (cotn -dz and (cotn 


in the complex z plane, on a contour consisting of the real axis. The 
integrals can be reduced to the infinite sums (39) and (41) and can 
also be expressed in terms of 


I I du 
J (cot u— olay) a and (cot u— 
u u u 


The substitution is made 


I 1 — en 
coth — — = lim + — 
Uu n 


m=1 


and the integration is carried out term by term. 


. 
: 
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DISCUSSION 


ABSTRACT 


Dr. Wolf's rigorous solution of the problem of a dipole oscillator on the surface of a 
homogeneous earth shows E, independent of frequency. Field measurements show Ey 
dependent on frequency in a manner agreeing with the solution of a dipole oscillator 
in a homogeneous medium. The disagreement of the measurements with the more 
rigorous theory, and the agreement with theory that neglects the air-earth boundary is 
paradoxical. 


The electromagnetic field in a homogeneous infinite medium can 
be completely described in terms of the Hertzian vector, Z which 
satisfies the equation 


Z= 4npoZ + 
Where p is the magnetic permeability, o is the electric conductivity, 
and x is the dielectric constant. All quantities are in electromagnetic 
units. 


In terms of Z, the electric intensity, Z, and the magnetic intensity 
A are: 


E = curl curl Z 
H = curl Z+ « curl Z 


When displacement currents, can be neglected the following solu- 
tion holds. It has the symmetry characteristics of a dipole oscillator 
located at the origin and oriented along the X-axis: 


Z = ie tei+t/r 
where 
and 


= 


This gives a magnetic field having only an azimuthal component. 
The electric field has radial and tangential components, but the azi- 
muthal component is everywhere zero. 

If one can neglect the stationary magnetic field in the air, this 
solution holds for a dipole oscillator on the surface of the earth. The 
vanishing of the azimuthal component of the electric field satisfies 
the condition that no current flow across the surface. 


| 
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The radial and tangential components of the electric field are 
found to be: 


E, 
Ey 


2 cos (1 + 
sin + + 


Where @ is the polar angle. 
From this it follows that on the earth’s surface: 


E, = e*[2 cos? 0(1 + ¢) — sin? 0(1 + ¢ + 


and 
E, = 3 cos @ sin + ¢ + 


A second solution derived by Wolf from formulae developed by 
Sommerfeld and Foster is presented in the foregoing paper. His equa- 
tions (23a) (23b) give : 


E,= [1 + (1 + Het sin? 
E, = 3 cos 6 sin 6e%*/r3, 


Both solutions reduce to the field of a static dipole as w approaches 
zero. The second solution indicates that Ey, is independent of fre- 
quency. The implications of this are startling: among others EZ, must 
be in phase with the dipole oscillator at all points. 

The conflict between the two solutions appears paradoxical. Wolf’s 
solution is rigorously derived from the now classic investigations of 
Sommerfeld. 

The solution offered here is derived directly from Maxwell’s equa- 
tions and can be found in essentially this form in standard text books, 
e.g.: Smythe: Static and Dynamic Electricity p. 438 et seq. 

Because nothing in the available literature indicated that E, had 
been investigated experimentally, time was taken during routine 
field work to investigate this point. The measurements! normally 
give the quadrature component of £,; by changing one of the probe 
electrodes to a position perpendicular to the line, E, was measured. 

The results of this run, together with a sketch of the electrode 
pattern, are shown in the accompanying figure. Both components 
are seen to be functions of frequency. The relative slopes agree well 
with the formulae here presented. 


1Wm. B. Lewis, Working Depths for Low Frequency Commutated Currents, 
Gropnysics, Vol. X, No. 1. (Jan. 1945), pp. 63-75. 
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Fic. 1. Variation with frequency of longitudinal and transverse components of 
electric field and quadrature components of potential per unit current measured parallel 
to and perpendicular to line of current electrodes and sketch of electrode pattern used 


for measurements. 


Several implications of the Sommerfeld solution that are at vari- 
ance with physical intuition when low frequencies are considered have 
been previously brought to the author’s notice by Dr. H. M. Evjen. 
Never before has such a clear cut example appeared as the case of 
E, being independent of frequency. It is hoped that this brief note 
will arouse sufficient interest in this problem so that some investigator 
will resolve the seeming paradox. 


Ws. BraDLey 
The Elflex Company, 
Houston, Texas 
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INDEX OF WELLS SHOT FOR VELOCITY* 


B. G. SWANT 


ABSTRACT 


The information listed on velocity surveys in 167 wells is in the nature of a supple- 
ment to the “Index of Wells Shot for Velocity” published in the October 1944 issue of 
Geopxysics. In addition to the new surveys listed, corrections and information missing 
on many of the surveys of the original Index are tabulated. 


INTRODUCTION 


At the time of publication of “Index of Wells Shot for Velocity” 
in the October 1944 issue of GEOPHYSICS, the statement was made that 
it was planned to bring the Index up to date at approximately six- 
month intervals. Upon further consideration by the Society’s officers, 
it was decided that this was too short an interval; therefore this is the 
first such revision to be published. Further revisions may be forth- 
coming at the discretion of the Society. 

In securing the information on additional surveys and corrections 
or additions to incomplete information on original listings, companies 
engaged in seismograph operations were solicited by letter, as in prep- 
aration of the original Index. The percentage of companies respond- 
ing was very high, and for such splendid cooperation the officers of 
your Society and myself wish to express our sincere appreciation. 

The first part of this revision lists the surveys which have been 
made since the original Index was assembled; also, several earlier 
wells which had been surveyed but for various reasons were not re- 
ported in 1944. In the latter portion of the revision, corrections and 
information missing on many of the surveys of the original Index are 
listed. The corrected or added information on each survey is in italics. 

New abbreviations which appear under the “Shot By” and “‘Spon- 
sored By” columns are: “K. V. I,” for Key Velocities, Inc., and 
“RMSWSA” which refers to ‘Rocky Mountain Seismic Well Shoot- 
ing Association.” All other abbreviations are the same as in the 1944 
Index. 


* A limited number of reprints of this Index are available to specially interested 
members of the Society, on request to the business manager. 

+ Continental Oil Co., Ponca City, Okla., serving as special editor to assemble and 
revise this Index. 
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PATENTS* 


GEOCHEMICAL PROSPECTING 


U. S. No. 2,398,580. F. W. Crawford. Iss. 4/16/46. App. 6/17/40. Assign. Phillips 
Petroleum Co. 
Method of Prospecting for Hydrocarbons. A method of extracting hydrocarbons from 
mud or cuttings by evacuating the sample container, connecting it to an evacuated con- 
densing chamber and then heating the sample and analysing the uncondensed effluent. 


U. S. No. 2,399,965. R. Weber. Iss. 5/7/46. App. 10/25/37. and 1/4/39. Vested in 
Alien Property Custodian. 

Method for Determining Combustible Gases in Gas Mixtures. A method of deter- 
mining hydrocarbons in soil gas by removing the gas from a shallow borehole, passing 
it over a heated filament whose resistance is measured in a bridge circuit during com- 
bustion, then passing the products and residual gas over the filament again and com- 
paring the resistance values. 


U. S. No. 2,400,420. L. Horvitz. Iss. 5/14/46. App. 7/20/42. 
Geochemical Prospecting. A method of soil analysis prospecting in which samples 
are collected from a uniform depth and analysed for a selected alkaline-earth metal. 


SEISMOGRAPH PROSPECTING 
U.S. No. 2,400,245. W. H. Mayne. Iss. 5/14/46. App. 1/21/44. Assign. Olive S. Petty. 


Seismic Surveying. A circuit for refraction seismograph operations having remotely 
recording amplifiers located at the detectors and in which control of the amplifier and 
its sensitivity and telephone communication are all effected over the same two con- 
ductors which are used for recording. 


WELL LOGGING AND SURVEYING 
U.S. No. 2,398,324. B. Pontecorvo. Iss. 4/9/46. App. 8/10/43. Assign. Well Surveys, 
Inc. 


Well Surveying. A method of locating cement behind well casing by mixing with 
the cement a neutron absorbing material and subsequently running a log using a source 
of neutrons and detector of secondary gamma radiation. 


U.S. No. 2,398,562. W. L. Russell. Iss. 4/16/46. App. 12/28/44. 

A pparatus for Well Logging. A mechanical well logging and dip measuring appara- 
tus having four styles which mechanically contact the formation wall at points go° 
apart on a horizontal plane and electrically record at the surface the variations in 
penetration of each style together with the azimuthal orientation of the apparatus and 
its depth. 


* Abstracts by O. F. Ritzmann, Gulf Oil Corporation Patent Department. 
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U.S. No. 2,398,761. C. B. Aiken. Iss. 4/23/46. App. 11/30/40. Assign. Schlumberger 

Well Surveying Corp. 

Method and Apparatus for Simultaneous Determination of Various Properties of 
the Subsoil. An electric logging system using a two wire cable for simultaneously 
measuring potentials at different electrode spacings by feeding down the cable several 
frequencies, generating difference frequencies by a balanced rectifier in the hole and 
applying these to the current electrodes, picking up potentials of various frequencies 
at differential potential electrode spacings and separately indicating them by means of 
filters at the surface. 


U.S. No. 2,398,800. J. W. Millington. Iss. 4/23/46. App. 7/6/40. Assign. Sperry-Sun 

Well Surveying Co. 

Electrical Prospecting Method and A pparatus. An electrical well logging system hav- 
ing a high frequency variable amplitude oscillator inductively coupled to a formation 
and in which the oscillator is modulated by a variable frequency audio oscillator also 
coupled to a formation, the output being detected and fed to the surface as an audio 
signal whose amplitude and frequency are measured. 


U. S. No. 2,398,868. R. W. Stuart and D. Silverman. Iss. 4/23/46. App. 6/26/43. 
Assign. Stanolind Oil and Gas Co. 
Apparatus for Electrical Well Logging. A device for electrically contacting a for- 
mation behind the casing which shoots into the formation a perforating bullet trailing 
an insulated wire. 


U.S. No. 2,400,046. F. E. Hummel. Iss. 5/7/46. App. 2/22/44. Assign. Shell Develop- 
ment Co. 
Mud Logging System. A continuous mud gas logging apparatus having an inclined 
pipe through which part of the mud is drawn into a separator by a water injection 
stream, the gas being pumped from the top of the separator to a recording analyzer. 


U. S. No. 2,400,170. D. Silverman. Iss. 5/14/46. App. 8/29/42. Assign. Stanolind Oil 
and Gas Co. 

Time Cycle Te:2metering. A telemetering system for transmitting well logging or 
other signals from the bottom of a well while drilling and which has a motor whose 
speed is controlled by the signal, a commutator on the motor to periodically connect 
an oscillator between the drill stem and an insulated drill bit, the oscillator signal 
being picked up by ground electrodes at the surface and its time duration recorded. 


U.S. No. 2,400,593. J. Neufeld. Iss. 5/21/46. App. 5/7/37. Assign. Well Surveys, Inc. 

Method of and Apparatus for Investigation of Cased Drill Holes. A method of elec- 
trically logging cased holes by passing current from the casing head of one well to the 
casing head of an adjacent well and measuring the vertical potential distribution inside 


the casing. 
U. S. No. 2,400,678. G. FE. Archie. Iss. 5/21/46. App. 11/27/43. Assign. Shell De- 
velopment Co. . 


Method and Apparatus for Electrically Logging Wells. A method of logging the 
salinity of formations in an empty or oil filled hole by measuring self potentials with 
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electrolytic electrodes containing different salts or different concentrations and re- 
cording the potentials. 


U. S. No. 2,401,280. J. E. Walstrom. Iss. 5/28/46. App. 6/18/42. Assign. Schlum- 
berger Well Surveying Corp. 

Method for Locating Permeable Formations Traversed by Boreholes. A method of 
locating permeable formations from the thickness of mud cake opposite them by 
adding a magnetic material to the mud during drilling and subsequently making a 
magnetic log of the borehole wall by measuring either retained or induced magnetism 
of the mud cake. 


U.S. No. 2,401,371. J. M. Pearson and C. R. Nichols. Iss. 6/4/46. App. 5/2/35. Assign. 

Sperry-Sun Well Surveying Co. 

Electrical Prospecting Method and A pparatus. An electrical well logging system with 
three well electrodes, two of which are supplied with a constant value of commutated 
d-c, and in which both rectified and unrectified potentials are measured between the 
third electrode and one at the surface through low pass filters. 


U.S. No. 2,401,576. R. R. Mason. Iss. 6/4/46. App. 9/22/43. 


Device for Measuring Length of a Line Run into or Removed from a Bore Hole. A 
cable measuring device having several measuring wheels interconnected by a chain 
and pressure wheels mounted on adjustable springs. 


U.S. No. 2,401,893. E. B. Williams, Jr. Iss. 6/11/46. App. 5/6/43. 


Side Wall Core Barrel. A side wall coring device having jaws which grip the side 
of the hole and set a whipstock to deflect the core barrel into the side wall. 


MISCELLANEOUS 


U. S. No. 2,397,560. R. C. Olesen. Iss. 4/2/46. App. 8/18/43. Assign. Consolidated 

Engineering Corp. 

Mass Spectrometry. An automatic beam sweeping circuit for a mass spectrometer 
whose accelerating voltage has a fixed main voltage and a varying voltage obtained 
from an auxiliary voltage supply which is regulated by the voltage on a discharging 
condenser. 


U.S. No. 2,397,661. D. G. C. Hare. Iss. 4/2/46. App. 9/16/44. Assign. The Texas Co- 


Radiation Detector. A Geiger-Muller counter type of gamma ray detector in which 
a number of anode wires pass through aligned holes in the cathode surfaces and the 
chamber is filled with a mixture of argon and petroleum ether. 


U. S. No. 2,398,238. E. M. McNatt. Iss. 4/9/46. App. 5/9/42. Assign. Standard Oil 
Development Co. 
Integrating Apparatus. An apparatus for integrating the equation f(t) =/op(x—y) 
- q(x) -dx by passing a chopped light beam to a photocell through screens outlining 
the functions / and g and counting the pulses. 
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U.S. No. 2,398,249. R. G. Richardson. Iss. 4/9/46. App. 2/23/43. Assign. George W. 
Borg Corp. 
Timing Apparatus. A watch rate timer in which the ticks are amplified and strobo- 
scopically illuminate a rotating scale driven from a constant frequency source so as to 
make one and one-half revolutions per tick. 


U.S. No. 2,398,372. C. A. Green. Iss. 4/16/46. App. 1/11/44. Assign. Curtiss-Wright 

Corp. 

Device and System for Pressure Measurement. An electrical pressure gauge in which 
the pressure is applied to the annular space between two concentric deformable thim- 
bles, the inner thimble having an electric resistance type strain gauge wound on its 
inner surface and the outer thimble having a similar gauge wound on its outer surface, 
the gauges being in a bridge circuit. 


U.S. No. 2,398,375. D. Heenan. Iss. 4/16/46. App. 7/5/43. Assign. F. H. Breen. 


Liquid Level Indicator. A liquid level indicator for a tank in which the liquid 
pressure is balanced by a mercury column, the height of mercury being detected by 
an electric circuit having contacts in the manometer tube or by a resistance unit in 


the tube. 
U.S. No. 2,398,401. D. Bancroft and R. W. Goranson. Iss. 4/16/46. App. 2/17/44. 


Dynamic Pressure Gauge. An underwater explosion pressure gauge made of a hol- 
low cylinder with end diaphragms and an axial internal support on which a filament 
type strain gauge is mounted. 


U.S. No. 2,398,488. T. Zuschlag. Iss. 4/16/46. App. 4/8/41. Assign. Magnetic Analysis 

Corp. 

Magnetic Analysis. An apparatus for magnetically testing ferromagnetic material 
for flaws or dimension irregularities having a long core wound with two end coils, the 
specimen being placed near one end of the core, a dummy piece of material near the 
other end of the core, and unbalance between the coils indicated in an a-c bridge. 


U.S. No. 2,398,666. R. E. Reason. Iss. 4/16/46. App. 4/27/42 and 5/25/43. Assign. 
Taylor, Taylor & Hobson, Ltd. 


Apparatus for Measuring or Indicating the Roughness or Undulations of a Surface. 
A detector for a roughness indicator in which a flat armature carrying the stylus is 
pivoted at the center leg of an E shaped field structure, the outer legs of the E having 
variable air gaps to the armature and coils connected in an a-c bridge circuit. 


U. S. No. 2,398,701. F. A. Firestone. Iss. 4/16/46. App. 6/29/42. Assign. United 
Aircraft Corp. 
Supersonic Inspection Device. A supersonic mechanical inspection apparatus in 
which short vibration wave trains are imparted to the material and reflections from 
flaws picked up, amplified, rectified and exhibited on a c-r tube screen. 


U.S. No. 2,398,815. E. FE. Turner, Jr. Iss. 4/23/46. App. 4/1/41. Assign. Submarine 
Signal Co. 


Submarine Signaling. A mounting for a ship’s submarine signaling apparatus in 
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which the transducer may be extended into a hydraulically inflated rubber blister on 
the ship’s bottom. 


U.S. No. 2,398,816. E. E. Turner, Jr. Iss. 4/23/45. App. 10/6/41. Assign. Submarine 
Signal Co. 
Submarine Signaling. A tank mounting for a ship’s submarine signaling apparatus 
in which there is interposed between the transducer and the external water a half wave 
length system made of a plate, a portion of the tank fluid and the skin of the ship. 


U.S. No. 2,398,857. F. H. Rolfes. Iss. 4/23/46. App. 12/28/42 and 10/26/43. 


Device for Recessing Blasting Cartridges for Insertion of Detonators and of Deto- 
nator-Retaining Clips. A tool for longitudinally piercing a dynamite cartridge to make 
a recess for the cap and which also may be used to make a transverse cut for a cap 
retaining clip. 


U.S. No. 2,398,934. D. G. C. Hare. Iss. 4/23/46. App. 9/27/41. Assign. The Texas Co. 


Radiation Detector. A Geiger-Muller counter type of gamma ray detector having 
a number of parallel plates with aligned holes which carry thin walled cathode cylinders 
having axial anode wires. 


U.S. No. 2,398,987. H. Ziebolz. Iss. 4/23/46. App. 5/23/42. Assign. Askania Regulator 

Co. 

Apparatus for Continuously Predicting a Trend in Observed Data. An automatic 
machine for measuring specimens and plotting the distribution curve of a number of 
measurements and which removes the first of the number of measurements when an 
additional measurement is made. 


U.S. No. 2,398,988. H. Ziebolz. Iss. 4/23/46. App. 3/11/43. Assign. Askania Regulator 

Co. 

Apparatus for Detecting Trends in Observed Data. An automatic machine which 
measures a characteristic of a number of specimens, classifies the measurements into 
groups and indicates the number of measurements in each group, so that from the 
change in distribution curve one may observe the trend of the production operation. 


U.S. No. 2,399,014. R. B. Foster. Iss. 4/23/46. App. 2/5/43. Assign. R. H. H. Foster 
and R. FE. Foster. 


Ground Speed Meter for Aircraft. A ground speed meter in which a false terrain on a 
motor driven endless belt is synchronized both in direction and speed with the image 
of the real ground terrain as observed through a telescope. 


U.S. No. 2,399,017. S. Goldman. Iss. 4/23/46. App. 7/26/43. Assign. General Electric 

Co. 

Locating Equipment. An echo ranging equipment in which a directional continuous- 
wave transmitter and receiver are continuously rotated with an adjustable angle be- 
tween them, the azimuth of reflection being indicated on a c-r tube screen and the 
range being determined by the angle between transmitter and receiver. 


U.S. No. 2,399,018. G. R. Greenslade. Iss. 4/23/46. App. 1/19/43. Assign. Flannery 
Bolt Co. 
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Testing Device. Apparatus for electrically testing hollow metal articles for flaws, 
having an insertable probe with two contacts which engage the inner wall of the article 
and two contacts which engage an adjacent piece of similar metal on the probe, the 
contacts being in a bridge circuit which indicates resistance variations. 


U. S. No. 2,399,034. D. D. Huyett and W. E. Fletcher. Iss. 4/23/46. App. 4/16/40. 
Assign. Atlas Powder Co. 
Electric Explosion Detonator. An electric blasting cap having a thermoplastic plug 

to which is bonded the thermoplastic insulation on the leg wires, the plug being bonded 

to either a metal or plastic shell by using a solvent. 


U.S. No. 2,399,212. R. C. Dawson, R. Couch, Jr., and V. A. Prange. Iss. 4/30/46. App. 
6/20/44. Assign. Syncro-Cut Service Co. 


Disintegrator Drill, A drill for metal having an arcing electrode mechanically oscil- 
lated in synchronism with the a-c supply and a rectifier in series with the arc circuit. 


U. S. No. 2,399,213. P. G. Edwards. Iss. 4/30/46. App. 7/18/44. Assign. Bell Tele- 
phone Laboratories, Inc. 
Timing Circuit. A system for reducing leakage in condenser circuits of long time- 
constant by connecting the condenser potential to the grid of a tube whose cathode re- 
sistor output maintains a guard circuit at the same potential as the condenser. 


U.S. No. 2,399,420. H. Ziebolz. Iss. 4/30/46. App. 7/24/42. Assign. Electronbeam, Ltd. 
Electronic Translating Device. An evacuated translating device in which the force de- 
flects a pivoted beam carrying symmetrically placed tube elements, the plate circuits 
being connected in a bridge circuit whose unbalance is amplified and used to counter- 
balance the force. 


U. S. No. 2,399,513. H. B. Shapiro. Iss. 4/30/46. App. 6/27/42. Assign. The Brush 
Development Co. 


Piezoelectric Transducer. A crystal microphone whose housing has a resiliently 
mounted partition which carries both the diaphragm and the crystal. 


U.S. No. 2,399,635. L. F. Hope. Iss. 5/7/46. App. 2/28/44. Assign. General Motors 
Corp. 
Torsiograph. Apparatus for observing torsional vibrations of a shaft by mounting 
a toothed wheel a-c generator on the shaft, combining its output in a mixer circuit 
with a synchronized saw-tooth signal and observing the distorted wave on a c-r tube. 
screen. 


U. S. No. 2,399,640. C. F. Kettering. Iss. 5/7/46. App. 4/13/43. Assign. General 

Motors Corp. 

Tem perature Sensitive Means. An infra-red detector having a mirror to focus the 
rays on a target made of four insulated blackened silver quadrants which are inter- 
connected by thermocouples so that vertical and horizontal temperature gradients 
may be measured separately. 
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U.S. No. 2,399,671. E. G. Gage. Iss. 5/7/46. App. 4/30/43. Assign. Two-thirds to 
Leon Ottinger. 
Tridimensional Radio Navigational System. A system for locating an airplane from 
the ground by using a transmitter on the plane and three or more direction finding 
receivers on the ground which indicate on two cathode ray tubes. 


U.S. No. 2,400,097. N. Brewer. Iss. 5/14/46. App. 4/18/44. Assign. Fischer & Porter 
Co. - 
Flowmeter. A rotameter type flowmeter in which the float has a female thread en- 
gaging a long pitch thread on a central guide rod. 


U.S. No. 2,400,102. W. M. Cobb, U. S. A. Iss. 5/14/46. App. 4/16/41. 


Detonator or Blasting Cap. An electric blasting cap of spherical shape and having 
a base charge detonated by heat of thermite which is ignited by nitrocellylose. 


U.S. No. 2,400,108. S. Elowson. Iss. 5/14/46. App. 9/4/42. Assign. Pittsburgh Equi- 
table Meter Co. 
Wide Range Flow Meter. A flow meter in which a piston is raised in a tapered tube 
and when the limit of travel is reached the tube is raised to uncover a variable orifice 
and thus extend the scale. 


U.S. No. 2,400,190. R. W. Clark. Iss. 5/14/46. App. 1/12/44. Assign. Radio Corp. of 
America. 


Multirange Meter. A multirange meter with a relay in the circuit which when actu- 
ated introduces a meter multiplier and mechanically indicates the multiplier factor. 


U. S. No. 2,400,305. A. Hazeltine. Iss. 5/14/46. App. 6/13/44. Assign. The Gray 

Manufacturing Co. 

Electromagnetic Sound Recorder and Reproducer. A phonograph record cutter or re- 
producer which is mounted on flexible springs for reproducing and clamped for record- 
ing. 

U. S. No. 2,400,308. E. C. Klepp. Iss. 5/14/46. App. 3/25/42 and 6/2/43. Assign. S. 

Smith & Sons (England), Ltd. 


Electric Ratiometer. An indicating ratiometer having a field magnet and moving 
system carrying two coils whose current ratio is measured by the equilibrium position 
in a variable field, the field variation being obtained by variation in the axial length 
of the pole piece. 


U.S. No. 2,400,326. A. Wolf. Iss. 5/14/46. App. 9/30/43. Assign. Geophysical Research 

Corp. 

Electrical Measuring Instrument. A circuit for measuring the ratio of two a-c volt- 
ages by changing the frequency of one, then passing both through a self-controlled 
amplifier, separating them with filters, rectifying the altered frequency, feeding it to 
the control element of the amplifier and measuring the unaltered frequency output. 


U.S. Re. 22,759 (Orig. No. 2,368,831). C. M. Hathaway. Iss. 5/21/46. App. 
9/19/42 and 1/10/46. 
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Oscillograph Galvanometer. A bifilar oil immersed oscillograph galvanometer having 
a permanent magnet and element cell with integral pole pieces and an oil expansion 
diaphragm. 


U. S. No. 2,400,552. H. Hoover, Jr. Iss. 5/21/46. App. 4/11/42. Assign. Consolidated 
Engineering Corp. 


Range Finder. A sound direction finder in which the difference in arrival times at 
the receivers is converted to potentials applied to plates of a c-r tube so that the result- 
ing spot deflection indicates the resultant direction of the source. 


U.S. No. 2,400,557. R. C. Lawlor. Iss. 5/21/46. App. 7/31/42. Assign. Consolidated 
Engineering Corp. 
Analytical System. A method of analysing mixtures with a mass spectrometer in 
which reference samples and the mixture are each ionized at various ionization voltages 
and the ion intensities compared. 


U.S. No. 2,400,564. B. F. McNamee. Iss. 5/21/46. App. 2/26/43. Assign. Consolidated 
Engineering Corp. 
Pickup. A transverse vibration pickup for airplane propellers in which the armature 
is mounted on leaf springs parallel to the propeller radius, one spring extending inward 
and one outward so that rotation does not affect the pickup characteristics. 


U. S. No. 2,400,571. R. C. Olesen. Iss. 5/21/46. App. 2/26/43. Assign. Consolidated 
Engineering Corp. 


Electrical Measurement. A modulated carrier type of strain recording system in 
which zero shift of the galvanometer when changing from velocity to acceleration 
recording is compensated by introducing an opposing battery current in the latter 
channel. 


U.S. No. 2,400,870. R. E. Bates. Iss. 5/28/46. App. 12/7/42. Assign. Submarine Signal 
Co. 


Submarine Signaling Apparatus. An under-water signaling device for ships which 
has a supersonic transducer inside a streamlined housing which may be swung into the 
water from a sea chest and inside of which the transducer may be turned into any 
azimuth. 


U.S. No. 2,400,923. W. H. Farr and W. F. Fagen. Iss. 5/28/46. App. 3/19/43. Assign. 
Stewart-Warner Corp. 
Gas Detection Apparatus. A carbon monoxide detector in which the gas is contacted 
with a heated filament whose change in resistance is observed with a Wheatstone 
bridge and in which the space around the filament is loosely packed with inert material, 


U.S. No. 2,400,940. H. J. De N. McCollum. Iss. 5/28/46. App. 3/22/43. 


Gas Detection Apparatus. A thermionic diode having a heated cathode of platinum 
tubing through which the gas being tested is passed, presence of hydrogen being indi- 
cated by an increase in diode current which operates a relay. 
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U. S. No. 2,401,019. F. Rieber. Iss. 5/28/46. App. 2/28/42. Assign. Interval Instru- 
ments, Inc. 

Time Interval Meter. A time interval meter in which charging of a condenser is ini- 
ated at the beginning and stopped at the end of the interval and the ratio of condenser 
potential to its charging source potential measured by a potentiometer with a vacuum 
tube null indicator. 


U. S. No. 2,401,049. B. J. Campbell and M. Hilton. Iss. 5/28/46. App. 12/13/43. 
Assign. Consolidated Vultee Aircraft Corp. 


Method of Making Electrical Strain Gauges. A method and jig for making a fila- 
ment type strain gauge which permits measuring and adjusting its resistance and sub- 
sequently cementing it onto cover strips in a clamping jig. 


U.S. No. 2,401,094. M. G. Nicholson, Jr. Iss. 5/28/46. App. 6/23/44. Assign. Colonial 

Radio Corp. 

Time Delay Apparatus. A time delay device having a Jong magneto-strictive rod 
with a magneto-striction pulse-generator at one end and a pick-up an adjustable dis- 
tance away, the delay being determined by the transmission time of the pulse along the 
rod. 

U. S. No. 2,401,140. C. H. Cordie. Iss. 5/28/46. App. 4/24/40. Assign. Hercules 

Powder Co. 

Explosive Cartridge Assembly. A device for assembling dynamite cartridges end to 
end and having a metal sleeve with resilient inward pointing teeth which can bite into 
the cartridges to keep them from coming apart. 


U. S. No. 2,401,163. R. Keller. Iss. 5/28/46. App. 10/29/42. and 12/3/43. Assign. 

Aktiengesellschaft Brown, Boveri & Cie. 

Device for Measuring Mechanical Acceleration by Electrical Means. A system for 
measuring or controlling acceleration of rotating machines by mounting an a-c gener- 
ator on the shaft, connecting it to parallel capacity and inductance whose currents 
operate a differential indicator or controller relay. 


U.S. No. 2,401,250. C. W. Kandle. Iss. 5/28/46. App. 9/30/43. 


Earth Drill, An earth auger which has a hinged trap over the auger opening and a 
reversible cutting blade to produce different depths of cut. 


U. S. No. 2,401,328. R. Black, Jr. Iss. 6/4/46. App. 1/16/43. Assign. Bell Telephone 
Laboratories, Inc. 


Sound Translating Device. A directional microphone having one side of the dia- 
phragm exposed to the atmosphere and the other side coupled to an acoustic phase 
shifting network having a tube with a series of path lengths to the atmosphere. 


U.S. No. 2,401,404. A. V. Bedford. Iss. 6/4/46. App. 1/12/44. Assign. Radio Corp. of 
America. 


Electrical Multiplying Circuit. A circuit for multiplying a speech signal S and a cod- 
ing signal K by generating the signals —K, —S, (S+K), (-S—K), (—S+K) and 
(S—K), squaring each one and combining the squared values to form the product SK. 
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U.S. No. 2,401,410. W. T. Brown. Iss. 6/4/46. App. 2/10/44. 


Detection of Sulphur in Coal. A method of detecting sulphur impurities in coal by 
exposing it to ultra-violet light and observing the fluorescence of impurities. 


U.S. No. 2,401,411. S. F. Carlisle, Jr., and A. B. Mundel. Iss. 6/4/46. App. 7/30/42 

Assign. Sonotone Corp. 

Frequency Response Testing System. An amplifier testing system for obtaining the 
response curve on a c-r tube screen in which the sweep oscillator also varies the nor- 
mally fixed frequency of a beat-frequency oscillator between desired limits through 
control of a variable reactance tube. 


U.S. No. 2,401,432. D. G. C. Luck. Iss. 6/4/46. App. 9/29/43. Assign. Radio Corp. of 
America. 
Radio Position Locator. A radio pulse-echo position locator in which the radial 
deflection of the c-r tube spot is produced by two deflecting coils operating alter- 
nately and rotated in azimuth in synchronism wth the transmitting antenna. 


U. S. No. 2,401,467. F. B. Doane. Iss. 6/4/46. App. 5/22/39. Assign. Magnaflux 

Corp. 

Method of and Apparatus for Magnetic Testing. A magnetic method of testing steel 
for flaws by establishing two magnetic fields whose relationship may be adjusted or 
varied continuously so that the resultant field will cross the flaw at an optimum angle 
and detecting flaws with paramagnetic particles. 


U.S. No. 2,401,527. A. W. Vance. Iss. 6/4/46. App. 3/30/43. Assign. Radio Corp. of 

America. 

Electromechanical Multiplying Device. An electrical signal multiplying device in 
which a voice coil drives one plate of a condenser microphone, one signal being applied 
as a-c to the voice coil and one signal being applied as d-c to the microphone, and the 
product being taken off as a-c output. 


U. S. No. 2,401,565. P. J. Holmes. Iss. 6/4/46. App. 5/28/41. Assign. One-half to 

R. R. Stoddart. 

Radio Navigation System. A radio navigation system for airplanes in which signals 
from two crossed loop antennae are fed into balanced modulators and combined with 
signals from a non-directional antenna so that the bearing to the source may be indi- 
cated on a phase indicator. 


U. S. No. 2,401 641. R. L. Hill. Iss. 6/4/46. App. 8/26/42. Assign. Atlas Powder Co. 


Explosive Assembly. A dynamite cartridge having on its wall a rigid plate with two 
channels into which caps or fuse may be placed. 


U. S. No. 2,401,712. F. Turrettini. Iss. 6/4/46. App. 7/22/41 and 6/5/42. Assign. 
Societe Genevoise d’Instruments de Physique. 


Device for Determining Electrically the Position of Objects. A magnifying scale 
reading device in which an optical glass oscillates an image of the scale past a similar 


| 


PATENTS 557 


grid and periodically excites a photocell whose amplified impulses flash a neon lamp 
illuminating an index rotating synchronously with the image oscillations. 


U. S. No. 2,401,723. J. H. Deming. Iss. 6/11/46. App. 2/12/43. Assign. Texaco De- 
velopment Corp. 
Method and Apparatus for Locating Objects. A method of making land mines loca- 
able by burying with them a small amount of radioactive material which may after- 
ward be found with a Geiger counter. 


U.S. No. 2,401,747. H. J. Dibblee. Iss. 6/11/46. App. 4/21/44. 

Timing Device. A method of measuring elapsed time between two events by charg- 
ing a condenser, initiating its discharge at the first event, interrupting the discharge at 
the second event and measuring the residual voltage. 


U.S. No. 2,401,876. A. L. Marker and E. W. Bergere. Iss. 6/11/46. App. 5/23/44. 

Tensiometer. A cable tensiometer having two blocks placed against one side of the 
cable and an arm on the opposite side which kinks the cable, the arm being rotated by 
a common torque wrench. 


U.S. No. 2,401, 929. J. H. Hammond, Jr. Iss. 6/11/46. App. 5/4/42. Assign. Radio 
Corp of America. ; 
Detection of Submerged Objects. A system for detecting submarines by towing a 


series of sources of water-penetrating rays kept at some depth by paravanes and a 
series of receivers at a shallower depth and observing when shadows or reflections occur. 


U. S. No. 2,402,235. L. A. Burrows and W. E. Lawson. Iss. 6/18/46. App. 4/6/40. 
Assign. E. I. du Pont de Nemours & Co. 
Blasting Initiator. A blasting cap in which the primary detonating compound con- 
tains finely divided water insoluble metal soap. 


U. S. No. 2,402,321. P. J. S. Enger and A. Parschin. Iss. 6/18/46. App. 6/16/43. 

Assign. Telefonaktiebolaget L. M. Ericsson. 

Device for the Transformation of Mathematical Functions. A device for producing 
mechanical motion or electrical signals following a symmetrical function of a driving 
motion in which the curve is reproduced on only one side of the point of symmetry and 
relays used to change the sign. 


U. S. No. 2,402,392. R. Goldschmidt. Iss. 6/18/46. App. 10/6/42. Assign. Radio 

Patents Corp. 

Electromagnetic Sound-Transmission Apparatus. A bone-conduction hearing aid 
transducer having an armature fastened to the skull and a permanent magnet field 
with wound pole pieces, a tuned reed being used as a magnetic shunt to provide fre 
quency correction. 


U. S. No. 2,402,403. C. M. Hines. Iss. 6/18/46. App. 1/10/42 and 5/27/42. Assign. 
Westinghouse Air Brake Co. 
Electric Decelerometer. An automatic railway car braking device having generators 
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on the car axles connected so that unbalance due to wheel slipping operates relays to 
control the applied airbrake pressure. 


U.S. No. 2,402,459. R. B. Smith. Iss. 6/18/46. App. 2/3/43. 


Submarine and Ship Detector Using Reflected Radio Waves. A system for automatic- 
ally locating ships or airplanes by sending out an ultra short wave radio signal, receiv- 
ing any reflection and thereupon sending a coded radio signal to a distant central 
station. 


U.S. No. 2,402,464. H. Suter. Iss. 6/18/46. App. 7/12/41. Assign. One-fourth to W. 
Ockrant. 


Velocity Measuring System. A method of measuring velocity by radiating an ultra 
high frequency electromagnetic weve, receiving the wave reflected from the moving 
body and observing the beat frequency between radiated and reflected waves. 


U.S. No. 2,402,531. D. R. Christian. Iss. 6/25/46. App. 6/4/42 and 3/28/44. Assign. 
The Brush Development Co. ; 


Transducer. An underwater crystal microphone of small size in which the dia- 
phragms are sealed by an electroplated metal coating. 


U. S. No. 2,402,666. A. Raspet. Iss. 6/25/46. App. 3/27/42. 


Helical Spring. A helical spring having a pitch angle for which the rotary magnifi- 
cation expressed as rotation of the spring per unit elongation is inversely proportional 
to the length of the spring. 
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“Radiation of a Sound in Water as Affected by Depth of Submersion,” by N. N. 
Andreyev, L. M. Brekhovskikh, and L. D. Rosenberg. (Vol. XLVII, No. 6, 
1945.) 
“Propagation of Sound Waves in a Liquid Layer Between Two Absorbing Half- 
Spaces,” by L. M. Brekhovskikh. (Vol. XLVIII, No. 6, 1945.) 
“The Sinevka Salt Dome and Its Geological Place,” by D. N. Sobolev and M. N. 
Shkabara. (Vol. XLVIII, No. 6, 1945.) 


PUBLICATIONS RECEIVED 559 


Economic Geology, Vol. XLI, No. 3 (May), No. 4 (June-July—in 2 parts), 1946. 
“Mineral Resources of China,” by V. C. Juan (June—July—part 2). 

Geofisica Pura e Applicata, Vol. VIII, Fasc. 3-4, Fasc. 5-6, 1946. 

Journal of Applied Physics, Vol. 17, No. 6 (June), No. 7 (July), 1946. 

Journal of the Institute of Petroleum, Vol. 32, No. 268 (April), No. 269 (May), No. 270 

(June), 1946. 

Nafta (“Petroleum”), Petroleum Institute of Poland, Krakow. Rok II, Nr 5 (May), 
and Rok II, Nr 6 (June), 1946. 

Review of Scientific Instruments, Vol. 17, No. 5 (May), No. 6 (June), No. 7 (July), 1946. 

Revista de Ciencias, Organo de la Facultad de Ciencias de la Universidad Nacional 
Mayor de San Marcos. Ano. XLVII, Numero 454, Diciembre de 1945. 

Terrestrial Magnetism and Atmospheric Electricity, Vol. 51, No. 2 (June 1946). 

World Petroleum, Vol. 17, No. 5 (May), No. 6 (June), and No. 7 (July), 1946. 


Reprints: 

Edouard Poldini, “Sur l’existence de courants électriques naturels liés aux gites 
d’anthracite valaisans.”” Compte Rendu des séances de la Société de physique 
et Vhistoire naturelle de Genéve. Vol. 60, No. 3, pp. 274-284 (August-Decem- 
ber 1943) 

Edouard Poldini, ‘“‘Contribution a ]’étude, par prospection électrique, des terrains 
anisotropes plissés.’”? Compte Rendu des séances de la Société de physique et 
histoire naturelle de Genéve. Vol. 60, No. 1, pp. 6-13 (January-March 
1943). 

Edouard Poldini, ““Quelques considérations sur l’interprétation géophysique.” Bull. 
Soc. Vaud. Sc. Nat. Vol. 62, No. 262, 1944, pp. 403-410. 

J. F. Enslin, D.Sc. ““A Comparison of the Gravimeter and Torsion Balance Methods 
of Geophysical Prospecting.’’ The Transactions of the Geological Society of 
S. Africa. Vol. XLVIT, 1944, pp. 283-299. 

B. D. Maree, M.Sc. “The Vredefort Structure as Revealed by a Gravimetric 
Survey.” The Transactions of the Geological Society of S. Africa. Vol. XLVTI, 
1944, pp. 183-196. 

B. D. Maree, M.Sc. ‘‘Notes on Some Magnetometric Surveys with Special Ref- 
erence to Anomalies on Karroo Dolerites.” The Transactions of the Geological 
Society of S. Africa. Vol. XLVI, 1943, pp. 107-118. 
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SPECIAL REVIEWS COMMITTEE 


Tuomas A. ELxins, Gulf Research & Development Co., Pittsburgh, Pa., Chairman 
R. A. GEYER, Humble Oil and Refining Co., Houston, Texas 
NorMAN RICKER, Carter Oil Co., Tulsa, Oklahoma 
D. C. SKEEts, Standard Oil Co. of New Jersey, New York, N. Y. 


Emanometric oil and gas prospecting,* by Sylvain J. Pirson, Petroleum Engineering, 

vol. 17, NO. 4, pp. 132, 134, 136, 138, 140, 142 (January 1946). 

Those geophysicists not specializing in geochemical prospecting who have watched 
with a great deal of interest the vicissitudes of the only method seriously claiming to 
discover oil directly will be interested in this frank account of difficulties and problems. 

Pirson’s method is based on the measurement of the rate of escape of ethane at 
the surface. That measurable amounts of a hydrocarbon or other gas can escape from 
an oil or gas deposit is a fundamental postulate of the method, and various objections 
that have been raised are discussed. The new name “emanometric” now replaces 
“geodynamic,” but merely indicates increased precision in measurement and new 
techniques in interpretation. 

These new interpretational techniques include (1) the recognition of a uniform 
hydrocarbon diffusion “background” field, explained by Pirson as the result of present- 
day generation of oil and gas; (2) the removal of artificial anomalies caused by altitude 
variations in the surface topography or by ground water movement; (3) the correction 
of readings for climatic, seasonal, and atmospheric effects using a system of repeated 
base stations; (4) the correction for deep geologic effects, such as shifts in the position 
of highs caused by the tilting of sedimentary beds; (5) the correction for the nature of 
the soil. This latter correction is difficult and should be eliminated by the proper selec- 
tion of station sites. Even this is far from easy and Pirson says: “The field operations 
are so complex, it appears that a geologist trained in soil science and who has a thorough 
understanding of the emanometric prospecting problem should be selected to conduct 
field operations.” 

An interpretation problem not yet solved is the separation of a deep leakage source 
from an overlying shallow one. Here the possibility of an electrolytic model, based on 
the analogy between d’Arcy’s law and Ohm’s law, is suggested. 

_ Since no data are presented and no indications are given as to the success or failure 
of the new method in practice, it must be assumed that it is still under experimental 


investigation. 
Tuomas A. ELKINS 


Ueber die Struktur und Figur der Erde (The structure and figure of the earth), by W. 
Heiskanen, Gerlands Beitraege zur Geophysik, vol. 57, pp. 132-170 (1941). 


Heiskanen’s paper constitutes a valuable contribution to the literature of geodesy, 
geology, and tectono-physics. It contains detailed summaries and analyses of the major 


* Some, but not all, of the material in this paper is summarized in Pirson’s recent 
review of nine years’ work in this field, ‘Disturbing factors in geochemical prospecting,” 
GEOPHYSICS, vol. 9, no. 3, pp. 312-320 (July 1946). 
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theories in these fields on the basis of information available up to 1940. The section 
dealing with the rapid anomaly changes along oceanic coast lines, marine gravity anoma- 
lies and similar anomalies on islands such as the Philippines, Hawaii, Cypress and others 
are of special interest because of the current activity in the geophysical exploration 
of continental shelf zones. The remainder of the paper deals with the evaluation of vari- 
ous gravity formulae, sources of error in gravity anomalies, convection current hy- 
potheses with special reference to their application as a mechanism for explaining moun- 
tain building cycles, theories dealing with down buckling of the earth’s crust and origin 
of major negative gravity anomaly zones. 

Heiskanen is in accord with Vening Meinesz’s down buckling theory and believes 
that it explains the origin of the negative anomaly zones of the East and West Indies 
and probably the one situated to the east of the Japanese Island Arc. He does not 
believe that any isostatic reduction can explain them and in addition the inhomo- 
geneities of mass responsible for these zones can not be situated in the near surface 
layers of the earth’s crust since the zones do not conform to topography. Hence, their 
origin must be sought for at much greater depths. The results of Kuenen’s laboratory 
experiments in which he used a mixture of paraffin, vasoline and oil for his ‘earth’s 
crust” are also in agreement with this theory. Heiskanen contends that Byerly’s results 
are also in harmony with the down buckling theory because Byerly was able to demon- 
strate from seismic data that the Sierra Nevada Mountains in California actually pos- 
sess a root. The width of this root varies from 40 to 70 kilometers and is situated at a 
depth of 30 kilometers (calculated from the floor of the normal earth’s crust). It appears 
from one statement that he makes on this point that he has misinterpreted to some ex- 
tent Byerly’s original paper. Heiskanen states that the existence of the root was dis- 
covered because the first of the earthquake waves whose focus was in northern or south- 
ern California was not received in the Owen’s Valley station. Actually the waves arrived 
at this station but the arrival time indicated that they had travelled at a much slower 
velocity indicative of a path through granite rather than basalt. 

The large positive isostatic gravity anomalies ranging from +30 to +100 milligals 
found in the oceanic areas can not be eliminated by any isostatic reduction, because 
even after applying all the possible reductions, the anomalies are still positive, It is 
therefore necessary to assume that an excess of mass exists in these regions. Heiskanen 
believes that the almost invariable rapid increase of isostatic gravity anomalies along 
the oceanic coasts is even more significant than the existence of the large positive 
anomalies found over the open ocean. The anomalies along the coast are generally 
negative and over the oceans near the coast are positive or less negative. The changes in 
anomalies after regional corrections have been applied may be as much as 78 milli- 
gals. The Hayford anomalies exhibit even a greater change. They indicate that this 
phenomenon can not be eliminated by merely assuming different compensation values. 
In attempting to arrive at some explanation for this phenomenon he discusses in some 
detail work by Bullard, Daly, Ewing, and Vening Meinesz. 

A detailed and critical evaluation of Grigg’s convection current theory leads Heis- 
kanen to the conclusion that if convection currents are rising below the continents and 
descending beneath the oceanic and coastal seas, it is possible to offer an explanation 
on this basis for the excess of gravity over the deep coastal areas and the rapid rise of 
gravity along the oceanic coasts, as well as the difference between continental and 


oceanic gravity values. 
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In his analysis of different gravity formulae he also discusses the validity of the 
tri-axial ellipsoid formula. The gravity characteristics of North America, Europe, the 
Mediterranean, India, Siberia, and a portion of the Indian and the Pacific Oceans, 
lend validity to the assumption of a tri-axial gravity formula. However, the large posi- 
tive gravity anomalies in the Gulf of Mexico, the Caribbean Sea and the negative anom- 
alies in East Africa, contradict the validity of the existence of the tri-axial formula. 
In addition, the gravity anomalies in the southern hemisphere do not seem to fit par- 
ticularly well into the tri-axial formula. If all, or nearly all of the latitude zones for the 
larger equatorial axis show almost similar directions, then the tri-axiality of the earth 
is most probably a fact. However, if the directions, \9, deviate appreciably from one 
another, then one must doubt the existence of tri-axiality. 

Heiskanen in his section on sources of error in gravity anomalies discusses in some 
detail the following factors: 

. The accuracy of the elevation determination of the station 

. The particular gravity formula used 

. The isostatic system used (Hayford, Airy-Heiskanen, Meinesz) 

. The assumed density of the earth’s crust 

. The accuracy of the available topographic and bathometric charts 

The accuracy of estimating elevations and depressions on the charts 

. The accuracy of the reduction tables used 

. The accuracy of the isostatic reduction charts 

. The fact that the indirect effect of the undulations of the geoid upon the gravity 
anomalies is taken into consideration 


He also points out that all American gravity anomalies must be increased by 6 milligals 
because of a new adjusted value in the tie between Potsdam and Washington, D. C. 
base stations. Care must also be used when including older gravity observations in 
any calculations because of the large differences in values observed when old stations 
have been reoccupied. For example, 47 gravity observations made between 1936 and 
1937 in Poland and Hungary differ on an average by 18 milligals from measurements 
made at the same stations between 1891 to 1901. Differences as large as 39, 33, and 
30 milligals were observed when stations were reoccupied at Bidermannsdorf, Szenic, 
and Lwow, respectively. 

Heiskanen concludes his paper with a plea for closer cooperation between geode- 
sists, geologists and geophysicists in order to avoid many mistakes and to reach more 
rapidly the desired goals of all. It is in order in concluding this review to call attention 
to Heiskanen’s “Catalogue of the Isostatically Reduced Gravity Stations” appear- 
ing in Annales Academiae Scientiarum Fennicae Ser. A. Tom. LI, N:0 10. This catalogue 
published in 1939 should be of particular value to anyone investigating gravity anoma- 
lies and corrections because it contains 4000 isostatically reduced gravity stations in- 
cluding several hundred marine gravity observations. 


CMY AN 


R. A. GEYER 


Recent work in the field of high pressures, by P. W. Bridgman. Reviews of Modern 
Physics, vol. 18, no. 1, pp. 1-93 (January 1946). 


This authoritative summary of work in the field of high pressures will be welcomed 
by geophysicists, since a knowledge of the physical properties of substances under the 
pressure of deep burial in the earth is of basic importance in their work. The summary 
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covers all work which has been published in accessible form since the publication of 
Professor Bridgman’s book The Physics of High Pressure (1930) up to the month of 
June 1945. The scope of the summary is indicated by the following section headings: 
Technique. 
1. Extension of pressure range. 
2. Measurement of pressure. 
3. Miscellanecus technique. 
Mechanical Effects of High Pressure. 
1. Volumetric effects in gases. 
. Volumetric effects in liquids. 
. Volumetric effects in solids. 
. Phase changes under pressure in one-component systems. 
. Phase changes under pressure in multi-component systems. 
. Effects of pressure on viscosity. 
. Effects of pressure on elastic constants. 
. Effect of pressure on plastic flow and similar phenomena. 
9. Miscellaneous mechanical effects of hydrostatic pressure. 
Thermal Effects of High Pressure. 
1. The effect of pressure on specific heats. 
2. Miscellaneous thermal effects of high pressure. 
Electrical Effects of High Pressure. 
1. The effect of pressure on the electrical resistance of solids. 
2. The effect of pressure on the electrical resistance of solutions. 
3. The effect of pressure on dielectric constant. 
4. Miscellaneous electrical effects of pressure. 
Magnetic Effects of High Pressure. 
Optical Effects of High Pressure. 
1. The effect of pressure on index of refraction. 
2. Photographic effects of high pressure. 
3. Miscellaneous optical effects of high pressure. 
Chemical Effects of High Pressure. 
Biological Effects of High Pressure. 
In evidence of the great amount of work which has been done in the field of high 
pressures during the past fifteen years, the bibliography lists 674 referneces. 
NorMAN RICKER 


on 


IRWIN ROMAN 


CONTRIBUTORS 


Irwin Roman received the degrees of 
A.B. from Washington University in 1913, 
M.A. from The University of Chicagoin 1915, 
and Ph.D. from The University of Chicago in 
1920. From 1925 to 1929 he was employed by 
the Geophysical Research Corporation. He 
has taught mathematics and physics at Wash- 
ington University, Northwestern University, 
Vanderbilt University, and Michigan College 
of Mining and Technology. He has been with 
the U. S. Department of the Interior since 
1935, in the Geological Survey and Bureau of 
Mines. His technical contributions to geo- 
physics have included various publications, a 
well explorer and a magnetic gradiometer. He 
is a fellow of the American Association for the 
Advancement of Science, a member of the 
Mathematical Association of America and 
Vice-President of the Section of Terrestrial 


Magnetism and Atmospheric Electricity in the American Geophysical Union. He is 
also a member of Sigma Xi and Phi Beta Kappa. 


D. H. WEINSTEIN received the degree of 
A.B. in mathematics from the University of 
California in 1930. He received the degree of 
Ph.D. in physics from the California Insti- 
tute of Technology in 1933. He has been em- 
ployed by The Superior Oil Company since 
1936. From 1942-1945 he did developmental 
work at the Signal Corps Laboratory, Fort 
Monmouth, N. J. 


D. H. WEINSTEIN 
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F. J. WrLttams received the degree of 
B.S. in Electrical Engineering from the Uni- 
versity of Oklahoma in 1926. He was a stu- 
dent engineer with Commonwealth Edison 
Company, Chicago, Illinois, from 1926 to 
1927, and engaged in engineering work from 
1927 to 1936. Since 1936, he has been engaged 
in seismograph exploration work with Geo- 
physical Service Inc., Dallas, Texas, except 
during 1941 and 1942, when he was employed 
as an engineer by Dallas Power and Light 
Company, Dallas, Texas. He is a member of 
the Society of Exploration Geophysicists. 


F. J. WILLIAMS 


Photographs and biographies of contributors to this issue not appearing above 
have been published previously, as follows: C. Hewitt Drx, Vol. VI, No. 4, p. 460 
(October 1941). THomas S. WEst, Vol. IX, No. 4, p. 573 (October 1944). CLARENCE 
C. BeacuaM, Vol. IX, No. 4, p. 574 (October 1944). G. W. Horton, Vol. VIII, No. 3, 
p. 318 (July 1943). ALFRED WOLF, Vol. VII, No. 4, p. 428 (October 1942). B. G. Swan, 
Vol. VII, No. 4, p. 425 (October 1942). 


: 
4 
4 Bee 
A 


THE SOCIETY ROUND TABLE 


APPLICATIONS FOR MEMBERSHIP APPROVED FOR PUBLICATION 


The Executive Committee has appro ved for publication the names of the following 
candidates for membership in the Society. This publication does not constitute an 
election, but places the names before the membership at large. If any member has 
information bearing on the qualifications of these nominees, he should send it to the 
Secretary within thirty days. (Names of sponsors are placed beneath the name of each 
nominee.) 

"ACTIVE 

Bastiaan Baars 

F. Goldstone, L. K. Mower, L. L. Nettleton 
Lloyal Orrin Bacon 

H. R. Joesting, S. J. Pirson, J. H. Swartz 
Dr. Perry Byerly 

J. B. Macelwane, V. J. Blum, Florence Robertson 
Francis Robert Coryn 

L. A. Scholl, Jr., R. L. Lay, C. R. Wallace 
Raymond Wesley Dudley 

A. R. Shugart, J. D. Perryman, W. W. Newton 
William Jesse Edmonston 

W. W. Newton, H. M. Houghton, Clyde L. Hayes 
Richard A. Geyer 

W. M. Rust, Jr. D. H. Gardner, D. C. Skeels 
Robert Howard Glover 

Frank Ittner, C. E. Reel, L. B. Trombla 

Paul Greenlee 

T. L. Allen, W. L. Crawford, P. E. Narvarte 
Walton Henry Hohag, Jr. 

M. J. Walczak, Harry L. Thomsen, Nevin R. Shade 
D. Mitchell Jennings 

Joseph A. Sharpe, G. A. Foster, A. F. Cline 
Byron K. Johnson 

H. B. Peacock, K. E. Burg, A. E. McKay 
Martin Victor Linwell 

Cecil H. Green, Floyd J. Williams, Martin C. Kelsey 
Hugh McCain, Jr. 

C. J. Harvey, M. J. Walczak, E. C. Livingston . 
Robert Rhea Phillips 

H. B. Peacock, E. D. Gaby, N. J. Smith 
Smiley Raborn, Jr. 

T. J. Bevan, Tom D. Mayes, Jack Peters 
Victor Pierre Ulrich 

L. L. Nettleton, F. Goldstone, L. K. Mower 
Klaas van der Weg 

N. D. Smith, F. A. van Melle, K. R. Weatherburn 


ASSOCIATE 
Myron Donald Adams 
W. W. Newton, H. M. Houghton, J. L. Morris 
(Continued on page 571) 
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1947 ANNUAL MEETING 


The 1947 annual meeting of the Society of Exploration Geophysicists will be held 
in Los Angeles, California, in conjunction with the annual meetings of the A.A.P.G. 
and S.E.P.M. Tentative arrangements have been made for headquarters at the Bilt- 
more Hotel during the period March 24 to 27. 

It is certainly not too early for ideas to crystalize into definite outlines for papers. 
We have all heard the comment in past years by prospective authors that they might 
have been able to turn out papers had the idea occurred to them earlier in the season. 
Then there are those particular subjects which require extended amounts of time in 
order to insure adequate treatment in a paper. Therefore, it is important that commit- 
ments for technical papers be made to the regional program committees as soon as pos- 
sible. The names of members serving on all the regional committees were announced in 
the July 1946 issue (Volume XI, No. 3). 

As announced in the July issue,‘ special Editor’s Committee is at work in an effort 
to obtain a sufficient number of new Case History type papers to complete Volume I 
of a projected Case Histories series in time for printing by the end of the current year. 
If possible the initial sale of this first volume will be made at the annual meeting. It is 
the feeling of the Executive Committee that new Case History type papers contributed 
to this volume could also be used in part at least for program material at the annual 
meeting. 

Prospective authors are referred to L. L. Nettleton’s paper entitled “Preparation 
of Manuscript for Publication in Gzopuysics” in Volume X, No. 3 of GEopuysics, for 
worthwhile suggestions on preparation of manuscript and illustrative material. 

Ceci, H. GREEN 
General chairman 
Program Committee for 1947 meeting. 


ANNOUNCEMENT OF PACIFIC COAST MEETING 


A local meeting of the Society of Exploration Geophysicists will be held in connec- 
tion with meetings of the American Association of Petroleum Geologists and the So- 
ciety of Economic Paleontologists and Mineralogists in Los Angeles on November 7 
and 8. 
The S. E. G. Technical Session will be held on the morning of November 7. The 
tentative program includes a symposium on Multiple Reflections led by T. P. Ells- 
worth, Richfield Oil Corp., with prepared short discussions and illustrative slides 
from about ten West Coast geophysicists; a general discussion of the Bikini tests by 
R. B. Morgan, Jr.; “Seismograph Correlation Grading,”’ by Phil Gaby, Consultant to 
Standard Oil Company of California. Other possible material for this meeting is not 
yet fully developed. 

Committee members for the fall meeting are: C. H. Dresbach, Western Gulf Oil 
Co., Chairman, Milton Born, Amerada Pet. Corp.; E. H. Vallat, Ohio Oil Co.; and 
Joseph C. Waterman Shell Oil Co. 


‘ 


PERSONAL ITEMS 


Lr. Cox. Ropert W. Jackson, Post Signal Officer, has been honored by award of 
the Legion of Merit for his contribution to victory in World War II. The medal was 
presented by Major General Clift Andrus, Commandant of the Field Artillery School, 
in a ceremony held recently in the Commandant’s office. 


Lt. Col. Robert W. Jackson, Post Signal Officer, Ft. Sill, Okla., was presented the 
Legion of Merit by Major General Clift Andrus, Commandant of the Field Artillery 
School. (U. S. Army Photo) 
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Colonel Jackson received the Legion of Merit for “outstanding services as Assistant 
Signal Officer, Luzon Force, on Bataan, Philippine Islands, from 12 March to 9 April 
1942. Despite determined enemy efforts to locate his installations, his ingenuity and 
prompt and energetic actions kept them constantly operating both in his own and 
lower echelons.” 

Captured by the Japanese in April 1942, Colonel Jackson lived through the Bataan 
Death March, after which he remained a prisoner of war until freed in September 1945. 
The early part of his imprisonment was spent in Camps O’Donnell and Cabanatuan 
in the Philippines, but he was later removed to Japan where he was finally liberated. He 
wears the American Defense ribbon with one battle star, the Asiatic-Pacific ribbon with 
star, Philippine Defense with star and the Victory ribbon. 

Prior to entering military service, Colonel Jackson was a geophysicist, employed 
by the Geophysical Service Inc. of Dallas, Texas. During five and a half years with this 
company he took part in domestic operations as well as in such foreign assignments as 
Canada, Colombia, and Sumatra. He attended North Texas State Teachers College 
and is a member of the Society of Exploration Geophysicists. 


R. G. Exms, formerly Party Chief with Geophysical Service, Inc., has taken the 
position of Geophysicist for The Union Sulphur Company of Sulphur, La. 


NorMaN A. HASKELL has accepted the position of Geophysicist with the United 
States Smelting Refining and Mining Company. His mailing address is now in care of 
the Geophysical Laboratory of this company at 1126 Boylston Street, Boston, Massa- 
chusetts. 


Rosert J. WatsON, formerly Geophysicist with the Carter Oil Company, Tulsa, 
is now Exploration Manager for the Standard Oil Company (Bahamas), Ltd. in Nassau. 


Witson G. SAVILLE, senior partner in Gravity Meter Exploration Company, 
Houston, and former United States Army District Engineer at Galveston, was sworn 
in as Chairman of the Houston-Harris County Port Commission this July after his 
appointment by the Houston City Council and Harris County Commissioners Court. 


R. F. ALDREDGE, Geophysicist with Westby Geophysical Corporation, has moved 
to Tampico, Tamps., Mexico where his address will be c/o Petroleos Mexicanos, De- 
partmento Exploracion. 


WarING BrRabD Ley, formerly with Geochemical Surveys, Dallas, has established 
offices as an independent geologist in Stamford, Texas. His mailing address is now Box 


28, Stamford, Texas. 


A. B. Hair has resigned his position as District Supervisor for National Geo- 
physical Co., Inc., of Dallas to organize Marine Exploration Company under his indi- 
vidual ownership. The new company will specialize in marine seismic surveys. His 
address is now 3200 West Dallas, Houston 6, Texas. 


ALEXANDER WOLF, formerly Director of the geophysical laboratory of The Texas 
Company enjoyed a month’s vacation while completing plans for opening an office in 
Houston as consulting physicist. 


Cuar_Es C. Bares returned from participating in the preliminary oceanographic 
and geologic survey of Bikini and neighboring atolls on July 1 to accept the position of 
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Oceanographer with the newly formed Division of Oceanography, Hydrographic Cffice, 
U. S. Navy. 

LEON FiscuER, Seismologist with Shell Oil Company, has transferred to Lamar, 
Colorado, Box 352. 


V. E. Cutip has moved to 687 31st Street, Oakland 9, California. 


Francis F. CAMPBELL has moved to Tulsa, Oklahoma, where his new mailing 
address will be Box 2040, Tulsa 2, Oklahoma. 


Davin H. Scort, of The Texas Company, has moved to Santa Barbara, where his 
office is located at 320 Benjamin Franklin Building. 


Dart WANTLAND has resigned his position with Seismograph Service Corporation 
to accept a position as Geophysicist for the U. S. Reclamation Bureau in Denver, 
Colorado. 

D. G. EGan is now with the Geophysical Department of International Petroleum 
Company, Apartado 1081, Lima, Peru. 


James R. HEArp, recently elected to associate membership, holds the position of 
Computer with United Geophysical Company and at the present time is located at 
2412 31st Street, Lubbock, Texas. 


PRIESTLY LEVERICH, former Party Chief of Mayes-Bevan Co., has accepted the 
position of Gravity Computer with W. E. Brogan Company. He can be reached in care 
of the company at 2711 Crocker, Houston 6, Texas. 


B. F. Owrncs is now Party Chief for Garrett Exploration Company, Box 267, 
Jacksboro, Texas. 


Cart L. Bryan transferred on July 1 from Gulf Research and Development Com- 
pany to Gulf Refining Company, Shreveport, La. His address is now P. O. Drawer 
1731, Shreveport, Louisiana. 


A. V. Hares, who recently advanced from Associate to Active membership, is 
now living in Birthright, Texas. 


Putt P. Gasy resigned his position with Geophysical Service, Inc. on May 1 of 
this year to take up new duties as consulting geophysicist for the Standard Oi] Company 
of California and subsidiaries. His office will be in the Standard Oil Building, San Fran- 
cisco, California. 


J. R. Maxey is now Senior Seismologist with the United Geophysical Company, 
Apartado 1085, Caracas, Venezuela. 


Upon his release from the United States Army, Louis C. PAKISER returned to his 
former position as geophysicist with Carter Oil Company, P. O. Box 801, Tulsa, 
Oklahoma. 

J. L. KEzE.ER resigned his position with Seismograph Service Corporation to be- 
come President of Key Velocities, Inc. of Tulsa, Okla. Mr. Kezeler formed the new 
company for the sole purpose of making velocity measurements for the seismograph 
companies. 
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GEORGE E. Waconer, for the past nine years Assistant Chief Geophysicist in 
charge of Geophysical Exploration for The Carter Oil Company, and the present 
Secretary-Treasurer of Society of Exploration Geophysicists, has taken up new duties 
as Division Exploration Manager for The Carter Oil Company in Shreveport, La. His 
mailing address is Drawer 1739, Shreveport, Louisiana. 


V. L. Jones has been appointed Professor of Geophysics at The University of 
Tulsa. After several years of industrial and consulting geophysics, Professor Jones has 
rejoined the University Staff and will head the newly created Department of Geo- 
physics. He will also continue his consulting practice in geophysics. 


ALEXANDER WOLF has established offices as Consulting Physicist at 4708 Caroline 
Street, Houston 4, Texas. 


Lt. Col. M. C. Bowsky has received his commission in the Regular Army, and is 
on active duty as Deputy Administrator, Office of the Zone Administrator, War Assets 
Administration, Bannister Road and Troost Ave., Kansas City 4, Missouri. 


P. E. NARVARTE, engaged for the past 12 years by the Petty Geophysical Engineer- 
ing Company in the interpretation of seismic data, has resigned the position of District 
Supervisor with that company to do independent geophysical consulting. His temporary 
address is 1422 West Hollywood Ave., San Antonio 1, Texas. 
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work, transformers sold by Harrison 
are hermetically sealed against sub- 
tropical humidity . . . give improved 
uniformity, minimum size and weight. 
Also available with plug in terminal 
block. Special-design transformers can 
be built to meet your demands if stan- 
dard sizes do not fit. 


GEOPHONE CABLE 


Harrison recently introduced a new, 
exclusive design multiple conductor 
geophone cable that meets the most 
exacting requirements for dependable 
performance. It eliminates cross-feeds 

. each conductor is insulated with 
real rubber latex color .coded .. . it 
contains no moisture absorbing mater- 
ials. It can be supplied from present 
stocks in 1, 8 and 13 pairs. 


A COMPLETE LINE OF GEOPHYSICAL 
AND ELECTRONIC SUPPLIES . 


ELECTRONIC - OIL EXPLORATION - SUPPLIES 


1422 SAN JACINTO 


HOUSTON, TEXAS 
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An A.A.P.G. Publication! 
Of The 


UNITED STATES 


Prepared under the Direction of the Committee on Tectonics, 
Division of Geology and Geography, National Research Council. 


CHESTER R. LONGWELL, Chairman, PHILIP B, KING, Vice-Chairman 
CHARLES H. BEHRE, WALTER H. BUCHER, EUGENE CALLAGHAN, D. F. 
HEWETT, G. MARSHALL KAY, ELEANORA B. KNOPF, A. I. LEVORSEN, 
: T. S. LOVERING, GEORGE R. MANSFIELD, WATSON H. MONROE, J. T. 
PARDEE, RALPH D. REED, GEORGE W. STOSE, W. T. THOM, JR., A. C. 
WATERS, ELDRED D. WILSON, A. O. WOODFORD 


A New Geologic Map of the United States and Adjacent 
Parts of Canada and Mexico 


7 Geologic structure, as evidenced and interpreted by a combination of out- 
cropping areas, bedrock, surface disturbance, and subsurface deformation, 
is indicated by colors, symbols, contours, and descriptive explanation. 
Igneous, metamorphic, and selected areas of sedimentary rock are mapped. 
Salt domes, crypto-volcanic disturbances, and submarine contours are shown. 


The base map shows state boundaries, rivers, a pattern of cities, and 1- 
degree lines of latitude and longitude. 


The scale is 1:2,500,000, or 1 inch = 40 miles. Printed in 7 colors, on 2 
sheets, each about 40 x 50 inches. Full map size is about 80 x 50 inches. 


PRICE, POSTPAID 
$2.00 rolled in mailing tube 
; $1.75 folded in manila envelope 
$1.50 in lots of 25, or more, rolled or folded 


The American Association of Petroleum Geologists 
Box 979, Tulsa 1, Oklahoma, U.S.A. 
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“a complete and well-rounded geophysical service. 
rating in all parts of the 
ice is available for surveysinany partof 
erica. Western service is available for surve) 
| United States and in South America. 
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| THE Correction 


GEOTECHNICAL | of Address 


: ! CO RPO R ATION Please note below the corrected mail- 


ing address for the new business 
office of the Society, located at 213 


ROLAND F. BEERS 
&g Ritz Building, Tulsa, Oklahoma. All 


President 
correspondence other than editorial 
matters should be addressed as fol- 

1702 Tower Petroleum Building lows, to 
Dallas, Texas C. C. Campbell, Business Manager 
| Telephone L D 101 SOCIETY OF EXPLORATION 
| 7 GEOPHYSICISTS 
Research Laboratory P. O. Box 1614 
Cambridge, Massachusetts Tulsa, Oklahoma 


C. H. FROST GRAVIMETRIC 
SURVEYS, INC. 


C. H. Frost, President Josera A. Suanpe, Vice-President 


GRAVIMETERS manufactured under license from Standard Oil 
Development Company 


GRAVIMETRIC AND MAGNETIC SURVEYS carefully con- 


ducted by competent personnel 


GEOLOGIC INTERPRETATION of the results of gravimetric 


and magnetic surveys 


1242 South Boston Avenue Tulsa 3, Oklahoma 
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Cyclotron Specialties Instruments are readily assembled into a 


of radioactivity. They combine simplicity and economy with all 
important requirements of the most exacting investigators 
and have won recognition of the leading workers in this field 
throughout the world...in educational institutions, industrial 
plants, laboratories and hospitals. 


For new workers, the use of a cathode ray oscillograph, giving 
an excellent visual pattern of the pulses, is recommended. The 
counter set is fitted with amplifier and external connections for 
this purpose. The Cyclotron 401-A Impulse Register is also 
widely used by itself for high speed impulse recording or me- 
chanical operation counting. - 


file Yh 

©, Le 


he 


lay 


yclotron Specialties Instruments are quickly assembled into 
complete radioactivity apparatus... 

@ sample to be tes’. ~ . placed in an ashing dish where its 
adioactivity is transmitted into the Geiger Tube through the 
t's thin mica window. Impulses from the tube are ampli- 
d and extinguished by a First Stage. This passes the im- 
pulses on to the Counter Set where they are amplified, scaled 
a scale of eight and indicated on the Impulse Register. 


IMPULSE REGISTER 

Accurately registers up to sixty im- 
pulses per second, and will solve your 
problems in high speed impulse or 


mechanical oper- 
ation counting. Hh 
Price, FOB factory 


GEIGER MULLER COUNTER 


An especially fine piece of equipment with $ : 
numerous applications in radioactivity 


work. Price, FOB factory... 


THIN MICA WINDOW GEIGER TUBE 
A compact, completely enclosed counter 


tube of sturdy construction with many im- $ 
portant features. Used extensively by lead- 
ing laboratories. Price, FOB factory =» « », 


Cyclotron Specialties radioactivity instruments are rug- 
ged and practically foolproof...made from the highest 
grade materials and desighed for dependable perform- 
ance‘under actual working conditions. Place your order 
now by mail. request 


Detailed specifications upon 


Cyclotron Specialties Com 


Moraga 8, California 
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PERFORMANCE 


SEISMOGRAPH 
RECORDING 
PAPER 


and humidity 


Because it has been especially designed to com- 
bine photographic excellence with resistance to 


abuse, Haloid Record delivers consistently fine 
recordings, even under the most adverse conditions 
of field and laboratory. It provides sharp contrast 
and legibility . . . it is plenty tough . . . it resists 
heat and withstands humidity. Wherever critical 
geophysicists demand superior seismograph re- 
cordings, you will find Haloid Record a popular 


choice. 


THE HALOID COMPANY 
465 Haloid St. Rochester 3, N. Y. 
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SIX 
DUAL 
AMPLIFIE! 


Heiland 24 Trace Seismic Equipment is the result of years of 
operating experience under all field conditions of weather, 
temperature, humidity, and terrain. Careful research and 
modern production facilities assure a precision product, sturdy, 
shock-mounted, with all major components quickly removable 
by means of plug and receptacle connections for inspection 
and maintenance. For exceptional accuracy, dependability, 
and operational ease, use Heiland Seismic Equipment. 


COMMUNICATION COMPANDER MIXER FORK AND FORK AMPLIFIER 


SWITCH 
} PANEL 


SPARE AMPLIFIER 


Please mention GEopHysiIcs when answering advertisers 


25 


DUAL 


H 
| 
arueys for You 
»_ On 
RECORDER 
| e e a ] C ) 
: 


; GEOPHYSICS the Journal of the Society of Exploration Geophysicists 


SEISMIC SURVEYS 


KEYSTONE EXPLORATION COMPANY 


OFFICES AND LABORATORY 
2813 WESTHEIMER ROAD 


HOUSTON @ TEXAS 
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The North American Gravity Meter is highly 
favored by both crews and operators for the many 
advantages it offers in gravity surveying. 


Crews prefer it for its extreme light weight and 
compactness. It reaches the zenith of portability 
in that one man can carry it on the back pack. 
It can be mounted in the back seat of standard 
sedan or small "jeep" and readings made in any 
location accessible by car in less than 2 minutes 
simply by extending tripod through car floor. 


This portability enables crews to cover consider- 
ably more area per day with the North American 
Gravity Meter. 


Its high degree of accuracy has been demon- 
strated many times by unbiased operators. In one 
instance a group of base stations checked within 
.01 of a milligal of the original values established 
a year earlier by another North American Gravity 


Meter. 
Many of these meters have never required servic- 
ing or even after ng of field opera- 
tion... a feature extremely valuable to operations 
le. in foreign fields or remote locations. 
NORTH AMERICAN GEOPHYSICAL CO. 
Gravity — Magnetic — Seismic Surveys Geophysical Apparatus 


636 Bankers Mortgage Bldg., Charter 4-3523 
Houston 2, Texas 
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STATES EXPLORATION 
COMPANY 


Gravity Meter Surveys 
H. L. Schiflett, Manager 


Service Interpretation 


P.O. Box 845 


Sherman, Texas 


GEOPHYSICAL 
SURVEYS 


EXPLORATION 
COMPANY 


2044 Richmond Rd. 
HOUSTON 6, TEXAS 


Paul Charrin 
President 
C. C. Hinson 
Vice President 


John Gilmore 
Vice President 


STANDARD PROCEDURE 


for Thermador Geophysical Transformers 


At Thermador there’s nothing unusual about 
measuring a transformer at 1 MV. For all Ther- 
mador input, interstage and low-level trans- 
formers are measured at 1 MV and 100 MV. 
Furthermore, this information is placed on 
each transformer, assuring uniformity in any 
quantity ordered. 


TOLERANCES HELD 

TO A REMARKABLE MINIMUM 

Thermador transformers are held to a maxi- 

mum tolerance of plus 5% to minus 5%. Spe- 

cial transformers with tolerances held from 

plus 3% to minus 3% are also available. 

SPECIAL STAFF... 

EXCLUSIVE TESTING EQUIPMENT 

A special geophysical staff and exclusive testing 

equipment assure the superiority of Thermador 

transformers. All inductance measurements are 

made with a specially designed Owens Bridge 

by applying certain measured voltages across 

the coil. 

HERMETICALLY SEALED 

Thermador Geophysical Transformers are her- 

meticaily sealed and vacuum impregnated to 

provide resistance against extreme conditions 

of humidity, heat and cold. 

es complete information, write Dept. GT to- 
ay. 


"Seven Leagues Ahead”’ 


THERMADOR 


ELECTRICAL MANUFACTURING CO. 
5119 District Boulevard « Los Angeles 22, Californic 
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During the past 19 years we have been privileged to con- 
duct Barret Magnetic Surveys for 73 clients, representing 
a good cross section of the exploration branch of the pe- 
troleum industry. 


It is a significant fact that most of our work in recent years 
has been done for some of our earliest clients. 


These clients and others continue to use our services for 
the sole reason that it has proved profitable to them in 
their search for new oil reserves. 


We will be pleased to explain without obligation how you 
may use Barret Magnetic Surveys to cut your exploration 
costs in terms of barrels of oil discovered. 


Ane 
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ADVANCED 


EXPLORATION 
COMPANY 


GEOPHYSICAL SURVEYING 


Spearpoink of the Petroleum Industry 
For ADVANCED 


Seismic Equipment and 
Technique 


CALL ADVECO F-8007 


622 FIRST NATIONAL BANK BUILDING 
HOUSTON 2, TEXAS 


Cc. W. BOCOCK, Ill GEO. D. MITCHELL, JR. 


JAMES L. SAULS, JR. 
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KING “FRONT-END’ WINCH 


for Civilian Willys Jeep 


The Model 100 King Winch for Civilian Willys Jeep incorporates the 
same dependable features that have been used for years on passenger cars 
and half ton pick-ups in the oil fields all over the world. These rugged 
and dependable features together with the latest improvements, such as 
centralized lubrication and improved drive-shaft suspension, assure a 
“front-end” winch that will give trouble-free service and long life. 


Power for the winch is taken from the front end of the engine crankshaft 
by means of a solid sheave and a patented sliding clutch using rubber con- 
tact blocks to absorb shock and misalignment. This clutch assures a positive 
drive, eliminates slippage, and can be engaged or disengaged at any time 
(even under load and with engine idling). 


The winch sets directly in the center on front of the Jeep and is easily and 
safely operated by one man. Recommended cable 150’ 5/16” 6 x 19 hemp 
center wire rope. 


Weight of complete installation 126 lbs. Speed ratio cable drum to engine 
72 to 1. 


Sold Exclusively Through Willys Dist:*butors and Dealers 


KOENIG IRON WORKS 


2214 Washington Ave. HOUSTON 10, TEXAS 
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Announcing— 
Geophysical Instrument Company's 


New and Radically Improved 


FIELD MAGNETOMETER 


When Geophysical Instrument Company set out to build its Schmidt 
type of magnetometer, it resolved to build the best and most advanced 
type of instrument possible with modern knowledge and modern tech- 
niques. It is believed that this aim has been attained and that the new 
G-I-C magnetometer is the best instrument of its type ever offered for 
sale. Among its features are: 


INSTRUMENT HEAD 


Improved moving system; better temperature compensation arrange- 
ment; better location of knife edge. 


Improved knife edge—more resistant to injury. 

Improved arresting system—virtually impossible to injure knife edge. 
Improved thermal insulation. 7 

Adjustment of moving system with telescope in place. 


TRIPOD 
All metal light-weight construction. 
Positive lock for tripod legs. 
Totally enclosed leveling screws. 
Slow-motion drive for auxiliary magnets. 


These are only a few of the improvements. There are others too numer- 
ous to detail here. In fact, there is hardly a feature of construction which 
has not been improved. It will pay you to learn more about this out- 
standing instrument. Further details gladly furnished. 


Geophysical Instrument Company 
Arlington, Va. 


Please mention GropHysics when answering advertisers 


34 
| 
a 


GEOPHYSICS the Journal of the Society of Exploration Geophysicists 


35 


SPIRALOK* 


the time-saving 
cartridge 
assembly 


special explosives 
for seismic work 


VIBROGELS™ 


VIBROCAP™ 


the no-lag electric 
blasting cop 


*REG. U.S. PAT. OFF. BY 
HERCULES POWDER CO. 


HERCULES 


Explosives Department HERCULES POWDER COMPANY 917 King St., Wilmington, Del. 
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/ CORRECT DATA IS 
VALUABLE INFORMATION. 


INCORRECT DATA 
WORSE THAN 


© Accuracy 


Foreign aad Domestic 
.  Geo-Magnetic 


Field Surveys 


4 MAGNETEX GEOPHYSICAL 


BURK BURNETT “BLDG. — FORT WORTH 2, 
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The true yardstick for measuring the value of any seismograph explora- 


tion service is accuracy of interpretation . .. either condemning or 
approving acreage. To those operators planning an exploration campaign 
we issue an invitation to investigate our record of operations, which 


provides unquestionable proof of the accuracy of our interpretations. 


SUBSURFACE SURVEYS 


Seismic Explorations. Inc., Gulf Building, Houston, Texas 
Established 1932 
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New Type Balanced Crown Block 
Easily Removable, All Steel Tabl 
ir Actuated Double Drum Draw 


VERYONE KNOWS that the 
_ way to get the ‘low down” on 


CHECK THE HIGH STAND- 
ARDS of all —built 
equipment with any driller who 
has operated it. He'll tell you 
there is no finer or more de 
pendable equipment anywhere 


‘DURING THE WAR DAYS, 


ment that would meet the advanced needs 
of the industry when peace came. Out of 
these dreams has evolved: 


BE SURE To INVESTIGATE hi nev 
drill before making any Coppapetegetts for 


Patented Floating Gear rea 

Hydraulic Throttles 

Mud Pump Unitized with Drill 

Engines Compounded through Fluid Drive 


EQUIPMENT 
MIDLAND, TEX, 


& drilling equipment is to hunt | 
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TYPE SGO-I 
RECORDING OSCILLOGRAPH 


A small and rugged, yet wholly complete, recording oscillograph designed pri- 
marily for portable and semi-portable use without sacrifice of properties of most 


larger units. 


@ Can be supplied with 12, 18 or 24 Traces 


© Utilizes rugged pencil type mirror Galvanometers for ease 
of removal or adjustment 


@ Timing Motor driven directly by self-contained 100 cps 
Tuning Fork 


@ Removeable Magazine for simplicity of operation 


@ Storage Magazine holds 150’ of 6” or any narrower width 
paper 


Small size, 8 9/6” x 14” x 121/.” deep. Total weight is 41 pounds. Other units 
to choose from if desired. 


Prices sent on Request 


SOUTHWESTERN INDUSTRIAL ELECTRONIC CO. 


Custom Built Geophysical Equipment 
2619 Milam Houston 6, Texas 
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TYPE GRU-2 
SEISMOGRAPH AMPLIFIER UNIT 


Complete self-contained unit—requires only 6-volt battery, geophones, and 
oscillograph to make a complete recording unit. Embodys the following features: 
@ Geophone continuity checking circuit 
@ Line balancing controls for “High Line” pickup elimination 
@ Paralleling switch for testing 
@ 35 cycle Oscillator for matching and testing 
@ Communication circuit for communication with shot point 
@ Electronic mixing (if desired) 
@ Automatic control circuit 
@ Battery compartment and built-in meter to check voltages 
@ 12 Geophysical amplifiers with all the latest features 
% Full automatic volume control 
%& Variable band pass filtter 
% Initial suppression control 
%& Hermetically sealed transformers 


Entire unit is rugged and dependable. Ideal size for truck mounting—43” x 22” 
x 16” deep. Other units to choose from if desired. 


Prices on Request 
SOUTHWESTERN INDUSTRIAL ELECTRONIC CO. 


Custom Built Geophysical Equipment 
2619 Milam | Houston 6, Texas 
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In oil exploration, as in a theater, a 
little light on the subject is desirable. 
Authentic data on which to base your 
decisions with no blind guessing 
raises your chance of success. 


The Mayes-Bevan organization make accurate gravimetric surveys and 
soundly interpret them for you, so, why be uncertain, let us help you with 
your oil exploration problems. 


KENNEDY BLDG TULSA, OKLAHOMA 
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@ WHY RISK the possibility of overlooking 
important producing horizons? Many new 
productive sands have been tested due to 
the results of Baroid Well Logging Service. 


@® KNOW YOUR WILDCAT! To obtain the 
most comprehensive data from your well, 
correlate the mud analysis log with the 
electric log and other well information. 


@ BAROID WELL LOGGING SERVICE supplies 
@ continuous analysis of the returning drill- 
ing mud and does not interfere with drilling. 


. @ Operators use it who desire immediate 
information on the occurrence of oil and 
gas as their wells are being drilled. 

@ Include BAROID WELL LOGGING SERVICE 
in your well program. 


Baroid Well Logging Units contain complete 

rd q D logging and mud testing equipment. Logging 

reveals the presence of oil and gas in the 

formations penetrated and practically elimi- 

cll nates the chance of passing any productive 
OGGUMG zones. By using this service to control a cor- 


ing program, only the necessary cores need 


be taken, resulting in a great saving in time 
& Pe RVI C ie and money in drilling a well. Mud tests are 
made to maintain a satisfactory drilling mud. 


BAROID SALES DIVISION 
NATIONAL LEAD COMPANY @ BAROID SALES OFFICES: LOS ANGELES 12, TULSA 3, HOUSTON 2 
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G. E.C. announces a NEW SERVICE—Aero- 
magnetic Exploration. For this service, G. E. C. offers the most advanced 
airborne detection equipment, independently engineered and built by 
Heiland Research Corporation. This equipment gives accurate, syn- 
chronized records of: 


Radar Triangulation (Ratran) + Radio Altitude (Terrain Clearance) 
Terrain Position * Time * Magnetic Anomalies + Barometric Altitude 


For exploration estimates on areas of 10,000 
square miles or over, write to 
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the world. 


‘This experience plus the policy of Pac impr 
through sustained research and development p 
rporation’s clients with the most effective: 
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